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Deutsche Kurzfassung 



Das theoretische Studium des Quark-Gluon Plasmas gewinnt immer mehr 
an Bedeutung seit Teilchenbeschleuniger wie das SPS, der RHIC, oder der 
sich im Bail befindliche LHC die erforderlichen hohen Energiedichten in 
Schwerionenkollisionen erreichen, die es erlauben, diesen neuen Materiezu- 
stand experimentell zu untersuchen. Einfache Anwendungen der Quanten- 
feldtheorie im Rahmen einer storungstheoretischen Entwicklung nach der 
Kopplungskonstante versagen bei hohen Temperaturen, und trotz eifriger 
Bemiihungen, die Situation in den Griff zu bekommen, haben wir quantita- 
tive theoretische Aussagen iiber den Phasentibergang nur von Gittersimula- 
tionen. Diese wiederum versagen fur ein Quark-Gluon Plasma bei hohem 
chemischem Potential und niedrigen Temperaturen, wie man es im Kern von 
dichten Sternen vermutet. Large- iV^-QCD - das ist Quantenchromodynamik 
(QCD) mit einer grofien Zahl von Quark-Sorten (number of quark flavors - 
Nf) - erlaubt es, Wechselwirkungseffekte von thermodynamischen Grofien 
wie dem thermischen Druck oder der Entropie exakt in der effektiven Kop- 
plung g^ s oc g 2 Nf fur alle Temperaturen T und chemische Potentiale \i q zu 
berechnen. Dies macht Large Nf QCD zu einem idealen Testwerkzeug fur 
verschiedene Naherungsmethoden. 

In der vorliegenden Arbeit prasentieren wir das exakte Large- Nf Resultat 
fiir den thermischen Wechselwirkungsdruck in der kompletten T-fi q -Ehene 
in einem Bereich, in dem der Einfluss durch den Landau-Pol numerisch ver- 
nachlassigt werden kann. Fiir kleine Werte der Kopplung vergleichen wir 
unser Resultat mit existierenden storungstheoretischen Ergebnissen in der 
Literatur, einschliefilich der aktuellen Berechnung des Drucks durch Vuori- 
nen fiir endliche Temperatur und chemisches Potential sowie einer alteren 
Rechnung von Freeman und McLerran fiir verschwindende Temperatur und 
hohes chemisches Potential. Unsere numerische Genauigkeit erlaubt uns, ex- 
istierende storungstheoretische Koemzienten zu verifizieren und zum Teil sog- 
ar zu verbessern, und auch storungstheoretische Koemzienten zur sechsten 
Ordnung in der Kopplung numerisch zu bestimmen, die analytisch bislang 
noch nicht berechnet wurden. Fiir verschwindendes chemisches Potential 
berechnen wir lineare und nicht-lineare Quarkzahl-Suszeptibilitaten. Wir 
zeigen, dass das moderate Skalierungsverhalten, das durch die Quarkzahl- 
Suszeptibilitaten nahegelegt wird, ziemlich abrupt bei fj, q > ttT zusammen- 
bricht, aber dass dieser nicht-pertubative Effekt in \i q immer noch in guter 
Naherung durch die Ergebnisse von Vuorinen bei kleinen Kopplungen und 
endlichem T beschrieben wird. Nur fiir T <C fi q versagt auch dieser Zugang, 
und wir kommen in den Bereich der sogenannten Non-Fermi-Fliissigkeit, 
die im Gegensatz zur klassischen Fermi-Fliissigkeit von langreichweitigen, 
quasistatischen transversalen Eichbosonen dominiert wird. In diesem Limes 
konnen wir nicht nur den bereits bekannten fiihrenden TlnT -1 Beitrag zur 
speziflschen Warme vervollstandigen, sondern auch iiber die fiihrende Ord- 
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nung eine storungstheoretische Reihe mit anomalen gebrochene Potenzen 
rp{3+2n)/3^ durch dynamische Abschirmung verursacht werden, angeben. 
Wir berechnen deren Koeffizienten analytisch bis zur Ordnung T 7//3 und find- 
en, dass diese tatsachlich das fiihrende anomale Verhalten der vollen QED 
und QCD bestimmen (also bei endlichem Nf). 
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Abstract 

The theoretical study of the quark gluon plasma gains increasing interest 
as particle accelerators like the SPS, RHIC, or the currently built LHC will 
reach sufficiently high energy densities in heavy ion collisions that allow us to 
probe this new state of matter experimentally. Straightforward application of 
quantum field theory at high temperatures fails in a perturbative expansion 
in the coupling constant, and despite some effort during the last decades to 
improve the situation, so far quantitative theoretical knowledge about the 
phase transition merely comes from lattice simulations. Lattice simulations 
on the other hand fail for a deconfined quark-gluon plasma at large quark 
chemical potential and small temperatures which is expected to be found in 
the core of dense stars. Large Nf QCD - that is quantum chromodynamics 
with large number of quark flavors - allows one to calculate thermodynamic 
properties like the interaction contribution of thermal pressure or entropy 
exactly in the effective coupling g^ s oc g 2 Nf for all temperatures T and 
chemical potentials [x q . This makes large Nf QCD an ideal testing ground 
for various approximation methods. 

In this work we present the exact large Nf NLO calculation of the ther- 
mal interaction pressure in the whole T-// g -plane where the presence of the 
Landau pole is negligible numerically. For small values of the coupling we 
compare our results to existing perturbative results in the literature, in par- 
ticular the recent calculation by Vuorinen for finite temperature and chemical 
potential or an older calculation by Freedman and McLerran for zero tem- 
perature and high chemical potential. Our numerical accuracy allows us to 
verify and even improve some of the existing perturbative coefficients, and 
to predict new coefficients to the sixth order in the coupling numerically 
that have not been calculated analytically yet. For larger couplings we de- 
termine where perturbation theory ceases to be applicable. At zero chemical 
potential we calculate linear and non-linear quark number susceptibilities. 
We show that the moderate scaling behavior suggested by the quark num- 
ber susceptibilities breaks down rather abruptly at fj, q > irT, but that this 
non-perturbative effect in \i q can still be reproduced well by the calculation 
by Vuorinen for small couplings and finite T. Only for T <C ^ q also this 
approach breaks down and we enter the range of a so-called non-Fermi liq- 
uid, which in contrast to a classical Fermi liquid is dominated by long-range 
quasi-static transverse gauge-boson interactions. In this limit, we complete 
the previously known leading TlnT -1 contribution to the specific heat, and 
also go beyond this order to find a series involving anomalous fractional pow- 
ers J 1 ( 3 + 2n )/ 3 caused by dynamical screening. We calculate their coefficients 
analytically up to order T 7//3 and find that these contributions indeed de- 
termine the leading anomalous contribution in full QED and QCD (i.e. at 
finite Nf). 
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1 xin ji chi bu liao re dou fu: Those of impatient heart can not eat hot tofu (bean curd). 
(traditional Chinese saying) 
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Prolog 



LHC/CERN in Swiss, 2007 2 

SWOOOOSH....!!! 

"Oh my god! What was that?" - "They are trying to smash us!" - "Oh 
mum! I wanna get out!". A bunch of hysterically screaming lead ions flying 
almost at the speed of light can barely fight the 8.36 Tesla magnetic fields 
(some 100,000 times the strength of Earth's magnetic field) that keep them 
on spiraling orbits in the Large Hadron Collider (LHC). 

SWOOOOSH....!!! 

"Aahhh!" - "That was VERY close!!! I can't believe they are doing this 
to us!". With only 27 km collider circumference, the ion bunches approach 
themselves some 11,000 times per second. 

SWOOOOSH...!!! 

"Ok, that's it! Let us oooouut!!" - "They already ripped off all electrons 
from me! What else do they wanna do?". The darkness in the collider is tan- 
talizing, as are the 1.9 K (some -270°C) coldness from the 5,000 surrounding 
super-conducting magnet coils. 

SWOOOOSH...!!! 

"I'm cold!!!" - "Brrrr... I'm freezing!!!". "You'll feel warmer soon", Old 
uncle Joe tried to soothe them. In his billions of years of life he had gone 
through many phases. So far he had survived everything. 

SWOOOOSH...!!! 

"What the..!" - "They are not... no!! Don't tell me they want to smash 
us into..." - "..into the Quark Gluon Plasma!!!" 
SWOOO - KABOOM!!!! 

Let me not go into too much detail - it's a horrible view: atoms are 
smashed head on head, their inner protons and neutrons are torn apart, 
even the formerly strictly confined up- and down-quarks end in a disgusting 
soup of unconfined quarks and gluons, whining and wincing only as shock 
fronts penetrate them. But this is science: Scientists have to take the loss of 
some lead atoms for the sake of knowledge about the Quark Gluon Plasma 
(QGP). 

2 Data in this prolog are taken from the LHC design performance website 2003 EQ. 
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Kapitel 1 

Introduction 



1.1 The experiment 

The study of the quark-gluon plasma in heavy ion collision experiments 
is one of the focal points of contemporary high energy physics. Histori- 
cally, scientists smashed elementary particles since the 1930's with Earnest 
Lawrence's invention of the cyclotron, which improved earlier attempts with 
linear accelerators tremendously. In a cyclotron, charged particles from elec- 
tron or ion sources are accelerated on a circular orbit before they hit a target. 
Even higher energies can be reached with synchrotrons, first constructed at 
the General Electric Research Laboratory for the University of California, 
Berkeley, in 1949. Contrary to the continuous beam of a cyclotron, the 
synchrotron works with beam pulses which ride on electro-magnetic waves. 

Modern particle accelerators like the Super Proton Synchrotron (SPS) 
at CERN or the Relativistic Heavy Ion Collider (RHIC) at Brookhaven Na- 
tional Laboratory are based on this concept of the synchrotron, relying on 
several steps of particle acceleration. For example, heavy ions at RHIC are 
produced in the Tandem Van de Graaff by static electricity. These particles 
are then carried to the Booster synchrotron where they are accelerated to 
37% the speed of light, before they are transfered to the Alternating Gra- 
dient Synchrotron where they reach 99.7% the speed of light. From there 
they will be injected into the first RHIC ring in clockwise direction or into 
the second RHIC ring in counter-clockwise direction, where particles will be 
collided into one another in one of the six interaction points of the ring [2]. 

It is interesting to note that for proton-anti-proton collisions like in the 
LHC one only needs one collider ring instead of two, as positively charged 
protons and negatively charged anti-protons can travel clockwise and counter- 
clockwise in the same magnetic field. However, the LHC also has two col- 
lider rings to also allow for lead-lead collisions, proton-lead collisions, proton- 
oxygen collisions, or collisions of protons with other light ions jS]. A summary 
of accelerators that are producing or will produce the quark-gluon plasma 
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Name 


Location 


Year 


Energy 


Beam 


Orbit 


SPS 


CERN 


1994 


160 GeV 


lead ions 


6.9 km 


RHIC 


Brookhaven 


2000 


200 GeV 


gold ions 


3.8 km 


LHC 


CERN 


2007 


7 TeV 


proton, lead 


26.7 km 



Tabelle 1.1: Particle accelerators that produce or will produce the quark 
gluon plasma (QGP). The year marks the date of starting operation for 
search of the QGP, center-of-mass energy is given per nucleon for SPS and 
RHIC. The abbreviations of the names mean Super Proton Synchrotron 
(SPS) |E], Relativistic Heavy Ion Collider (RHIC) (2], and Large Hadron 
Collider (LHC) [T]. 



is given in table 11.11 Besides hadron colliders (using ions or protons) which 
will probe the quark-gluon plasma, future collider plans also include lepton 
colliders as CERN's proposed electron linear accelerator CLIC (Compact 
Linear Collider), or muon colliders. Since leptons are fundamental particles 
(contrary to hadrons which consist of quarks and gluons), these colliders are 
especially apt for precision measurements of particle properties. 0] 

The press releases about evidence for the production of the quark-gluon 
plasma in heavy ion collisions at the SPS were first announced in 2000 
Only indirect observations are possible since the QGP formed in the initial 
stage is quickly turned into a system of hadrons - a process called "hadroniza- 
tion". Detection of single quarks is not possible due to color confinement, 
a property of strong interactions at low energies, according to which quarks 
must always combine to color-neutral hadrons before being able to travel 
to the detector. A number of experiments was necessary before scientists 
dared to announce the discovery of the QGP: The theoretical analysis of the 
measured hadron abundances resembles a state of "chemical equilibrium" 
at a temperature of about 170 MeV which marks the quark-hadron transi- 
tion. In particular, there is an observed enhancement of hadrons containing 
strange quarks by a factor of 2 to 15 (depending on the hadron), relative 
to proton-induced collision, and an observed suppression of the charmonium 
states J ftp and ip' [Zj which contain charm quarks with masses of about 1.2 
GeV, much higher than the transition temperature of about 170 MeV. Direct 
observations of the QGP via electromagnetic radiation at the SPS are diffi- 
cult due to high backgrounds from other sources. Higher energies at RHIC 
and the LHC may allow for better observation of the plasma radiation and 
enable detailed studies of the early thermalization processes and dynamical 
evolution of the quark gluon plasma. These experiments might also help 
to study the order of the phase transition and to locate the tricritical point 
within the QCD phase diagram, where the first-order transition changes into 
a crossover behavior. 
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T 



LHC, 
RHIC 



SPS 



quark-gluon plasma 



170 MeV 




10 MeV 



308 MeV 



Abbildung 1.1: Schematic QCD phase diagram of strongly interacting mat- 
ter. 



Depending on temperature and particle density, matter exists in various 
phases. Just like water can exist in the familiar phases of solid, liquid, or 
gaseous state with phase transitions between them, so will it enter a new 
phase of a plasma (ionized gas) at high enough temperature to break up 
chemical bonds and to ionize hydrogen and oxygen atoms. If the tempera- 
ture still increases, matter undergoes another phase transition where quarks 
and gluons, which are the constituents of protons and neutrons, leave their 
nuclear confinement to form a quark gluon plasma. On the lower end of 
the temperature scale, matter may exist in the phase of a Bose-Einstein 
condensate, superfluids, or superconductors. 

The fundamental theory describing quarks and gluons is the theory of 
strong interactions, quantum chromo dynamics (QCD). Unlike electromag- 
netism, which becomes stronger as the mutual distance decreases, strong 
interaction gets weaker. As a consequence, quarks, the fundamental parti- 
cles of QCD, are confined in packages of three (called "baryons", like protons 
and neutrons) or two ("mesons", e.g. pions) at low temperatures. Any at- 
tempt to separate confined quarks rather produces a new quark-antiquark 
pair from the gluon field in-between that combines with the separated quarks 
so that no single quark or net color charge can be observed. 

In figure fTTTI we see a schematic view of the QCD phase diagram. Strongly 
interacting matter can exist in three distinct phases, depending on temper- 
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ature T and chemical potential fi: the hadronic phase (confined to hadrons 
as a hadron gas or nuclear matter), the quark-gluon plasma, and the color- 
superconducting quark-matter. Our world is located on the lower line of the 
diagram around the phase transition of 308 to 313 MeV, in the form of hadron 
droplets (nuclei). Temperature T is given in units of MeV which can be con- 
verted to Kelvin using the Boltzmann constant ks = 1.3807 x 10 _23 J/K = 
8.617 x 10~ 5 eV/K. For a nice day of 27°C we get 27°C w 300 K w 25.9 meV. 
The core of our sun has a temperature of about 1.6 x 10 7 K « 1.35 keV. To 
get to the quark-gluon plasma we still need to increase temperature by five 
orders of magnitude to at least 170 MeV. 

The other axis of the diagram shows the quark chemical potential fx. The 
term "chemical potential" may be misleading since we deal with temperatures 
and densities far from ordinary chemistry, but it denotes a useful quantity 
that is applicable to any thermodynamic system: The chemical potential is 
the change in the energy of the system when an additional constituent par- 
ticle is introduced, with the entropy and volume held fixed (or equivalently 
the change in the Helmholtz free energy with temperature and volume held 
fixed). For systems containing different species of particles, there is a sepa- 
rate chemical potential associated with each species. If the particles can be 
transformed into one another, particle species with higher chemical potential 
will transform into a species with lower chemical potential, releasing heat. 
Therefore in an equilibrium state of several species, the chemical potentials 
of all species must equal each other. At zero temperature, (infinite) nuclear 
matter has a ground state at a quark chemical potential of fj, = 308 MeV. 
This is the nucleon mass ~ 939 MeV minus the binding energy 16 MeV 
divided by three for the average energy per quark or quark chemical poten- 
tial. In our everyday life we find nuclear matter density only in atomic nuclei 
- the droplets in which quarks are confined. Only if we increase the pres- 
sure of the system so that all droplets start touching and overlapping each 
other, we can build a phase of infinite nuclear matter. Up to this point, the 
pressure of the system at zero temperature is zero, because any compression 
just means decrease of the void space between the droplets. If the density 
is increased beyond the point where the droplets start to touch, pressure 
will start to increase, too. This is the so-called liquid-gas phase transition 
between the hadron gas and nuclear matter. 

The pressure always changes continuously if we change from one phase to 
another, but if its first derivative with respect to T is discontinuous, we call 
this a phase transition of first order. If the first derivative of the pressure 
is continuous, but its second derivative is discontinuous, we have a second 
order phase transition. Finally, if the pressure is continuous to all orders, but 
there is a rapid change in pressure, we can talk of a crossover. Usually, points 
of phase transition in the diagram form a line that starts at one of the axis 
of the diagram and may terminate at a critical point. Also the first-order 
liquid-gas phase transition has a critical point for a temperature of about 10 
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MeV where the transition becomes second order. Above this temperature 
one cannot distinguish between the gaseous and the liquid phase. 

Below a temperature of 170 MeV and chemical potential of about 350 
MeV, quarks stay confined in hadrons. If we further increase temperature 
and/or quark density, we arrive at the quark-hadron phase transition at 
which quarks cease to be bound in individual nuclei and become either a 
deconfined quark-gluon plasma or a color superconductor. For a range of 
medium quark chemical potential, the phase transition is of first order, but 
for higher temperatures and smaller chemical potentials there is a critical 
point below which we only have a crossover, as suggested by lattice calcu- 
lations. Particle detectors like SPS, RHIC or LHC probe will probe the 
quark-gluon plasma in a range around this tricritical point. This is also the 
path that our universe took in early times, some 20-30 [is (« 170 MeV) after 
the big-bang, starting with high temperature and a quark chemical poten- 
tial of about 10" 6 MeV passing the quark-gluon phase transition down to our 
universe temperature of 2.725 K [H] and an average quark chemical potential 
of about 313.6 MeV [Hj. 

For small temperatures and high chemical potential we enter the phase 
of a color superconductor ^QJ^^- This phase has been studied extensively 
during the last decade, and its theoretical foundations include concepts taken 
from the theory of ordinary superconductors: In ordinary superconductors, 
electrons are bound together by phonon interactions into Cooper pairs. Color 
superconductivity occurs because in the anti-triplet channel there is an at- 
tractive interaction between two quarks at the Fermi surface ^2jE3- These 
quarks then condense in the new ground state of the system by forming 
Cooper pairs. Depending on T, color superconductivity can again be classi- 
fied in different phases: If only up and down flavors are involved, we have a 
2-flavor color superconductor (2SC), for three flavors including the strange 
quark we have a color-flavor-locked (CFL) phase. Other phases studied in- 
volve a color-spin-locked (CSL) phase and a polar phase. Matter may further 
occur in the form of crystalline LOFF (named after Larkin, Ovchinnikov, 
Fulde and Ferrell) where Cooper pairs with nonzero total momentum are 
favored Color superconductivity may exist up to a temperature of 6 to 
60 MeV above which we enter the quark-gluon-plasma phase. 

For small to medium temperature ranges, the QGP may be described 
as a non-Fermi liquid. It is called "Non-Fermi", because its behavior cannot 
entirely be described by the classical (Landau-) Fermi liquid theory of non- 
interacting cold fermions. This is due to transverse gauge bosons (transverse 
gluons) which are not screened at small temperatures due to the properties 
of dynamical screening. We have long-range, quasi-static interactions which 
give anomalous contributions to thermodynamic quantities like the specific 
heat. Such non- Fermi liquid behavior might have influence on astrophysical 
calculations, for example the cooling rate of proto-neutron stars ^Hj- Non- 
Fermi liquid behavior also appears in solid state physics experiments, for 
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example in the specific heat of the YbRli2Si2 crystal [TTj I18j. 

1.3 Outline 

This thesis is organized as follows: 

In chapter 2 we will introduce large Nf QCD where the number of quark 
flavors Nf is sent to infinity while the effective coupling g e d oc g 2 Nf stays 
of order 1. We will see that this theory is exactly solvable in the effective 
coupling at next-to-leading order of the 1/Nf expansion. Also, the scale 
dependence is completely solvable to this order. Therefore this theory is 
particularly suitable for testing other perturbative and non-perturbative ap- 
proaches against this limit. In contrast to lattice gauge theory, large Nf can 
be readily extended to finite fi. 

In chapter 3 we will explore the non-Fermi-liquid regime of the large 
Nf theory. We will derive the pressure and the specific heat perturbatively 
and show that beyond the leading logarithm there are anomalous fractional 
powers and we calculate their coefficients up to order T 7 / 3 in the specific 
heat. We will see that fractional powers appear from the transverse gluon 
propagator whose calculation is carried out in careful detail, but both, longi- 
tudinal and transverse contributions to the pressure are needed to complete 
the leading logarithm of the specific heat. 

The final chapter 4 gives a summary and an outlook. 



Kapitel 2 

Large Nf 



2.1 Introduction to Large Nf QCD 
2.1.1 Why Large N f l 

Large Nf QCD is quantum chromodynamics (QCD) with the number of 
quark flavors Nf assumed to be large (Nf — > oo) and the strong coupling 
a s = 5 2 /(47r) — ► assumed to be small so that the effective coupling (which 
we will define later) g^ s ~ 9 2 Nf ~ a s Nf ~ 0(1) stays of order 1. Also, the 
number of colors N c stays small N c ~ 0(1). 

Why do we want to study Nf — ► oo? After all, in our everyday world we 
only experience two quark flavors, that are the "up" and "down" quarks with 
masses of the order of 1-10 MeV. For the quark gluon phase transition at a 
temperature of about 170 MeV, only one more quark flavor, the "strange" 
quark with a mass of about 80-155 MeV jTHI, will play a considerable role 
in the thermal state. Even if we consider the other known quark types 
"charm", "bottom" (also called "beauty"), and "top" (which is barely called 
"truth" anymore) with masses beyond a GeV, we end up with a total of 
Nf < 6 quark flavors. Clearly, the aim is not to realistically model QCD, 
but to study the theory in a well defined limit that turns out to be exactly 
solvable up to next-to-leading order in an 1/Nf expansion while containing 
essential physics of the full theory. 

Ideally, we would like to solve full QCD, but so far only approximate nu- 
merical results from lattice simulations are available (at least for small fJ,/T), 
but with not completely understood systematic uncertainties. For tempera- 
tures far above the intrinsic scale of QCD, T 3> Aqcd, we expect the coupling 
g to get small (20] . Due to this property of strong interaction, called asymp- 
totic freedom, we can expand the theory in terms of small g. It turns out 
that a strict perturbative expansion of thermodynamic parameters like the 
pressure in terms of the small coupling g would only work for ridiculously 
high temperatures and fails to work for the interesting region of tempera- 
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ture values slightly above the phase transition T c or the QCD scale Aqcd: 
Successive orders in perturbation theory show very poor convergence. But 
this is not the only problem: Hot QCD has an intrinsically nonperturbative 
scale in the magnetostatic sector at wave numbers k ~ g 2 T. For the pressure 
this means that perturbation theory ceases to work at order g 6 T 4 [2T1I22] . 
For other observables, this might even mean stricter restrictions to an ex- 
pansion in g: Nonperturbative corrections are suppressed only by a power 
of g 2 for photon emissivity [22], a power of g for the Debye screening length 
[231 ES) or even not suppressed at all but appear at leading order in the case 
of the baryon number violation rate |26| I27j. But even in theories without 
the problem of a nonperturbative magnetostatic scale like in hot QED, the 
poor convergence of the perturbative series remains [2*% 1 l2li | l37) | l3~ H 132] . 

There have been attempts to improve the situation by reorganizations or 
partial resummations of the perturbative expansions of hot QCD. In scalar 
models where one already meets similar difficulties of poor convergence |33j . 
Karsch et al. [B3] proposed to keep a screening mass unexpanded at any 
given order of the loop expansion and to fix this mass by a stationarity 
principle, which technically means subtracting a mass term from the bare 
Lagrangian and adding it to the interaction part. This successful approach 
was extended to QCD by Andersen et al. [HH1 ESI EH EH ESI EH] by replacing 
the simple mass term by the gauge-invariant hard-thermal-lop (HTL) action 
|411 142j . a method which they termed HTL perturbation theory. Another 
approximation method is based on ^-derivable approximations ^3 EH for 
2PI skeleton diagrams advocated by Blaizot et al. @S1 ESI EH EE1 EH] Here, 
the starting point is an expression of the thermodynamic potential in terms 
of dressed propagators, where bare propagators are functionals of the full 
propagators that have to satisfy a stationarity condition. Truncation in this 
scheme does not happen in terms of the coupling g, but one resums diagrams 
to a given loop order of the 2PI skeleton diagrams. 

All of these approaches show improved convergence properties, but it is 
difficult to tell to what extent these resummations can predict the correct 
behavior of full QCD. Apart from the necessary condition that these theories 
should coincide with perturbative QCD at very small couplings (correspond- 
ing to temperatures far beyond the scope of current collision experiments), 
only lattice calculations (SOI EH could give truly independent indications as 
to which resummation scheme to trust. Lattice theory on the other hand 
works more like a "black box" for theorists and provides only limited insight 
into the actual underlying physics, like screening effects or collective modes. 
One has therefore been looking for simpler theories that could be solved ex- 
actly, so that other resummation schemes could be tested within that simpler 
framework. One such theory is large- N scalar field theory with 4> A interac- 
tion (S2HSH] with an exactly solvable N — > oo limit. Another limit that has 
been studied in the past is the large N limit of a scalar field theory in 6 
dimensions with cubic interactions [HE] that mimics QED with N flavors. It 
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is richer than large N (/> 4 -theory in that self-energies vary with momentum 
and wave-function renormalization is needed just like in full QCD. Yet, this 
theory involves instabilities which prevents one from direct comparisons to 
QCD. A more appealing test bed is large Nf QCD. 

Besides being exactly solvable up to next-to-leading (NLO) order in 1/Nf, 
large- iV/ QCD contains the same complicated frequency and momentum de- 
pendent gauge field screening and damping as full QCD and its perturbative 
series shows similarly poor behavior as full QCD. One should keep in mind 
though that since large Nf QCD only resembles part of the full QCD theory, 
we do not expect large Nf to give any useful predictions for a small number 
of flavors. In fact, with a running coupling that leads to a Landau pole (like 
in QED) and excludes the possibility of confinement, large Nf at NLO is 
quite different from full QCD. Also, the leading order contributions to full 
QCD in an expansion in the coupling g are numerically dominated by sub- 
leading powers of 1/Nf for small Nf, as we will see in the next section. Yet, 
large Nf can be regarded as an important test theory, and any resummation 
scheme had better get close to the exact NLO result in the large Nf limit. 
On the other hand, even if a theory turns out to reproduce the right large 
Nf limit, that does not guarantee that it will correctly predict the full QCD 
result. After all, large Nf is a test theory, but contains a lot of the physics 
of full QCD, so if we can do this test, let's better do it. 

2.1.2 Comparison to strict perturbative expansion in g 

It is instructive to see which terms of perturbation theory survive in the 
large Nf limit. Strict perturbation theory has a long tradition. An overview 
over the known coefficients is given in figure 12.11 The leading order result 
to the pressure of perturbative QCD is just the Stefan-Boltzmann pressure 
of non-interacting quarks and gluons. The first non-trivial contributions to 
this result proportional to a s were calculated independently by Shuryak and 
Chin in 1978 (H3JEH] by calculating one- and two-loop vacuum graphs. The 
next order already turned out to be non-analytic in a s and required the 
resummation of an infinite subset of diagrams. This was done by Kapusta 
to order o?J 2 (HE] and Toimela to order a 2 s \na s [57] . Subsequent orders 
demanded new tools for evaluating sum-integrals by Arnold, Zhai and Kas- 
tening (23123 EH], or the introduction of dimensional reduction [59l 160} IfiT] 
by Braaten and Nieto [H21- The last perturbative series coefficient that has 
be calculated so far for full QCD is the lna s coefficient by Kajantie, Laine, 
Rummukainen, and Schroder (HSHM]- The order coefficient is completely 
non-perturbative as pointed out by Linde (2^ and Gross, Pisarski, and Yaffe 
(221- Recently the perturbative result has been extended to finite chemical 
potential by Vuorinen [65] . 

In this perturbative series the effect of large Nf is straightforward to be 
seen: Since we keep the quantity a s Nf ~ 0(1) of order one while we send 
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Perturbative calculations of QCD at high T/fi 
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+(1 + ir)^ {-799.2 - 21.96JV) 



1.93JVJ 
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{ 1139.8 + 65.89N/ + 7.653iV| 
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Planck 1900 
Shuryak/Chinl978 

Kapusta 1979 
Toimelal983 
Arnold&Zhai 1995 

Zhai&Kastening 1995 
Braaten&Nieto 1996 

Kajantie, Laine, 
Rummukainen & 
Schroder 2002 



{ ?!+?A r /+?iV 2 + C (fi)Nf} (^) 3 + 0(^)1} 



* extension to 7^ : Vuorinen 2003 

* C(/2) = 0.040. ..for /2 = vrT (A.I.&A.Rebhan 2003, see section [MM 

* ?!: completely non-perturbative (Linde 1980; Gross, Pisarski & Yaffe 1980) 

Figure 2.1: Known coefficients for the pressure of perturbative QCD by 
2003. The coefficient marked as ?! denotes the breakdown of perturbation 
theory in the order a 3 . Only a few coefficients of this series contribute 
to the NLO contribution of the large-iVj limit, namely those of the form 
a™NJ. Particularly, terms containing a logarithm of the coupling lna s do 
not survive the large- Nf limit. The large- Nf coefficient C(p) of order Nja^ 
has been extracted numerically from the exact NLO result. 
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Nf — > oo and a s — ► 0, the leading order (LO) contribution to the pressure 
is just given by the contribution proportional to Nf. P = 7ir 2 NfT 4 /60. 
Actually, this contribution is infinitely large for Nf — ► oo, but we can always 
subtract this known leading order contribution to obtain a finite correction 
at next-to-leading order (NLO). From fignre l2~Tl we see that now all orders in 
the coupling a s contribute to the large Nf NLO result which are of the form 
a p s N p with p = 0,l,|,2,|. A graphical comparison of the partial sums 
of the series against the exact large Nf NLO calculation will be presented 
below in figure f2T8l 

2.1.3 Outline 

The pressure in large Nf QCD was first calculated by Moore for zero chemical 
potential [HSI 1 and successively extended to finite chemical potential |68| . 
In the following sections we will present the exact large Nf result for the 
thermal pressure and some derived quantities thereof like the quark number 
susceptibilities or the entropy. The theory inevitably contains a Landau 
pole. We will study its effect on the result and show for which range of 
the coupling the effect of the Landau pole is negligible. At zero chemical 
potential we compare the exact large Nf result of the pressure and the quark 
number susceptibilities with known results from thermal perturbation theory 
(SHI EE1 EH] obtained at small chemical potential where dimensional reduction 
j^nilZnil^ is applicable. We verify the recent three-loop result of Vuorinen 
[7T] on quark number susceptibilities numerically in the large Nf limit as well 
as a numerical coefficient in the pressure at zero temperature obtained long 
ago by Freedman and McLerran [221 HE] ■ For small values of the coupling 
g, our numerical accuracy allows us to extract a number of perturbative 
coefficients at order g 6 that are not yet known from analytical calculations. 

It is remarkable that large Nf can be safely calculated all the way down 
from large temperatures and zero chemical potential to large chemical po- 
tential and zero temperature for effective couplings g^ s < 20, which is a 
main advantage compared to lattice gauge theory. For a long time, the so- 
called "sign problem" prohibited calculations of fermions with finite chemical 
potential on the lattice. Only recently there has been important progress re- 
garding calculation of thermodynamic quantities within lattice gauge theory 
[TH 175} IT6l 1771 1781 179] . Using the large Nf limit, we can test the scaling 
behavior noticed in lattice calculations by Fodor, Katz, and Szabo [25] and 
find that it breaks down rather abruptly at \x q > irT. We can also test the 
range of applicability of dimensional reduction by comparison to perturba- 
tive calculations at finite temperature and chemical potential by Vuorinen 
|65] . where we study the effect of the choice of the renormalization scale. 



1 The first published version contained an unfortunate coding error which was revealed 
and corrected by an independent calculation involving the author of this thesis |(i7j . 



14 KAPITEL 2. LARGE N F 

o 

(A) 

oJXo cvTLo rJXo 0JX9 

i Jri b + i = i b 

^opj %jxT V^T %jxT 

(B) 

Figure 2.2: Diagrams contributing to the large Nf limit. At leading order 
(LO) of the order O(Nf) there is only the fermion loop (A) contributing. 
At next-to-leading order (NLO) of order O(l) there is an infinite number of 
diagrams contributing (B) which can be resummed by the Schwinger-Dyson 
method. 

At small temperatures T < /i we enter the region of a non-Fermi-liquid 
behavior, which will be studied in detail in the next chapter. 

Further discussions on the large Nf limit can be found in the literature 
by Peshier jHH] on the quasiparticle picture (HB E2J EE], who points out 
the large differences between large and small Nf QCD with respect to the 
different strong coupling behavior and questions the immediate significance 
of a comparison within large Nf NLO. A detailed discussion of the HTL- 
quasiparticle picture of QCD with respect to large Nf and its implications 
on 2PI ^-derivable approximations is given by Rebhan 

2.2 Pressure at Large-iV/ QCD 

In the following we want to calculate the thermal pressure in the large Nf 
limit. As we already mentioned there are indeed two quantities that are in- 
volved in the limit: The number of flavors Nf and the coupling g. By forming 
the limit Nf — > oo we actually keep the effective coupling g 2 s ~ g 2 Nf of order 
1. In this sense, a suppression by a factor 1/Nf is equivalent to a suppression 
by a factor g 2 . We will see how this provides a very strong ordering principle 
for diagrams in the large Nf expansion. There are no special assumptions 
on the number of colors. We will calculate the pressure at next-to-leading 
order (NLO) where the theory can be solved exactly, that is to all orders in 
the effective coupling, apart from possible difficulties introduced by a Lan- 
dau pole that we will discuss later. We will follow the original approach of 
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references [Su IHZj and define the effective coupling as 



5 2 iY/, QED . 



where Cp is the quadratic Casimir of the representation containing the 
fermions and dp and a! a are the dimensions of the fermionic and adjoint 
representations. This allows us to treat the two cases of massless QCD and 
ultrarelativistic QED the theory at the same time. Since the coupling runs in 
these theories, we need to know the (3 function. The one-loop (3 function for 
QCD, that is a SU(N) gauge theory with Nf fermions in the fundamental 
representation, is given by 



udg ud(g 2 ) n . 9 , 2o 4 /ll 2 \ , 

For sufficiently small Nf, that is Nf < ^N, the /3 function is negative, which 
means that QCD is asymptotically free. In the limit of Nf — ► oo we see that 
this behavior changes as the sign of (5(g) gets positive and large Nf therefore 
is no more an asymptotically free theory. The (5 function then becomes 

This relation turns out to be exact at leading order in 1/Nf, even if we 
add higher loop beta function contributions (see for example reference |85| ) 
to the one-loop result (|2,2,2|) because they are diagrammatically suppressed 
at least by a factor of 1/Nf oc g 2 . We can solve the renormalization scale 
dependence exactly, giving 

+ ^Bf±. (2.2.4) 



This theory contains a Landau pole of the order of Al ~ ^exp(67r 2 /^ ff ). 
Following reference |HE] we define the Landau scale Al such that the vacuum 
gauge field propagator diverges at Q 2 = A^. The vacuum gauge field prop- 
agator is defined as D~^. = q 2 — q$ + Il vac with II vac from equation (|2.2.9|) 
below, which leads to 

A L = / UMse 5/6 e 67r2 / 9 'ff^ MS ) . (2.2.5) 

This Landau singularity means that the exact definition of the theory is 
ambiguous unless the UV completion is specified. But that should not hinder 
us from doing finite temperature calculations since the ambiguity for the 
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thermal pressure at NLO is suppressed by a factor (max(T, /i)/AL) 4 . As 
long as we stay well below the Landau pole with temperature T and chemical 
potential fx, we can expect to get results which are not afflicted by the UV 
incompleteness of the large Nf theory. 

The pressure can be calculated by computing the diagrams contributing 
to the free energy, which is the trace of the sum of all 1PI vacuum bubble 
graphs. At leading order in the 1/Nf expansion we only have one diagram - 
that is the bare fermion loop in figure l2~2l To obtain the thermal pressure, 
we subtract off the vacuum part of the diagrams. The leading order pressure 
is thus given by 

Ifc, = N,Tr ( T £ / In - / 1„ @) . (2.2.6) 

\ n,odd I 

It is well-known how to calculate this contribution for finite temperature and 
chemical potential and the pressure is given by 

/7tt 2 T 4 u 2 T 2 li 4 \ 
P LO = NN f ( + ^— ~\~ — — J . (2.2.7) 

u ; V 180 6 12vr 2 J v ' 

The leading order contribution in the pressure is thus of order Nf and strictly 
speaking infinitely large in the limit of Nf — > oo. Of course, we can always 
think of this as the leading contribution to an expansion in terms of 1/Nf 
with large but finite Nf. 

The next-to-leading order (NLO) of this expansion is proportional to 
N® = 1. We can calculate it from the gauge boson loop with an arbitrary 
number of fermion loop insertions (plus corresponding counterterm inser- 
tions), which are all of order N®. For each fermion loop insertion we get a 
factor Nf and a factor g 2 from the two vertices, thus just our effective cou- 
pling g 2 Nf oc g 2 s ~ O(l). Any additional bosonic or ghost insertion would 
give factors of g 2 or higher, without introducing new fermion loop factors of 
Nf, and are thus suppressed by at least a factor of g 2 oc 1/Nf. Therefore 
this gauge boson loop is all there is to NLO. It can be resummed by the 
standard Schwinger-Dyson resummation. Its contribution to the pressure is 
(again subtracting the vacuum part) 



NLO 



-\^[ T E / ^H[D 1 r(Q) + ^(Q)) 

\ n,even * ' 

- / Hin^rwi+cw))) (2-2.8) 

where the trace runs over group and Lorentz indices. The vacuum part of 
the gauge-boson self energy is given by 

iCc(Q) = - Jj (v^Q 2 - Q'Q") ^ " |) . (2-2-9) 
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and the one-loop bosonic self-energy can be calculated via 

n G (Q) = f^n^CQ) - [vac] = 2U T (Q) + YL L (Q), (2.2.10) 
U H (Q) = n 00 (Q)-[vac] = ^U L (Q) (2-2.11) 

where Hq and Hh will be calculated in the next section in equations (|2,3,21|) . 
(I2.H.22H . and H2.H.2HI) . These cannot be given in closed form except for 
their imaginary parts |67j . but are represented by one-dimensional integrals 
involving the fermionic distribution function nj. Using this decomposition 
of the self energy, the trace of the resummed propagator can be rewritten as 

Trln([^Y"(Q)+n t Al h(Q)) = 21n(Q 2 +n r +n vac )+ln(Q 2 +n L +n vac )-ln(Q 2 ) 

(2.2.12) 

if we choose an appropriate gauge in which Dq 1 factorizes in the same way as 
II, as for example in the Feynman gauge [D^ 1 ]^ (Q) = rj^Q 2 . Any gauge 
dependence drops out to our order of interest, since we only regard the 
difference between thermal and vacuum contributions. We can now perform 
the sum over Matsubara frequencies in the usual way of replacing it by a 
contour integral that is deformed appropriately. The result is then given by 

m 

-^Imln(g 2 - ql + II vac ) 

+ ([n b + -]Imln( _ , ) - -Imln( g2 _ ^ ) 

with the bosonic distribution function nb(co) coming; di- 

rectly from the sum over even frequencies qo = 2mrT. When evaluating the 
integrals above exactly by numerical means, we can safely integrate parts 
proportional to in Minkowski space, since those are exponentially ultra- 
violet safe. For terms without ra& more care is required. We will refer to the 
former contributions as "n^'-parts and the latter as "non-?V contributions. 
The "non-rif," contributions are potentially logarithmically divergent, unless 
a Euclidean invariant cutoff is introduced [00] ■ We will discuss this point in 
more detail in section l2~4l 

For non- vanishing chemical potential /i we use the fermionic distribution 
function 



n fV°> T >ri = \ (TO + e (*+d/r + 1 ] 



which enters via the gauge boson self-energy expressions Ihr and Ul- 
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2.3 Leading-order gauge boson self-energy 

The leading-order gauge boson self-energy in the large Nf limit is diagram- 
matically represented by a gauge boson propagator containing one fermion 
loop insertion. In the imaginary time formalism it can be written as 

ttTP^p) = 2^V/^/^/^Po(^o,k)p (go,q) 

x " (fco) - n ^°) (kPtf + gUtf _ g^k a q a + g^m 2 ) (2.3.1) 
ko - Qo - iu n 

with q = k — p. The constants g and tt can be adjusted to fit various au- 
thor's conventions (see also the conventions used in appendix El . We follow 
Weldon's conventions [86j by setting g 2 /tt = —1. The fermionic Boltzmann- 
factor is given by n(ko) = l/(e k °/ T + 1) and the spectral function is given 
by p (^o,k) = 27re(k )6(k 2 - k 2 - to 2 ) = ^(S(k - e k ) - 5(k + e k )) with 

e(ko) = ko/\ko\ and e k = Vk 2 + to, 2 . We can find a spectral form of the self 
energy by using 

n^p)=r^^^) ( 2 .3.2) 
J-oo 2ir z-p 

which is valid for any complex z [HH eq (3.1.36)]. The spectral density is 
then given by 



/d 3 k f°° dko f°° dqo 
J-oo^ J-oo 2^ P ° (fc0 ' k)Po(90 ' q) 

x (n(fco) - n(q )) 2vr5(p ~ k + go)/'* (2.3.3) 

with 

If" = (k^q u + q»k v - g» v k a q a + g^m 2 ) . 

The spectral density is manifestly real which will simplify our calculations. 
Since the spectral density will play an important role, we will try to push 
its calculation as far as we can. After integrating over go, we are left with 
with an expression depending only on W = (fco,k) and p^ = Q?o>p)- A 
transformation of variables k^ — > —k'^ + p^ in the part proportional to 
n(go) = n(ko — po) leaves all other parts invariant up to minus signs. It 
basically replaces n(ko — po) by —n(—k' Q ) = n(k' ) — 1 so that we can replace 
(n(ko) — n(qo)) by (2ra(fco) — 1) in equation l|2.3.3p . If we use z = cos# as 
the angle between pk = pkz, we can express 

(v l. \ o (i \ 1 xi k 2 +p 2 - (k - po) 2 
Po{ko ~ Po, k - p) = 2vre(/c - Po)-^o(z ^ ) (2-3.4) 
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and integrate over z in J ?|Mp = t^w J-i dz J °° /c 2 d£;. Analogous to Weldon 
(HSl we can then define 

la = g' MU I^ = -k 2 + k 2 + P 2 - P 2 , (2.3.5) 
I H = u"u u I^ = ^(k 2 -p 2 + {k -po)(3k - Po )) 

with = (1,0,0,0) and integrate over ko- This is the last integral that we 
can do analytically. The spectral density turns out to be antisymmetric in 
ko meaning 7r(fco,k) = — ir(— A;o,k). The final result for the spectral density 
of the self-energy can then be written as a sum of two parts that only differ 
by the sign of po 

^x(P0,P) =T^x(P0:P) ~ ^x{~P^P), (2.3.6) 

xe(k - po)0{\k - p\ < \k - p \ < \k+p\) 

with X = G or H as in 

I G = I G {k = k) = pl-p 2 , (2.3.8) 
I H = lH{k = k) = ^(2k+p-p )(2k-p-po). 

The ^-function stems from the integration over z with — 1 < z < 1 and its 
usage here means 6>(true expression) = 1 and 9 (false expression) = 0. 

Self-energy in the complex energy plane 

Starting from the spectral form l|2.3.2|l we can calculate the self-energy in 
Minkowski or in Euclidean metric or anywhere between in the complex plane. 
To compare our results to the real-time formalism we start with the self- 
energy with Feynman prescription 

Hf(po,p) = U(p + ip e,p). (2.3.9) 

Separating = ~ — iir5(x) with P denoting the principal value, and 

knowing that 7r(A;o, k) is a real function, we can easily separate this equation 
into real and imaginary parts 

~ f°° dv' P 
Ren F (po,p)=/ ^(p' ,p) r , (2.3.10) 

J-oo 27r P0 ~ Po 

Imn F (p ,p) = -^e(Po)7r(p ,p)- (2.3.11) 
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In the general case we can write 

II(a-Mo,p) = / — tt(po,p) — — — 
J-oo 2tt a + ib 



Po 



I 



dp ( a-p . b 

2vr \(a-p y + b 2 (a - p y + b\' 



where we separated real- and imaginary parts. By using the antisymmetry 
property of the spectral density l|2.3.6|l we can finally write 

U(a + ib,p) = —^(po,p)[- 2 2 - - — ■ 2 2 

J 2vr \(a-p y + b 2 (a + p ) 2 + b l ) 

-i I ^-ir + (p 0l p) ( - \ - - — ■ — b ) (2.3.13) 

J 2vr \(a-p y + b 2 (a + p ) z + b l ) 

where all integrals to be evaluated are real. Particularly for the case a — > 
and i = w / we recover the result for the Euclidean metric 

U Eucl (iLo,p) = U{iuj + e,p), (2.3.14) 
Ren Eud (tu;,p) = / — ir^(p ,p) r i2j _„2 

J — oo ^ 



-oo 

(•OO 



LO z +p% 



= l^^^fpj^ (2 - 3 ' 15) 

lmU Eucl {iio,p) =0. (2.3.16) 

Minkowski space results 

Plugging our self-energy l|2.3.6p in the expressions (|2.3.10|) above we get the 
Minkowski space results. First we note that we can write the integral over 
Po that we have to deal with in a simple way since p and k are positive: 

dp e(k -p )9(\k - p\ < \k - p \ < \k + p\)f(p ) 



/ rp r2k+p 

6(k-p)[ I dpofipo)- / dpof(p 

\J~p J2k-p 

/ r2k—p p2k+p 

+6{p -k)i^j dp f(p ) - j dpofipo) 

/p i-2k+p 
dpofipo) - / dpofipo). (2.3.17) 
-p J2k—p 
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The resulting integral for the real part can then be written as 
ReU x (po,p) = 



dp' - / dp' 

p Jlk—p 



(2.3.18) 



Po ~ Po Po + Po 



T J I x(k,p' ,p). 



Using 



/ dp' 



Po - Po Po + Po 



T ) lG(k,p' ,p) 



and 



/«*( 



P o ~ (Po ~ p) (log(bo - Po\) + log(b' +po\)) 



T ) lH{k,p' ,p) 



(2.3.19) 



Po ~ Po Po + Po 



= 2 \- 8kp 'o ~ p 'o ~ ( 2k ~ P ~ Po)( 2k + P ~ Po)log(\p -pol) 

_(2/, - p + Po )(2k + p + Po ) log(|p + po|)] (2.3.20) 

we finally arrive at the real part of the self-energy that is stated in Weldon's 
paper [86] 



Ren G (p ,p) 



2p 



2k + po + p 2k - p + p 



2k + po — p 2k — po — p 



and 

ReIl H (po,p) 



g 2 } g 2 N f 

7T / 2vr 2 

1 - Q log 
P 



o 



dk I n{k) 



Po+P 



Po-P 



(2k + po + p){2k + po - p) 
4p 

(2k -po - p)(2k - po+p) 
4p 



log 



log 



2k + po + P 



2k + po - p 
2k-po-p 



2k-po + p 



(2.3.21) 



(2.3.22) 



Note that Weldon already subtracted the T = vacuum part and replaced 
n(k) — \ by n(k). For the imaginary part no further calculation is necessary 
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and we can simply relate the imaginary part of the self-energy to its spectral 
density (l2~Tfil) 



Imn x (po,p) = ~2 e (Po) {^xiPOiP) -^x(-POiP)) 



(2.3.23) 



We are left with integrals over k that can be solved analytically. We rewrite 
the integral over k for the sign- and the ^-function as 



dke(k -p )6(\k - p\ < \k - p \ < \k+p\)f(k) 



oo 

e(j?-p$) I + f{k)dk-e{p -p) 

I PQ+P 
2 



2 



PQ-P 
2 



f(k)dk. 



(2.3.24) 



The function f(k) is one of the following f(k) = n(k), kn(k), or k 2 n(k) if 
we already subtract the vacuum part from n(k) — |. This leaves us with the 
following integrals 



+ Tlog(e x/T + l), 



Ftix) = 


poo 

/ n(k)dk = - 

J X 


F 2 (x) = 


poo 

/ kn(k)dk 

X 


F 3 (x) = 


/•CO 

/ k 2 n{k)dk = 

J X 



n 2 T 2 



X 



+ xTlog(e^ T + 1) + I a Li 2 (-e B/T ), 



x 



+ x 2 Tlog(e x/T + l) 



+2xT 2 U 2 {-e x l T ) - 2T 3 Li 3 (-e x / T ) (2.3.25) 

with Li n (x) being the polylogarithm function. The functions Fi(x) have the 
property that Fi(x — > oo) — » and F.- L (x > 0, T — > 0) — > 0. The thermal 
contribution to the imaginary part of the self-energy can then be written as 



ImIlx(po,p) 



-Apo 



9 2 \ 9 z N f 



ir ) 2irp 



2 2\l„+,Po+P\ 



+0(po-p)[F$( 



rPO+P, 



(2.3.26) 



with X = G or H and 



F^(x) 



n(k)I G dk = {p 2 -p 2 )F l {x), 



(2.3.27) 
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vl -p 2 



[•CO 

F±(x) = / n(k)I H (k,± Po ,p)dk = T 2 Po F 2 (x) + 2F 3 (x). 

(2.3.28) 

Using symmetric and antisymmetric versions of these functions F% = (F^ + 
F^)/2 and = -F^)/2 we can rewrite equation 112.3.2KB as 



ImIIx(po,p) 



1 , ^W^/ 
^l( |P0+Pl )-^l( b °~ Pl 



(2.3.29) 



, / , \jrAf\P0+P\s / > A ,p -P 



This way of writing is especially convenient for calculating the G-part of the 
self- energy, since Fq = vanishes. 



Euclidean space results 

To obtain the results in Euclidean space we start from equation 1)2.3. 15j) with 
7r + . As in the case for the Minkowski metric, we can simplify the integral 
over po with sign- and ^-function according to equation (|2.3.17|1 and we are 
left with real integrals of the basic form 



/ 



d P0 ,. 2 ^\ lG{k,P0,p) = -Po + (uJ 2 +p 2 )log(uj 2 +p 2 ) (2.3.30) 



u A +p z Q 

and 



/ 



dpo 9 2P ° 2 lH{k,Po,p) = l;Po(8k - p ) - 4/cwarctan (— 
uj z + Pq 2 \ uj 

-±(Ak 2 -p 2 -u J 2 )\og{uo 2 +p 2 ). (2.3.31) 



For uj 2 + p^ > there is no danger of running into singularities. Using the 
integral limits for po from — p to p and 2k — p to 2k +p as in 1)2.3.17)) we can 
write the self-energy in the Euclidean metric as 



ReII G (^, p) = I dk[n(k)--) (2.3.32) 



X(4fc+ 2p YOg ^ + {2k+pY 
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ReU H (iuj,p) = [jt^t^L J dk L(k)-^\ (2.3.33) 

„, w 2 +p 2 -4£; 2 Lu 2 + (2k-p) 2 

2k H -. log — = —-. 

4p u) 2 + {2k+p) 2 

2kuj ( 2k - p p 2k +p 
arctan h 2 arctan arctan 

p \ U U> OJ 

This is in principle the result given by reference |66j up to the appearance 
of the arc tangents in the last line. The expression given there can be fully 
recovered by using arctan(z) = izlog(l — iz) — i«log(l + iz) for complex z 
and gathering the sum of logarithms into one logarithmic expression 

2k — p p 2k + p i 1 + tu)-ip\ 2 
arctan h 2 arctan arctan ~ log -r- — . (2.3.34) 

However this expression should be used with care. If one takes the principal 
branch of the logarithmic expression on the r.h.s., one obtains results that 
might differ from the arc tangent expression by some multiple of it. Solving 
for the branch-cut, that is where the expression inside the logarithm equals 
— 1, one finds that expression for real u, p, and k can be written as 

1 + 4fe2 2 2 



2 l+(^F 2v^^ 

(2.3.35) 

where one should take the principal value of the arctan and log expressions 
involved. As shown in equation (|2.3.16|) the imaginary part of the self-energy 
in the Euclidean space vanishes. 



Complex Result 

The result in the complex plane allows us to link results in Minkowski space 
and Euclidean space by great arcs. The safest way to obtain the result 
in the complex plane is to start with equation (|2.3.13|) . where all integrals 
involved are real, and the result for the self-energy is separated into real and 
imaginary part. Starting point are the following indefinite integrals 

F§ e ( P0 ) - IdpJ . a ~, P 2 \ h2 - ( "\ P 2 \ h2 )lo (2.3.36) 
J \{a-p y + b 2 (a+p ) 2 + b 2 J 

= — Pq — 2ab ( arctan — + arctan - ^° ) 

\ b b ) 

-\{a 2 - b 2 +p 2 ) {\og(b 2 + (a-po) 2 ) + log(6 2 + (a + Po ) 2 )) , 
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FlTiPo) = d P0 [- CT - , - , 2 , J / G (2.3.37) 

= —(a 2 — b 2 — p 2 ) ^arctan — + arctan — ~J~^j 
+ab (log(6 2 + (a - p ) 2 ) + log(6 2 + (a + p ) 2 )) ■ 

Analogous results for F§ e and Fjf 1 can be easily derived with any symbolic 
integrator. We combine these integrals according to equation l|2.3.17jl to 
form the result in general complex plane (X = G or H): 

(o 2 \ a 2 N f'°° / 1 \ 1 

-§-) ^/ J dk ln(k) - - j — (2.3.38) 



x 



F^(P) ~ F% e {-p) - F^{2k + p) + F# e (2fc - p)) 



Imn x (a + i6,p) = -(-Cj?J!l dk(n(k)-±)± (2.3.39) 
x (Fk m (p) - F}T(-P) ~ F I x m {2k + p) + FjT(2k - p)) . 

For b ^ the expressions in the logarithms are positive, and the arc tangent 
functions stay finite, so that there is no danger of crossing a branch cut or 
encountering other singularities. Therefore these expressions are safe to be 
used in calculating the arc between Minkowski and Euclidean space. 
Again, to subtract the vacuum part, just replace (n(k) — i) by n(k). 



2.4 Landau pole ambiguity 

In the pressure formula (12.2.1 HH that we intend to integrate, we had two 
different kinds of contributions that we called "rife" and "non-?V. We already 
mentioned that the former are safe to be calculated in Minkowski space: for 
large qo the bosonic distribution function suppresses the integrand exponen- 
tially by 0(exp(— qo/T)), while for large q but moderate qo, the thermal self- 
energy has an exponentially small imaginary part Imll ~ exp(— (q — qo)/2T) 
as can be seen from the exact result for the imaginary part of the self energy 
in Minkowski space (|2.3.26|) . knowing that the functions Fi from Ij2.3.25|) 
vanish exponentially as their argument increases. For the "non-nf," parts in 
(|2.2.13|) the situation is more subtle. 

The "non-n;," parts are not exponentially suppressed for large qo. Also, 
in the region where qo ~ q, the self-energy functions are only suppressed by 
1/q so that they might cause UV divergencies that should be studied. It is 
best to go into Euclidean space and study the relevant terms of the pressure 
function (|2.2.13|) . In the limit of small T 2 /Q 2 , the logarithms involved can 
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be expanded, giving 



2 In 



Q 2 + n vac + n T , . Q 2 + n vac + n L 



Q 2 ~\~ n vac 
2n T + n L 



+ In 



Q' 1 



2iil + n 2 



Q 2 + n vac 2(Q 2 + n 



vac 



+ ... 



(2.4.1) 



The functions Ut and ITl can be expanded using 1)2.2.10(1 and (|2.2.1ip and 
the expressions for RellG^iu;, p) 1)2.3.32(1 and ReIL}j(iu,p) 1)2.3.33)) . One can 
expand the two functions in the limit T 2 <C Q 2 by series expansion of the 
expression within the brackets in terms of small k and then performing the 
k integration for each order of k separately. The result is given by |6l)] 



Relic (iuj, p) 



off 



7ir 2 (3tu 2 - q 2 ) t4 248vr 4 (5u; 



<2>T 



45(Q 2 ) 5 



315(Q 2 ) 



2\4 



+ 



-T 6 +0(T 8 ), 
(2.4.2) 



ReHniiuJ, p) 



cff 



7^V , 2487r 4 g 2 (5u; 2 -g 2 



Q>>T 45(Q 2 ) 2 945(Q 2 ) 4 



-T 6 + 0(T 8 ). (2.4.3) 



The integration over the first term of the series (|2.4.1j) with the vacuum 
self-energy from l|2.2.9p then gives 



1 f d^Q 2n T +n L 77r 2 T 4 5 2 ff f d 4 Q 3w 



2 J (2tt) 4 Q 2 + n vac " 45 7 (2vr) 4 (Q 2 )3 1- 5e 2 ff ln(Q 2 /^ s )/127r 2 ' 

(2.4.4) 

By simple power counting we see that this integral is potentially logarith- 
micly divergent. However, if we perform the angular integral over the 4- 
sphere first, we see that the contribution f d^4(3uj 2 — q 2 ) =0 vanishes for 
any distance \Q\. Since the rest of 1)2.4.4(1 is a function of Q 2 only, any 
regularization or cutoff which respects Euclidean invariance will be UV well 
behaved. It is important to average over the Euclidean four spheres first, be- 
fore integrating over the radius \Q\. If we first perform one of the integrations 
uj or q first, and then integrate over the other variable without respecting 
Euclidean invariance, we might get spurious logarithmic divergencies. Of 
course, as long as we choose the upper integration limits such that in the 
end we integrate over the complete 4-sphere, nothing should happen (except 
for numerical cancellations to become more demanding). Integrating over a 
4-dimensional cylinder in the uj/q space, or cutting out a four- dimensional 
cube might bring in potentially logarithmic contributions for a large cutoff 
Q. But numerically it is easiest to integrate over the surface of the 4-sphere 
for a fixed Q, and then integrate over the radius up to the cutoff. 

Introducing a cutoff might cause gauge fixing dependence. For Lorentz 
and Coulomb gauges this turns out not to be the case because the self-energy 
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is gauge invariant as large Nf is basically an Abelian theory up to NLO. By 
the very same procedure as above, it can furthermore be shown that the 
introduction of an Euclidean invariant cutoff can be applied safely also to 
finite chemical potential. 

As in references [Ml EH we apply a cutoff and stop the d 4 Q integration 
at Q 2 = aA^, varying the value of a between 1/4 and 1/2 to estimate the 
irreducible ambiguity. 

In reference [HE], the "n&" terms from equation (|2.2.13|> were calculated 
in Minkowski space. Terms without n& were computed along a complex 
frequency contour which ran up the Minkowski axis to uJ m3jX < ALandauV^ 
for some a < 1, then along the great arc to Euclidean space, and back 
down to qo = \J q^ax ~ I 2 '■> finally, a Euclidean integration of the free 
term was performed over 4-spheres in Euclidean space up to Q 2 < AL andau a. 
In reference [HZ] a simpler way of calculating the integral was pursued: all 
pieces linear in were calculated in Minkowski space, and all terms without 
nf, were calculated in Euclidean space. By actually performing calculations 
in both ways, we had a rather non-trivial numerical check on the result 
which finally helped to reveal a coding error in the original result published 
|67| I66j . In the numerical implementation both ways turned out to agree 
within numerical errors of about 10 -5 . The term after the vanishing leading 
term in equation (|2~4~T|l is of the order (n|+n|) 2 /((5 2 + n vac ) 2 ~ T 8 so that 
the ambiguity introduced by the Landau pole is of the order 0(T 8 /AL andau ). 

2.5 Results and discussion 

In Fig. l2.3l we display our exact results 2 for the interaction pressure P — Po oc 
A^9, where the ideal-gas limit 



has been subtracted, for the entire fi-T plane (but reasonably below the scale 
Landau pole). For this we introduce an angle cj> = arctan ^ and encode the 
magnitudes T/Al and ///Al through the running coupling <7 2 ff (/iMs) with 
fills = + V 2 according to (|2.2.5|l . 

We found that the ambiguity arising from the presence of a Landau pole 
reaches the percent level for g 2 s ^ 28, where Al/ y ~k 2 T 2 + fi 2 S 19. At 
larger coupling (corresponding to larger T and/or fi), this ambiguity grows 
rapidly and will be shown in the two-dimensional plots below by a (tiny) red 
area. 

In the following we shall compare the exact large- Nf result with known 
results from perturbation theory at high temperature and small chemical 

2 Tabulated results are available on-line at |http : / /hep . itp . tuwien. ac . at/~ipp/data/| 




(2.5.1) 
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Figure 2.3: Exact result for the large-iVj interaction pressure P — P$ nor- 
malized to Ng(7r 2 T 2 + fj 2 ) 2 as a function of g% s (p,Ms) with p,^ = ir 2 T 2 + fj, 2 , 
which is the radial coordinate, and (b = arctan 

potential, where dimensional reduction is an effective organizing principle, 
and with results at zero temperature, where dimensional reduction does not 
apply. We also investigate to what extent quark number susceptibilities 
at vanishing chemical potential determine the behavior at larger chemical 
potential. 

2.5.1 Comparison to dimensional reduction 

Dimensional reduction is a method of expressing the static properties of a 
(3+l)-dimensional field theory at high temperature in terms of an effective 
field theory in 3 space dimensions j^EO] ■ In the imaginary time formalism, 
all nonstatic modes involve the temperature scale, so integrating them out 
to obtain an effective field theory on scales <C T leaves only the static modes 
with zero Matsubara frequency. The dimensionally reduced theory involves 
the color-electric screening scale gT and a dimensionful coupling g 2 T which 
is also the scale of color-magnetic screening. The perturbative result for the 
thermal pressure of hot gauge theories with fermions has been obtained by 
using this approach in reference (HSj at zero chemical potential to order g 5 . 
The analytic result of dimensional reduction through order g 5 with complete 
analytic dependence on arbitrary T and \i was calculated by jHS] and in the 
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large Nf limit this reduces to 



DR. 



-LT— rr, 3 + g eff - 

+T 12^ m£+ (48vr)2 r 



aE3 + 0G? c Yn (2-5.2) 



where the effective field theory parameter m 2 ^ is given by 



m E 
2^2 



1 + 



2^2 



7T 



ffeff 



(6tt) : 



2 In 



4^T 



1 - 



+ 0(«&) (2.5.3) 



and the coefficient of order g 4 is analytic in T and \i and given by 



OE3 



12 



2rp2 , ^ , 



7T- 



7T J 



In 



4vrT 



+4T 4 
,u 2 T 2 



1 16 C'(-3) 8 C 7 (-l) 

12 7 15 C(-3) 3 C(-l) 



7T- 



14 + 24 7 - 32 



C'(-i) 
C(-i) 



+ ^(43 + 36 7 ) 



7i 



-96T 4 {3H(3, 1) + 8N(3, z) + 3N(3, 2z) - 2N(1, *)} 
48i/iT 3 



96/T 2 T 2 



{N(0,z) - 12H(2,z) - 12K(2,2z)} 



7T 

with the definitions from 

Z E 

C'foy) = 
^(n, to) = 

N(io) E 



{4N(1, z) + 3N(1, 2z)} + U4 ^ 3r H(0, z) (2.5.4) 



7T° 



2 Z 27rT' 



C\-n,w) + (-l) n+1 C'(-n,w*), 
V(w) + *(to*). 



(2.5.5) 

(2.5.6) 
(2.5.7) 
(2.5.8) 



Here, n is assumed to be a non-negative integer and w a general complex 
number. ( is the Riemann zeta function and \& is the digamma function 
i$>(w) ee T' (w)/T(w). Note that despite the appearance of complex quantities 
in l|2.5.4p the coefficient «E3 is real as it has to be. 

For later usage we also present the effective field theory parameter of the 
coupling from dimensional reduction 



9 c sE = T 



2 , 5efT 
Sett + 



3(2^) 5 



2 ]n 

4vrT v ; 



+ 0(g 6 eS ). 



(2.5.9) 
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The expansion of the pressure in terms of small fj,/T is given by 
P-Po 



v + V^ 2 + V 

— Z ,77-^ 



7T- 



7T J 



32 



1 

12w" 



+ —Tm% + 



/20, 



24. 



2. .2 



In 



/%[S 

4vrT 



,1 88, 8C'(-3) 16 C'(-l) 

+ In 2 + 47 — H — - 

1 3 5 3 £(_ 3 ) ^ 3 C (_i) 



26 + 321n2- 247 r2 2 



+ 



12/i 4 



7T J 



7T 



hi— — +7+C4 



4vrT 



+ 



(48tt) 



+ 0( 5c 6 fr T 4 ), (2.5.10) 



where the terms oc g^ s and involving \i have been first obtained by Vuorinen 
|71| . The contribution to order g 5 arises from the NLO correction to the 
effective-field-theory parameter m E , which at finite \i was first computed in 
reference |69j . An alternative way of writing l|2.5.3p is given by 



T 2 \ 3 vr 2 T 2 



with the function 



<?eff <! 1 



6vr 2 



ln ?72^ + 2 P( ^r ) 



+ o(<?, 6 



eg; 

(2.5.11) 



£>(x) = -2 7 -41n2-2Re*(-(l + ia;)). 



(2.5.12) 



For small x this function can be expanded as 

n=l 



P(x) = 4 f)(-l)» (l - C(2n + 1) x 2 " , (2.5.13) 



with a radius of convergence of 1, which corresponds to fj, = irT. The 
only nonanalytic terms in <? 2 ff to the pressure in the dimensionally reduced 
effective field theory at large Nf stem from the "plasmon term" cx m E . As can 
be seen from the expansion in figure l2~T1 of the perturbative QCD pressure, 
the logarithms ln(g) do not survive the large Nf limit. 

In figure EU we plot the result obtained from dimensional reduction by 
Vuorinen |65j in the large Nf limit. In this plot we choose a naive renor- 
malization scale of /xms = y ( 7r ^) 2 + ^ ■ By comparison to figure l2~3l we see 
that the small coupling range is well reproduced for couplings g 2 ff < 9 for 
smaller \xjT but deviates considerably for larger fi/T. This is partly due to 
the choice of the renormalization scale which apparently is not the best. The 
exact result is independent of the choice of the initial renormalization scale. 
More precisely the scale dependence of the exact result enters only via the 
coupling g c s(p>Ms) and results can be easily converted from one renormal- 
ization scale to another by applying the renormalization scale dependence 
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Abbildung 2.4: Dimensional reduction (DR) result of the pressure (-Pdr — 
Po)/(Ng(TT 2 T 2 + /J?) 2 ) in the large Nf limit in the same range as in figure 
12 ..31 The renormalization scale is chosen as /2ms = f° r = 

/xms sin 4> and ji = /xms cos 4>. 
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relation (|2.2.4p . In the perturbative result this holds true only up to the 
order of perturbative accuracy. Numerically, due to truncation to a certain 
order in the coupling, the pressure not only depends on the renormalization 
scale indirectly via the coupling g e s {Pms), but also directly, since higher 
order contributions that would fix the renormalization scale dependence are 
not taken into account. The result will therefore in general differ whether we 
started from a renormalization scale p,\ and rescaled the resulting pressure 
according to l|2.2.4p to a renormalization scale fii — > p>2, or if we directly 
calculated the perturbative result with the latter renormalization scale p,2- 
Therefore, for perturbative results there is this additional question of which 
renormalization scale to choose. 

A possible optimization is provided by the FAC (fastest apparent conver- 
gence) scale, which is derived by demanding that the effective-field-theory 
parameter m 2 E from 1)2.5.11(1 vanishes to order g^ s and all higher orders. This 
means solving In 1/2^ T + 2^(7^) so ^ na ^ we obtain 



For /i = this expression reduces to //Ms(FAC-m) = 7rTe^~ yE while for 
T = it is given by /Sms (FAC-m) = 2 fie 2 . 

In figure 12.51 we see the difference between the pressure obtained from 
dimensional reduction to the exact large Nj result for this FAC-m renor- 
malization scale (HE]- As in the three-dimensional plots we vary temperature 
T and chemical potential \i depending on the angle 4> = arctan 2 ^. The 
result is given for values of the coupling < < 24 where the ambigu- 
ity introduced by the Landau pole is negligible numerically. In the plot we 
rescaled the dimensionally reduced pressure by the exact renormalization 
scale dependence l)2.2.4p to the scale /2ms = \J ( n T) 2 + ^ 2 to allow for uni- 
fied comparison. The lines enclose those areas where the difference is less 
than zblO -5 , ±10 -4 , or ±10 -3 in units of the normalized pressure. 1CP 3 is 
about the size of the NLO results, and thus marks the breakdown of the 
agreement between the two approaches (100% deviation). Notice that there 
is always a factor 10 of increase of the inaccuracy between neighboring lines 
of the contour plot. We have very good agreement for g 2 s < 7, but this 
value slowly decreases as the ratio fj,/T is increased. Only for small angles 
4> < 20° the behavior changes considerably and the ranges of good agreement 
shrink faster with decreasing fi/T. For the FAC-m scale dimensional reduc- 
tion gives lower values of pressure than the exact large Nf result for 4> ^ 25° 
and higher predictions for smaller fi/T. In between there is a transition from 
negative to positive differences which of course is of no real relevance to the 
error ranges in the plot. 

Another natural choice for scale might be to apply FAC to g\. This 



//ms (FAC-m) = vrTexp - 



-IE ~ 




(2.5.14) 
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Abbildung 2.5: Difference between the pressure obtained from dimensional 
reduction to the exact large Nf result for the FAC (fastest apparent conver- 
gence) renormalization scale /zMs(FAC-m) = 7rTexp(i — je — \^{^r)) f° r 
different values of <fi = arctan ^ and <? 2 ff , normalized to N g (7r 2 T 2 + /i 2 ) 2 . 
The result of the dimensionally reduced pressure is rescaled by the exact 
renormalization scale dependence (|2,2,4|) to the scale fiMS — y ( 7r ^) 2 + ^ 2 
to allow for unified comparison. 
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Abbildung 2.6: Same as figure l2~5l except that the scale is now FAC-g instead 
of FAC-m, that is applying fastest apparent convergence to g 2 E . Contrary to 
the FAC-m scale, the error is now always positive, meaning that dimensional 
reduction at this renormalization scale always overestimates the exact result. 
Apart from that the error of the FAC-g scale is comparable to the one of the 
FAC-m scale. 
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produces in the Nf — > oo limit 



MMs(FAC-g) = vrTexp - lE - -V(-^) 



(2.5.15) 



which curiously is just a factor l/\/e smaller than the FAC-m scale but shows 
otherwise exactly the same functional dependence on \i and T. 

Figure EU shows the difference of the pressures now for the FAC-g scale. 
It is conspicuous that this time there is no sign change in the plot. Indeed, 
contrary to what we saw in figure 12.51 for the FAC-m scale, the difference 
between Pqr and P cxa .ct = -Pnlo is now always positive, meaning that di- 
mensional reduction at this renormalization scale always overestimates the 
exact result. But apart from that, the absolute value of the errors of the 
FAC-g scale plot is comparable to the one of the FAC-m scale. One should 
not be irritated by the horizontal line in figure [231 at 4> « 25° where the 
range of small error is extended to larger values of the coupling g^ s : This 
merely indicates the change of sign of the error. 

Having such good agreement for a wide range of 4> and g^ s , it seems 
tempting to try to improve the result by averaging over the two scales: Tak- 
ing the arithmetic mean of the two scales FAC-m and FAC-g we obtain a scale 
that is a factor (1 + l/\/e)/2 ~ 0.803265.. smaller than the FAC-m scale. 
Remember that the FAC-m and FAC-g scale are directly proportional to 
each other with /2 M s (FAC-g) = p, M s (FAC-m) /y/e « 0.606531/2 M s (FAC-m) 
so that the arithmetic mean scale is also proportional to them. 

Figure 12.71 shows the remarkable agreement that follows from this choice 
of scale: The thermal pressure of dimensional reduction agrees to the exact 
large Nf result for <p > 65° over the whole range of couplings. If we look at 
smaller couplings g^ s < 10, the two results agree perfectly down to <f> > 35°. 
Only below cj) < 30° the spreading of errors does not differ much from that of 
the scales FAC-m and FAC-g separately. One should note however, that this 
might be a mere coincidence in the large Nf limit and that no real conclusion 
can be inferred from this result for finite, smaller Nf. After all, for finite Nf 
FAC-m and FAC-g are not necessarily proportional to each other (except for 
fj, = where both are simply proportional to the temperature T). 

2.5.2 Pressure at zero chemical potential 

At zero chemical potential we shall study in more detail the convergence 
properties and renormalization scale dependences of perturbation theory on 
the one hand, and the Landau pole ambiguity of our large Nf result on the 
other hand. In figiire l2~8l we give the numerical result of the thermal pressure 
as a function of g^^ifi = 7re _7B T). The perturbative results depend on the 
value of the renormalization point p, which we vary between irT and AirT. 
In order to compare the different perturbative results with the exact large 
Nf result, we use the exact running coupling 1)2 .2 .4(1 to rescale everything 
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9lM 

Abbildung 2.7: Same as figures 12,51 and 12,61 except that the scale is now 
chosen as the arithmetic mean of the two scales FAC-m and FAC-g. We see 
remarkable agreement for 4> > 65° over the whole range of couplings down 
to 4> > 35° for g^ s < 10. Below 4> < 30° the spreading of errors does not 
differ much from that of the scales FAC-m and FAC-g separately. 
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Figure 2.8: Exact result for -Pnlo /Pfree as a function of g^gifi = 7re~ 7B T), 
rendered with an abscissa linear in g e s, in comparison with the perturbative 
results through order g 4 and g 5 with a renormalization point chosen within a 
power of 2 of 2ttT. The line marked "FAC-m" corresponds to Jl = 7re2~ 7E T 
where the perturbative result to order g 4 coincides with the one to order 
g 5 . The line marked "FAC-g" corresponds to the renormalization point p, = 
■ne~ lE T and shows the result through order g 5 . 

to p, = ixe~~ m T . We see that the perturbative result through order g 5 gives 
reliable results only for effective couplings <^ ff < 4. If the perturbative result 
to order g 5 is optimized by fastest apparent convergence of m E (FAC-m) 
from equation (|2.5.14f) which for ^ = means the scale p, = 7re2~ 7B T, the 
agreement with perturbation theory is improved and extends to g^ s m 7. 
Figure EH shows for comparison also the FAC-g scale. 

Figure 12.91 shows the same thermal pressure as figure I2.8| but with the 
perturbative results varied between fife and fie where ji = 7re _7E T which is 
in fact the FAC-g scale. The central line of the g 5 result is just the FAC-g 
line from the previous plot. The central line of the g 4 result also corresponds 
to the FAC-g scale. Different from the FAC-m scale, the FAC-g scale stays 
not the same from order g 4 to order g 5 , but instead changes noticeably. 
Remember that the FAC-m scale was just chosen such that the g 4 and the g 5 
results lie on top of each other. Figure l2~9l also includes the arithmetic mean 
of the two FAC scales which is \ (FAC-m + FAC-g) = |(1 + l/y^e) FAC-m = 
\{\fe+ l)FAC-g. It turns out that this line is almost indistinguishable from 
the exact result. Only for larger values of g^ s can one see a tiny deviation 
from the exact result. 

The exact large Nf result flattens out for higher values of g^ s and reaches 
a minimum at g^ s ~ 12. For this coupling, AL an dau — 480T and the ambigu- 
ity introduced by the Landau singularity is completely negligible. For still 
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Figure 2.9: Same as figure l2~%l only with renormalization point chosen within 
a power of e of fj, = ne~' yE T. The FAC-g result is not explicitly marked, but 
coincides with the scale choice of fi in this plot. Therefore the central line 
of the # 5 result shows the FAC-g up to order g 5 , and the central line of 
the g 4 result shows the FAC-g result up to order g 4 . Note that unlike the 
FAC-m result, the FAC-g result differs considerably from order g 4 to order 
g 5 . Plotted in this graph is also the arithmetic mean of the two FAC scales 
as described in the text. This line is almost indistinguishable from the exact 
result. 




Figure 2.10: The result for PNLo/Pfree up to g^ = 35 and with different 
cutoffs. 
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Figure 2.11: Exact result for the interaction pressure at zero chemical po- 
tential as in references (Ml EI] but as a function of g^ s (ftMS = kT) or, al- 
ternatively, log 10 ( Al /vrT) . The dashed line is the perturbative result when 
the latter is evaluated with renormalization scale /2ms = P-fac = vre 1 / 2 ~ 7 T; 
the dash-dotted lines include the numerically determined coefficient to or- 
der g® s (with its estimated error) at the same renormalization scale. The 
result marked u g^ s = g® s " corresponds to choosing /2ms such that the order- 
ing coefficient vanishes and retaining all higher-order terms contained in the 
plasmon term oc m E . In this and the following plots the (tiny) red band 
appearing around the exact result at large coupling displays the effect of 
varying the cut-off from 50% to 70% of the Landau scale Al- 



higher values the pressure rises and starts to exceed the free theory pressure 
at g^ s > 28, as can be seen in figure l2~TT)l This occurs at a coupling for which 
ALandau/r < 34. While this still seems to be a reasonably large number, the 
numerical result starts to become sensitive to the cutoff just where the pres- 
sure approaches the free one. The four curves displayed in figure l2~T0l show 
the result of varying the parameter a in the UV cutoff v^-^Landau in the 
Minkowski and Euclidean parts of the calculation (clm and oe respectively) 
from a = 1/4 to a = 1/2. The numerical result is rather insensitive to this 
below <7g ff w 25, but very sensitive in the region where the pressure starts to 
exceed the free one. 

In figure l2~TT1 we show the pressure at \i = and compare with even higher 
orders of perturbation theory. Here the result is presented as a function 
of g^ s (fiMS = ^T) and not divided by the free NLO contribution from the 
pressure, but with the ideal-gas limit from equation (j2.5.1|) subtracted. In the 
plots we also give the corresponding values of the Landau pole log 10 ( Al/7tT) 
to different values of the effective coupling g^ s . This plot also contains the 
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result up to g®g which is not known analytically yet, but obtained from 
our numerical results at small coupling. For small coupling the agreement 
of our numerical results with perturbation theory is sufficiently accurate so 
that we can numerically extract coefficients to this order as follows: From 
the pressure in dimensional reduction (|2.5.1 Oj) we can infer by applying the 
exact running coupling of large Nj (|2.2.4p that the g^-term in the pressure 
at fi = has the form 



a 6 



iAt =o \ Air J 



(M)V[c 6 + f01n 2 ^| (2 .5.i 6) 



_16*V 1 + 127 _ 4641n2 g(-32 + 1 ^V/%s 
81 V 5 C(-3) C(-l)y vrT 

Notice the appearance of the (as yet) unknown constant Cq. For small values 
of <7g ff , perturbation theory is expected to give excellent agreement with the 
exact result, so that we can actually use our (numerically obtained) exact 
result to determine unknown coefficients in the perturbative series. By least- 
square fitting we obtain numerically the estimate Cq = +20(2). 

In real, finite- Nf QCD this result corresponds to the Njg 6 T 4 coefficient 
of the pressure. As we saw in figure I2.1| the purely gluonic contribution 
oc N®g 6 T 4 of the same coupling order g G is completely nonperturbative. 

In the figures 12.81 and 12.91 the agreement between exact result and per- 
turbative result up to order g^ s was limited to a range of about g^ s < 9. In 
ngure l2~TTl we show the improved result that remains accurate up to g^ s ~ 16 
by including our numerical estimate of the (^-coefficient and using /2fac- 
The agreement with the exact result can even be further improved by fixing 
the renormalization point such that the g® s coefficient vanishes and keeping 
all orders of the odd terms in g e g by leaving the plasmon term oc w? E unex- 
panded in g% s . The result of this procedure is indicated by the gray area in 
figure l2~TTl Obviously the result now gets very close to the exact result. Of 
course, this procedure does by no means guarantee that an analogous scheme 
for real QCD with finite Nf will produce similarly convincing results, but 
we can take it as another indication of the observation that keeping the pa- 
rameters of the dimensionally reduced theory unexpanded can improve the 
convergence of thermal perturbation theory |89j . 

Figure 12.121 shows a comparison of the ^-derivable 2-loop result in the 
HTL approximation (full lines) and in the next-to-leading approximation 
(full lines ending in dashed lines) as performed by Rebhan [84] . In the 
<l?-derivable approach of reference |46j . a fermionic gap equation had been 
assumed that turned out not to be compatible with the limit of large Nf. As 
is shown in |84j . the originally quadratic gap equation for the fermions does 
not comply with the Casimir scaling in the large Nf limit which would favor 
a linear dependence on the asymptotic mass Moo. The new gap equation 
(16) in reference (Hi] reproduces qualitatively the nonmonotonic behavior in 
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P/P (large- Nj limit) 






m D /T 

Abbildung 2.12: Comparison (taken from reference (SH) of the ^-derivable 
2-loop result in the HTL approximation (full lines) and in the next-to-leading 
approximation (full lines ending in dashed lines), for fx = T and 4ttT, and the 
exact result for the pressure in the limit of large Nf. The gray lines denote 
the next-to-leading approximation with quadratic fermionic gap equation 
considered in reference (Hj , but which Rebhan argues in reference [HI] needs 
to be replaced by equation (16) therein. The Debye mass appearing on 
the lower axis of the plot is related to the effective coupling via rhu/T = 
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Abbildung 2.13: Comparison (taken from reference (HD]) of the NLO HTLpt 
prediction for the O(Nj) contribution to the free energy, the exact numerical 
in the large Nf limit, and the perturbative prediction accurate to g 5 as a 
function of g c fr(^DR) = ^/sJg(jj, DR ) = 27r v / s / a s (/i DR ) where /x D R = vre~ 7 T. 
Dots indicate the point at which there is no longer a real- valued solution to 
the gap equation for mo. The renormalization scale \i is varied by a factor 
of e around /^dr and the perturbative g 5 result is evaluated at the central 
scale. 



5g ff of the exact result. By a curious coincidence, for N = 3 and Nf = 3 
the revised gap equation has exactly the same solutions as the uncoupled 
quadratic gap equations that have been previously in use. Only for Nf > 3 
there exists a coupling where the fermionic mass ceases to grow monotonicly 
with g. Because of this coincidence, the numerical changes in the previous 
results of are almost completely negligible. 

Figures f2. 131 and 12.141 show the HTLpt (Hard Thermal Loop - perturba- 
tion theory) comparison to the large Nf limit by Andersen et al. (HOI- The 
plots show the NLO HTLpt prediction for the O(Nj) contribution to the 
free energy, the perturbative prediction accurate to g 5 , and the exact result 
up to NLO in the large Nf limit. It should be noted, however, that both the 
HTLpt and the 2-loop ^-derivable results are perturbatively accurate only 
up to and including order g^ ff , where the perturbative result is rather ill- 
behaved. As noticed in (SJJ], the HTLpt predictions for both the free energy 
and the Debye mass (which was also discussed in [20]) seemed to diverge 
from the exact result around g c ff ~ 2 regardless of the scale that was chosen; 
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Abbildung 2.14: Same as figure l2~T3l only with the renormalization scale /i 
varied by a factor of 2 around 2-kT . The perturbative g 5 result is again 
evaluated at the central scale. [90] 



however, for both quantities, choosing the scale to be fx = [idr = 7re _7 T 
seemed to reasonably reproduce the exact results. The authors of (HD] con- 
cluded that their result is comparable to the performance of the ^-derivable 
approach in the large Nt limit in figure l2~T2l 

2.5.3 Quark number susceptibilities 
Linear quark number susceptibility 

The (linear) quark number susceptibility is defined as the first derivate of 
the quark number density N = N/V with respect to chemical potential (see 
also appendix IC.2p . 

dM d 2 P 

x = T^ = W (2 ' 5 ' 17) 



Figure l27T5l displa,vs the exact large- iV/ result for the interaction part of x 
at zero chemical potential as a function of g e R (or alternatively log 10 (AL/7rT)). 
Similar to the thermal pressure, the result is nonmonotonic, but the mini- 
mum already occurs at g^ s (7rT) 8.6, and the free-theory value is recovered 
at g^ ff (7rT) « 22.5, where the Landau ambiguity is still well under control 
since Al/T 100 at that coupling. 

The perturbative (dimensional reduction) result can be read from the 
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Figure 2.15: The interaction part of the quark number susceptibility at /i = 
compared with strict perturbation theory to order g^ s and g® s , respectively, 
with renormalization scale varied about ttT by a factor of 2. 



linear term in /x 2 of (|2.5.1()|) and gives 



X ~ Xo 
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n ^^T + l8 C(3) 



(2.5.18) 



The coefficient of g^ s has only recently been obtained in [JJ in a three-loop 
calculation. We can confirm its closed- form value by a numerical fit, which 
gives agreement with an accuracy of 2 x 10~ 4 , thus providing a good check 
on both our numerics and the analytical calculations of [2^. This level of 
accuracy allows us to also extract the order-g 6 term as (for /ims = KT) 



N g T* 



-4.55(9) x (M) 



(2.5.19) 



In Fig. 12.151 we show the perturbative results to order g^ s and g® s , vary- 
ing the renormalization scale about ttT by a factor of 2 (now without the 
improvement of keeping effective-theory parameters unexpanded) . The value 
Ams = is in fact close to //fac where it makes no difference whether m? E 
is kept unexpanded or not. We find that the quality of the perturbative 
result for the susceptibility is comparable to that observed in the pressure, 
with /2ms = n T being close to the optimal choice. 
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Figure 2.16: The interaction part of the higher-order quark number suscep- 
tibility d 2 (x — Xo) at /i = compared with strict perturbation theory to 
order g^ and g® s , respectively, with renormalization scale varied about ttT 
by a factor of 2. The coloured bands of the ^-results cover the estimated 
error of the numerically extracted perturbative coefficients. 



Higher-order quark susceptibility 

We also computed explicitly the higher-order susceptibility d 4 P/dfi 4 \^=o 
(which has been investigated in lattice QCD with Nf = 2 in reference |75J). 



Our exact result in the large- Nf limit is shown in figure 12.161 In this 
quantity, we find that the nonmonotonic behavior observed above in the 
pressure and the linear susceptibility is much more pronounced. The min- 
imum now occurs at g 2 s ~ 3.7, where perturbation theory is still in good 
shape, and the free-theory value is exceeded for g 2 s <J 9. Using (|2.5.10|) we 
find to order g 5 
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C(5) 



(2.5.20) 



+ 0(g 6 cS ). 



In the original publication of these results [HE] the coefficient appearing at 
order g 4 s has been numerically extracted as C4 = —7.02(3). In the meantime, 
the complete \x dependence of the dimensional reduction result to order g 4 s 
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has been worked out in reference |65j from where one can obtain the exact 
result 



The complete agreement with reference [EE] provided on the one hand an 
independent check on the correctness of the 3-loop calculations of [HE! and 
on the other hand a check on the accuracy of our numerical analysis. 

Using (|2.5.21j) we can extract the term of order g® s in equation (|2,5,2()|) 
as — 39(l)^ ff /(1287r 8 ) for = irT. The perturbative results to order g^ s and 
to order g® s are compared with the exact result in figure 12.161 This shows 
that the accuracy of the perturbative result again improves by going from 
order g^ s to order g® s , but the renormalization scale dependence increases 
sharply at large coupling. 

Pressure at larger chemical potential from susceptibilities 

With regard to the attempts to explore QCD at finite chemical potential by 
means of lattice gauge theory [23 EH1 HH HEj ; it is of interest how well the 
pressure at larger chemical potential can be approximated by the first few 
terms of a Taylor series in . 

In reference (ZH] it has been observed that the ratio of AP = P(T, fj.) — 
P(T, fj, = 0) over the corresponding free-theory quantity APo is practically 
independent of [i for the range of chemical potentials explored. This is 
also realized when quasi-particle models are used for a phenomenological 
extrapolation of lattice data in a method introduced by Peshier 

et al. [HS]- In figure 12.171 we show the deviation from this "scaling" at 
higher values of /z/T by considering the quantity 5P = P(T, fj.) — P(T, 0) — 
±xU=o(/^ 2 + Ai 4 /(2vr 2 T 2 )) divided by P^ e Q = N g ir 2 T A /4b. The combination 
(fj? + fj, 4 / (2tt 2 T 2 )) appearing therein is such that a replacement of P and 
X by their interaction-free values Po and xo makes 5P vanish identically. 
(As can be seen from the above perturbative results, 5P also vanishes for 
the leading-order interaction parts oc g 2 s ) In the exact large- Nf results of 
figure \TT7\ we observe that for coupling g 2 s < 4 the deviation 5P is at most 
a few percent of -P^lo f° r ~ n > ^ rapidly grows for fi/T > tt. This 
is in fact also nicely illustrated by the 3-dimensional plot of the pressure in 
figure 12. 3| which has a rather conspicuous kink at 4> = 45° corresponding to 
p = irT. 

It turns out that including the exact higher-order susceptibility at fi = 
does not lead to a better approximation of the pressure at larger chemical 
potential. The dashed lines in Fig. 12.171 correspond to 5P = P(T, fi) — 
P(T,0) — 4rxU=0 — ^T7^U=0- While this slightly improves matters at 
small fi/T, it results into even quicker deviations for larger fi/T. 

It is of course impossible to say whether this behavior would also appear 
in real QCD, but since it occurs already at comparatively small g e s in the 
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= -6.9986997796998 . . . 



(2.5.21) 
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Figure 2.17: Deviation from the scaling observed in Ref. (ZH] in lattice QCD 
for small chemical potential in the quantity 5P = P(T, /i) — P(T, 0) — 
|xU=o(^ 2 + V A /(2tt 2 T 2 )) (full lines) and in 5P = P(T,n) - P(T,0) - 
^|-xU=o — 41"^? 1/^=0 (dashed lines), both normalized to -Pj^LO = A r 9 7r 2 T 4 /45, 
for g 2 cS (nT) ="1,4,9,16. 

large- Nf limit, where the peculiar nonmonotonic behavior of the pressure as a 
function of g e R does not yet arise (the minimum in the normalized interaction 
pressure occurs at g 2 s (nT) m 14), it may be taken as an indication that 
extrapolations of lattice data on the equation of state from small chemical 
potential to large fi/T are generally problematic. If anything, real QCD 
should be more complicated because of the existence of phase transitions 
which are absent at NLO in the large- iVj limit. 

2.5.4 Pressure at zero temperature 

Our exact result for the thermal pressure at zero temperature and finite 
chemical potential is given in figure l2~T8l as a function of g 2 s (jlMS — I 1 )- ^ n 
contrast to the pressure at zero chemical potential and finite temperature, 
the interaction pressure divided by /i 4 is monotonically decreasing essentially 
all the way up to the point where the Landau ambiguity becomes noticeable. 

The thermal pressure at zero temperature and large chemical potential 
for QED and QCD has been obtained to order g long ago by Freedman and 
McLerran (221 EH|- At this order, there is a non-analytic zero-temperature 
plasmon term cx <? 4 ln(g), whose prefactor is known exactly, but the constant 
under the logarithm only numerically. The transposition of their result, 
which has been obtained in a particular momentum-subtraction scheme, to 
the gauge-independent MS scheme can be found in references ^SlEE]- The 
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Figure 2.18: The interaction part of the pressure at zero temperature and 
finite chemical potential as a function of g^(p-MS — m) or > alternatively, 
log 10 (AL//x), compared with the perturbative result of Freedman and McLer- 
ran [22 ESI to order g* s , and our numerically extracted order-^g result, both 
with renormalization scale in the perturbative results varied around p>MS = V 
by a factor of 2. The coloured bands of the ^-results cover the error of the 
numerically extracted perturbative coefficients. 



large- Nf limit of this result reads 
P-P 
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(2.5.22) 

and involves one of the numerical constants computed in reference [72] . 



79 7T 2 7 log (2) 



C 4 = 

18 



2b 

— « 0.536, 



(2.5.23) 



where b has an integral representation, given in equation (II.3.25) 3 of ref- 
erence (Z2], that apparently cannot be evaluated in closed form. With bet- 
ter computer equipment, b can however be evaluated numerically to higher 
accuracy than that given in [22] as b = —1.581231511..., which leads to 
C A = 0.5358316747 .... 

The accuracy of our numerical results is sufficiently high to confirm the 
correctness of the result for C4 with an accuracy of ~ 2 x 10 -4 . With the 
knowledge of the exact value of C4 we can also extract, with lower precision, 



3 Note that there is a typo in the equation mentioned: Comparison with equation 
(II. 3. 24) shows that there is a missing exponent 2 after the second set of large round 
parenthesis in equation (II. 3. 25) of as pointed out by |t)8| . 
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the next coefficients at order g^ s , which again involve a logarithmic term: 
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In figure 12.181 we also study the renormalization scale dependence and 
apparent convergence of the perturbative result. We have varied /zms about 
fj, by a factor of 2, and it emerges that the larger values are somewhat favored. 

At low temperature T fx, dimensional reduction does not occur. If one 
nevertheless considers the effective-field-theory parameter m? E of (|2.5. 1 1|) in 
this limit, one finds that the function V{x) therein approaches —2 (]n2x+ r y), 
so that the T 



limit of m 2 E exists and reads 



2 2^ff /-, S%K 
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2/i 



+ 0(g, 



cS) 



(2.5.25) 



Fastest apparent convergence applied to this quantity would suggest /2ms = 
2ea/t « 3.3/x. This turns out to be not as good as the choice of 2/i, though 
slightly better than /2ms = M- 

The region of low temperature T <C /t might be readily explored in the 
large limit, where it gives rise to anomalous contributions of the pressure 
in the context of non-Fermi liquid. We will study this region in the next 
chapter. 
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Kapitel 3 

Non-Fermi Liquid 



3.1 Introduction 

3.1.1 From ideal gas to non- Fermi liquid 

What is a non-Fermi liquid? The straightforward answer is to start from the 
ideal gas and extend the model layer by layer: The ideal gas is a statistical 
description of non-interacting point-like particles. If we go to low temper- 
atures, we cannot neglect quantum mechanics and we have the choice of 
describing either bosons, which leads to the Bose-gas and to Bose-Einstein 
condensation at very low temperatures, or fermions, called (ideal) Fermi 
gas. Still, the Fermi gas is an accumulation of non-interacting fermions. If 
we turn on short-range interaction we get to the description of a (Landau-) 
Fermi-liquid. Thermodynamic quantities will mostly show the same order 
of dependence on the temperature T as in a Fermi gas, but their magnitude 
might change drastically. Finally, if the interactions we introduce change the 
dependence on T qualitatively, which is the case for long-range quasi-static 
interactions, the system cannot be described by the quasi-particle picture of 
the Landau-Fermi-liquid theory, and is therefore called a non-Fermi-liquid. 
Thus, a non-Fermi liquid is a cold gas/liquid of fermions with long-range in- 
teractions that change the behavior of a Landau-Fermi liquid qualitatively. 

The ideal gas of not too dense, not too cold, and "non-interacting" (ac- 
tually only interacting by elastic impact for thermalization) particles obeys 
the ideal gas law pV = NksT. While simple and universal, this law is not 
sufficient to calculate all thermodynamic quantities like energy or entropy of 
the system. To calculate the energy for a fixed number of particles E(S, V), 
we need an additional quantity, for example the specific heat Cy = T 
which depends on the inner degrees of the ideal gas. Each degree of freedom 
in every particle of a system will contain the same energy on average so that 
the specific heat is given by 

C V (T) = ^Nk B (ideal gas) (3.1.1) 
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for / degrees of freedom per particle and the Boltzmann constant k B . For 
pointlike particles in three dimensions / = 3, for molecules of two atoms 
for example there is one additional degree of freedom for rotation and one 
for the change of bond length so that / = 5. The specific heat is an ex- 
tensive quantity, but we can divide by the volume Cy = Cy /V = ^nk B 
which now depends on the particle density n = N/V . Note that the specific 
heat for an ideal gas is independent of temperature. The entropy of the 
system can be calculated for example by integrating (fy)y = Cy /T and the 

Maxwell relation (§y) T = [J}tJ = Nk B /V which for the ideal gas gives 

S(T, V) = Cy ln(T/T ) + Nk B ln(V/V&) + S(T , V ). Statistical distribution 
of quantities like the mean velocity is given by the Maxwell-Boltzmann dis- 
tribution function, which is not valid anymore for small temperatures where 
quantum mechanics starts to play a role. Here we have to use Bose-Einstein 
statistics for bosons or Fermi-Dirac statistics for fermions. 

The ideal Fermi gas is a description of cold, non-interacting fermions 
OlESJ. Examples include electrons in metals and semiconductors (when 
the Coulomb interaction between them is neglected), as well as neutrons 
in a neutron star (again when neglecting the interaction). Fermions obey 
the Pauli principle: no two fermions with the same quantum numbers can 
be in the same energy state. As a consequence, fermions will fill up the 
Fermi sphere in momentum space - if one possible momentum is occupied 
by a fermion, the next fermion has to go to the next higher momentum. (If 
they differ by a quantum number, e.g. the spin, then they might occupy 
the same energy level.) The particle density corresponding to a completely 
filled Fermi sphere of radius pp in momentum space for a spin 1/2 fermion 
gas is given by n = N/V = p F / '(3tt 2 h?) at zero temperature. The number- 
density distribution for T = is a step-function: 1 for energy states below 
the Fermi energy, and for energy states above. For non-zero temperature, 
the number- density distribution is given by the Fermi-Dirac distribution n = 
(g^-cl/fo 7 ) _j_ 1) _1 5 where the chemical potential is just the energy on the 
Fermi surface at zero temperature h\t=o = £f = e(pf)- The specific heat 
can be calculated from the energy change due to the Fermi-Dirac distribution 
at non-zero temperature. It is given by 

7T 2 

Cy{T) ta —k 2 B D{e F )T + <3(T 3 ) (ideal Fermi gas) 
= V^k* B T + 0(T*) 

= yiV^-^ + 0(T 3 ) (3.1.2) 

with a density of states factor D(e) = ■^F(m//i 2 ) 3 / 2 e 1 / 2 and the Fermi 
energy e_p = t^(3tt 2 N/V) 2 / 3 = p_p/(2?n) = k B T F . We see that the specific 
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heat of the ideal Fermi gas grows linearly in T for small temperatures. This 
is in contrast to the specific heat of an ideal gas (|3.1.1j) that is independent 
of temperature. Another striking difference to the ideal gas is the pressure: 
Since the fermions fill up the Fermi sphere, only a small fraction of fermions 
can be in the ground state. As a consequence, the pressure of a Fermi gas 
is nonzero even at zero temperature. This is in contrast to the pressure of 
an ideal gas that would be zero (pV = NksT). This so-called degeneracy 
pressure stabilizes a neutron star (Fermi gas of neutrons) or a White Dwarf 
star (Fermi gas of electrons) against the inward pull of gravity. 

If we include (short-range) interactions between the fermions, we call 
the system a Fermi liquid. A description of Fermi liquid was first given in 
1956 by Landau [HEIHEI and is thus commonly dubbed Landau- Fermi-liquid 
theory. It can be applied to liquid |He (which contains 2 protons, 1 neu- 
tron and 2 electrons - an odd number of fermions per atom - such that the 
atom itself is a fermion 1 ), as well as to the electrons in a normal metal. 
Landau-Fermi-liquid theory is qualitatively very similar to the theory of an 
ideal Fermi gas in which quasi-particles take the role of the non-interacting 
particles of the Fermi gas. Quasi-particles are collective excitations of the 
macroscopic system with certain energy e and momentum p. The properties 
of quasiparticles are mainly characterized by the dispersion relation e(p). 
Quasi-particles might have finite lifetime which gets sufficiently long in the 
vicinity of the Fermi surface to allow for a description similar to the par- 
ticles in a Fermi gas. Strictly speaking, only certain kinds of interaction 
lead to the Landau-Fermi liquid, for example quasi-particles should still be 
fermions. (If there are attractive forces that favor a pairing of the fermions 
into bosons, the description of the Landau-Fermi-liquid is not applicable 
anymore.) Since the quasi-particles are fermions, they also fill up the Fermi 
sphere in momentum space. Quasi-particles always have a spin 1/2 spec- 
trum. If the underlying particles had spins different from 1/2, this would 
lead to a degeneracy of the energies of the quasi-particles so that each branch 
corresponds to a spin 1/2 quasiparticle. The velocity of a quasiparticle is 
defined as vp = (de/dp)\ p=PF . In the non-relativistic case one can introduce 
the effective mass of the quasiparticle as m* = pp/vf- With this definition, 
the specific heat can be derived as 

C V (T) » V^^k 2 B T + <3(T 3 InT). (Landau-Fermi liquid) (3.1.3) 

The only difference to the formula (J3.1.2|) in the leading order contribution 
is that the particle mass m is replaced by the effective mass m* of the 

1 Note that liquid ^He with an even number of fermions per atom forms a boson gas and 
is a superfluid below 2.17 K. At this temperature, ^He behaves like a Landau-Fermi liquid. 
Only at much smaller temperatures, around 2.7 mK, two atoms of fHe align themselves 
to give an overall spin s = 1 and angular momentum I = 1 and also form a superfluid. 
This effect was found in 1971 by Osheroff, Lee, and Richardson who in 1996 received the 
Nobel prize for their work. 
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quasi-particle. We still have the same linear dependence on temperature 2 
|98j . Note also that the next order contribution changes from 0(T 3 ) for 
the Fermi gas to 0(T 3 In T) for the Landau-Fermi liquid |99j . Not only the 
specific heat, but also other thermodynamic quantities like compressibility 
or spin-susceptibility show the same qualitative temperature dependence as 
in the Fermi gas, but might have different coefficients. The classical Landau 
Fermi liquid theory was extended to relativistic Fermi systems, enabling the 
study of high density matter with weak interaction via scalar and vector 
meson exchange [100] . As mentioned, not every kind of interaction leads to 
a Landau-Fermi liquid: If the interaction favors a pairing of the fermions, 
we effectively deal with a Bose gas. (In fact, in a Fermi gas arbitrarily small 
attractive forces lead to Cooper pairing of fermions.) If we have long-range, 
quasi-static interactions, the quasi-particle description is also not valid any 
more, and we obtain non-Fermi-liquid behavior. 

Anomalous contributions to the specific heat were first calculated by 
Holstein et al. |101j who found a deviation from the linear dependence on 
T by a term proportional to TlnT -1 . For very small temperatures T, the 
temperature dependence will actually be dominated by the TlnT -1 term. 
For QCD, quasi-static transverse gauge boson interactions further lead to 
a series with fractional powers of T . The specific heat of this non-Fermi 
liquid is given by (now with the usual quantum-mechanical natural units of 

k B = n=i) 

C V = C V /V ps [j%T——£- + ..Tln..T _1 (non-Fermi liquid) 

+..T 5 / 3 + ..T 7 / 3 + 0(T 3 InT) (3.1.4) 

where the coefficients were first calculated in 2003 |102j and will be presented 
in the following sections. 

3.1,2 Non-Fermi liquid 

A "non-Fermi liquid" is in principle any thermodynamic system of fermions 
that goes beyond the classical Landau-Fermi liquid description (including 
for example Cooper pairing of fermions and color superconductivity). Here 
we will narrow its usage to a Landau-Fermi liquid that is changed quali- 
tatively by the introduction of long-range interactions. It was noticed by 
Holstein et al. |1 1| in 1973 when studying the de Haas-van Alphen effect 3 

2 For a Bose liquid the specific heat grows proportional to the cube of the temperature 
Cv = V r 2 7 r 2 r 3 /(15(ftu) 3 ) with the sound velocity u = e/p El 

3 The de Haas- van Alphen effect states that magnetization shows an oscillatory depen- 
dence on the inverse magnetic field if the magnetic field is strong enough, T < [ibB <C /i 
with the Bohr magneton /ib = eh/(2m e ) times the magnetic induction B still smaller than 
the chemical potential |07l llO.'-ij . A determination of oscillations in the magnetization 
M as a function of the inverse magnetic field 1/H (the de Haas- van Alphen effect) for 
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that long-range interactions lead to a deviation from the linear temperature 
dependence of the specific heat for small T. In order to explain the special de 
Haas-van Alphen oscillations in some noble metals in a slight modification 
of the original model, a current-current interaction, which in contrast to the 
Coulomb potential is not screened at zero frequency, leads to an anomalous 
contribution to the specific heat of the form TlnT -1 . The coefficient of this 
term is of the order ~ ol{vf /c)(m* /m) 2 ~ 10~ 5 smaller than that of the 
dominant term linear in T for non-relativistic QED applications, so origi- 
nally it was argued that this anomalous contribution cannot be the cause for 
measurable effects. (A factor of 4 that was missing in the original calculation 
was corrected by |104j so that the leading log correction to the specific heat 
is Cy = {gQp 2 F /?>&-K 2 )T\riT~ l ). While this contribution implies that the en- 
tropy of the non-Fermi liquid exceeds the entropy of a Landau-Fermi liquid 
below some temperature, it is not possible to infer at which temperature this 
happens without complete knowledge of the argument of the logarithm: As 
long as we take the logarithm of a dimensionful quantity In T, we are missing 
the information on the relevant scale for the appearance of the anomalous 
behavior. 

Long-range interactions are generally screened in the presence of a large 
Fermi sea. The reason why this does not happen for transverse gauge bosons 
like photons is that gauge invariance prevents them from acquiring a mass, 
unless gauge invariance is spontaneously broken (as in superconductors). An 
electron gas interacting via transverse gauge bosons shows non-Fermi liquid 
behavior depending on the number of space dimensions D. The specific 
heat is proportional to Cy ~ TlnT -1 only in D = 3 dimensions. For 
D < 3 dimensions, Cy ~ T D / 3 while for D > 3 the system behaves like an 
ordinary Landau- Fermi liquid |l()4(ll()5] . In this sense D = 3 corresponds to 
a quantum critical point with respect to variation of the dimension. 

A recent calculation of non-Fermi-liquid corrections using renormaliza- 
tion group resummation techniques in reference |106j suggests a different 
leading nonanalytic behavior of the specific heat proportional to T 3 lnT, 
which actually is of the same order as ordinary Landau-Fermi liquid correc- 
tions coming for example from electron-phonon interactions [HHl 11071 11081 
ll()9| . However, the starting point of reference |1()6| already neglected dia- 
grammatic contributions which would give TlnT -1 corrections. 

It is sometimes assumed that the anomalous contributions to entropy 
or specific heat will hardly play any role at all, since color superconductiv- 
ity (CSC) will dominate dense quark matter long before non-Fermi liquid 
behavior becomes effective |110| II 11] . The argument is based on simple di- 
mensional analysis of the energy scale, and compares the critical temperature 
of color superconductivity to the scale where the leading log correction be- 



different orientations of the field has been the most successful method for mapping out 
the shape of the Fermi surface of metals [lO.'-ij . 
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comes as large as the leading contribution, that is where non-Fermi liquid 
becomes nonperturbative. Using the known formulae for the CSC critical 
temperature and our complete leading logarithm contribution presented in 
the following, we can calculate the perturbative correction stemming from 
non- Fermi liquid effects. As we will see in figure l3~6l of section l3.6.2| pertur- 
bative NFL corrections to the specific heat are of the order of 10%-20% at 
the CSC critical temperature for g ~ 0(1) and of the order of g/(3\/2) for 
small g. Moreover, quark components which do not participate in the for- 
mation of diquark condensates will produce non-Fermi-liquid behavior even 
in the color superconducting phase. 

Non-Fermi liquid behavior may therefore occur in astrophysical situa- 
tions, for example the cooling rate of proto-neutron stars |16j if they involve 
a normal (non-superconducting) degenerate quark matter component. 

In solid state physics experiments, non-Fermi-liquid behavior has been 
measured in so-called heavy- fermion metals |112| . for example in the specific 
heat of the recently analyzed YbRli2Si2 crystal There, TlnT" 1 has 

been demonstrated experimentally over more than an order of magnitude 
in the vicinity of the quantum critical point for non-Fermi liquid behavior. 
At ambient pressure there are only a few undoped compounds that show 
non-Fermi-liquid behavior, like UBei3 CeNi2Ge2 |114| . or CeCu2Si2 

|115j . Other heavy fermion metals can be tuned to a quantum critical point 
by varying for example doping, pressure, or magnetic field. 

In the following we will show how the specific heat for ultradegenerate 
QED and QCD can be calculated from the pressure integral we used for 
calculating the large Nf limit. Non-analytic contributions in T come from 
transverse gauge boson contributions. We will complete the leading logarith- 
mic contribution TlnT -1 in that we calculate the argument of the logarithm. 
Beyond this contribution, dynamical screening gives rise to anomalous frac- 
tional powers 7 1 ( 3 + 2n )/ 3 i f or which we calculate the coefficients up to and 
including order T 7 / 3 . 

3.2 Entropy at small temperatures 

We want to study thermodynamic quantities in the region of small temper- 
ature T <C (J,. The anomalous contributions that we would like to calculate, 
can be located in the entropy density 

from which the specific heat can be derived. One might expect that the 
"rib" contributions containing the Bose-Einstein distribution in the pres- 
sure H2.2.13J) are negligible compared to the "non-raj," contributions, because 
nb{u}) = l/(e w / T — 1) will vanish exponentially with small T for each (fixed) 
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<jj. However, we will see that "rift" contributions actually cannot be neglected 
for a region with uj < T and that they are the source of anomalous terms in 
the low-temperature series. 

Looking at the "non-rift" contributions first (which are those parts of the 
calculation that we have to integrate in Euclidean space and which require 
the introduction of a cutoff due to the Landau pole), the corresponding part 
of the entropy is given by 



We see that the coefficient at order g^ s is the same as the g^ part of the 
strictly perturbative result for the pressure. The latter is also known as the 
exchange term j2H] and coincides with the g^ s part of 1)2.5.101) . Here we have 
also numerically extracted the order-# 4 ln(g) corrections from the exact large 
Nf result of S Q0Q . nb at small coupling. 

The dash-dotted lines in figure 13.11 show the exact (NLO large Nf) re- 
sult for S QOn - nb for different couplings <?g ff (/%[S = m) = 1> 4, and 9 and small 
temperatures < T/fj, < 0.15. In this range of temperatures, S non . nb is well 
approximated by the linear term 1)3.2.2)1 . The dashed lines in figure 13.11 show 
the "n&" contributions to the specific heat in the large Nf limit evaluated 
numerically. For small T we see that these parts cannot be neglected. For 
sufficiently small T/fj,, S nb is positive and even dominates so that the total 
result for the entropy turns out to exceed its free-theory value for a certain 
range of T/fj,. The size of the range where the the anomalous contribu- 
tion dominates in the entropy gets larger with increasing g^ s . The largest 
part of the positive and nonlinear contributions at small T j [i in fact comes 
from the transverse vector-boson modes. These are only weakly dynamically 
screened at small frequencies and completely unscreened in the static limit 
because of gauge invariance. We will see that it is the transverse part of the 
"non-nb" term that will fully contain the temperature-dependent anomalous 
contribution. The large Nf limit will allow us to calculate the anomalous 
contributions in a very straightforward manner, so let us start by looking 
more closely at the transverse contribution. 

3.3 Transverse contribution 

We started our large Nf calculation from equation ()2.2.13p which gives the 
pressure associated with the gauge boson loop with a resummed fermion loop 
insertion. All interesting non- Fermi liquid (TlnT, T 5 / 3 , and T 7 / 3 ) behavior 
to the specific heat arise from the thermal part, the "non-rtft" part, of the 




(3.2.2) 



0.328(1) x In |4 + 0.462(5) + 0(g 6 cS lng cS ) + 0(T 3 ). 
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s-s< 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 ^ 



Abbildung 3.1: The interaction part of the entropy at small T/fj, for 
<?efT (Ams = A*) = 1, 4, and 9. The "non-?V contributions (dash-dotted lines) 
are negative and approximately linear in T with a coefficient agreeing with 
the exchange term cx g^ s in the pressure at small coupling; the "ri&" contri- 
butions (dashed lines) are positive and nonlinear in T such that the total 
entropy exceeds the free-theory value at sufficiently small T/fx. Transverse 
gauge boson modes dominate the anomalous contribution in the "non-nb" 
part . 
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transverse propagator. In the following, we will calculate the contribution 
to the pressure coming from the transverse part 

= -/ -^l°°^2n b IrnHQ 2 -1 2 o+^T + U vac ). (3.3.1) 

We will start with a motivating calculation for the leading logarithmic con- 
tribution to the pressure and systematically extend this calculation to finally 

obtain the anomalous Fermi liquid contributions up to and including order 
T 7/3_ 



3.3.1 Leading log contribution 

To calculate the pressure at small T, we start from equation (|3.3.ip where we 
can write the argument function as Imlnrr = arctan(Im x/Re x) as long as 
Re x is positive. Keeping this in mind, we use the following low-temperature 
approximations: We first neglect the real part of Re(Hr + n vac ) ~ + 0(g 2 ) 
so that the main contribution comes from q 2 . For the imaginary part we only 
go through order 0(q) (the complete results for real and imaginary parts of 
TLt are given in appendix El s ° that we use the following expression 

lining - <Jo + n T + II vac ) ~ arctan — . (3.3.2) 

We first apply the g-integration of the form 

dq q 2 arctan ^ = - ( nq'L — 2q"L arctan — + go In ( 1 + 

q 6 6 \ qo V % 

ig (l + ln^+O(g 3 ) (3.3.3) 
3 V Qo J 

for a maximum upper bound q m . From the imaginary part we see that 
q m = 2/i. We are left with the following qo integration which gives 



PT,n b 


4tt f°° 


dq 




8^ Jo 


7T 






Mn 




72vr 2 






[in 




72n 2 



327ru 6 
~2~^ + IE ~ - 

9 c s T n 



+ 7fi-Z2C'( 2 )) (3-3.4) 



with 7£ — %('(2) ~ 1.14718. One does not have to worry about first ex- 
panding in small qo and then integrating qo from to oo: The bosonic 
distribution function ra& = l/(exp((/o/T) — 1) ensures that for small T only 
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small values of qo < T are sampled in the integrand - larger qo contribu- 
tions are suppressed exponentially. We find the contribution to the entropy 
density S = S/V = {8P/dT) v ^ as 

S T ,n, _ g 2 c ^ 2 T / 327^ 1 6 



N g 36vr 2 



We can also form the specific heat Cy = Cy /V at constant volume and 
number density which for unit-volume is given by |116j 



*- T \{w)A")X*) r r (3 ' 36) 

where J\f is the number density M = N/V = (dP/d/xjTy- The application 
of this formula including the leading-order contribution is given in equation 
(|C,1,5|) from appendix O We find that to leading order the coefficient of the 
logarithm remains the same for entropy and specific heat. There is only a 
shift by —1 in the sublogarithmic term 

Cv,T,n b _ 9 2 e ^ 2 T ( 32tth 



1 + 7*- J»C(2)). ( 3 " 3 - 7 ) 



N g 36vr 2 \~g*sT 

This first version of the result already includes the correct pre-factor of the 
leading logarithmic contribution of the specific heat as given in the literature. 
It confirms the calculation of reference |1()4| who also found the leading log 
contribution Cy = {g^p 2 ? / ?>&it 2 )T \riT~ 1 as the specific heat correction, and 
even goes beyond this known result in that dimensional dependence and the 
constants "under" the logarithm are determined. But the term proportional 
to T (the argument of the logarithm) still calls for improved approximations 
and should not be trusted yet: It will change by inclusion of higher order 
terms as we will see in the next section. 



3.3.2 Straightforward improvements 

It is possible to retain another term and replace A/j, 2 by (4/i 2 + q 2 ) in (|3.3.2jl . 
This will enable us to obtain the correct leading contribution to the q- 
integration even for larger q. The integral over the arc tangent involves 
roots of some cubic equation, but it is still solvable. After performing the 
small-go expansion, the result is almost the same, up to an additional term. 

(3.3.8) 

Using q m = 2//, this term only contributes by a shift of +3/2 inside the 
braces. The second integral is straightforward and leads to the following 
pressure 



2 7T Z J 



N g 727T 2 V 9 2 cS T 
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and its corresponding entropy 

%f^('^ + '+- ">))• (3-3-10) 

It is also possible to add the vacuum term which also contains a part 
proportional to q 2 : Instead of q 2 in the denominator of (|3.3.2|) we could 
start with 

q 2 - q 2 + II T + II vac ~ q 2 + In ^ + O(, 2 ) + 0(g 3 ). (3.3.11) 

The calculation will be almost unaffected, only in the final step we need to 

2 _ 

change g 2 s — > g 2 s / (1 + In It is clear that this only gives corrections 

subleading in g 2 s . In fact, this substitution is nothing else than a change in 
scale as given by the exact scale dependence of g e s in (|2.2.4|) . 

Another extension that we should consider is the T-dependence of the 
self- energies. Up to now we tacitly assumed constant self-energies Ht(T) = 
IIt(O) at small T, but this is a simplification which might affect the next- 
to-leading order in temperature. We find the following temperature depen- 
dencies: 

2 

lmU T (q ,q,T) = lmU T (q , q, 0) - T 2 ^q 9{2(i - q) + O(gg), (3.3.12) 

2 2 

ReU T (q ,q,T) = ReU T (q , q, 0) + T 2 ^^ + O(q 2 ) + 0(q 3 ). (3.3.13) 

The leading order contributions turn out to have the same order as the T = 
terms, so that our calculation does not change qualitatively. Introducing new 
variables 

TT 2 T 2 

X = 1 + ^r^, (3.3.14) 



3// 

Y = l + #?"||ln^ (3-3.15) 
36/^ 67H /XMS 



we can write the resulting pressure as 



Pr,n b gj^T 2 ( 327T/1 Y 6 3U 

= r ~9^T X 1E ~^ Q {) + 2Xj T (3 " 3 - 16) 

effectively adding T 4 contributions to the pressure (or equivalently T 3 con- 
tributions to the entropy or the specific heat). This expression reduces to 
(J3~3~5|) if we set T -> in X and Y. 

We should also consider consistently including higher order terms of q 
and qo in the arctan arguments of (|3.3.2p or l|3.3.8|) . In the next section 
we will see how this will eventually lead to an integral with a tenth degree 
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Abbildung 3.2: Integrand for the (/-integration (/ 2 Imln((/ 2 — q$ + 11^ + II vac ) 
for /i = /2ms/2 = 1, qo = 0.2, g 2 s = 1. The solid line shows the exact 
result that follows from the full T = self energy expressions, the dashed 
line shows the result with Re Hp expanded through order O((?o) and 0{q 2 ) 
and 0(2// — \q±qo\) replaced by 9(2/1 — q). The parameter qo = 0.2 is chosen 
this large as to clearly show the three different ranges. We will see that the 
error introduced by changing regions I and III by our approximation only 
contribute at order 0(q^) or higher, whereas the main contribution for lower 
orders only comes from region II, if the dashed line is integrated from q$ to 



polynomial in the denominator that cannot be readily solved. We have 
to introduce approximations in order to factorize the polynomial to obtain 
doable integrals, and we shall carefully check that we don't omit any vital 
contributions while applying our approximations. 

3.3.3 Approximations to order T 3 

We want to carefully examine the full integrand of ()3,3,1)) and study which 
kind of approximations can be applied consistently up to order T 3 . For 
T = an exact solution of the real and imaginary parts of the gluon self 
energy Ut and IT^ can be given (see appendix El equations (|D. 1.8)1 . (|D.1.9)1 . 
(|D.1.15p . and HD.1.16)) ). These solutions contain expressions like In \q— qo\ or 
6(2/j, — \qo ±q\) so that the g-integration naturally splits into three regions. 
The full line in figure 13.21 shows the exact integrand. The dashed line in 
the figure shows the integrand that we want to use: We resolve all absolute 
values in the integral in a region qo < q < 2/i, and expand the real part 
Re(Hr + n vac ) through order O(<?o) an d 0(q 2 ). Also, as before, we replace 
9(2/1— \qo±q\) by 9(2/i — q) so that our inverse transverse propagator D^, 1 = 



2/i. 
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q 2 - ql + II T + II vac for T = reads 



ReD 



T 



2 feff i 2 M 
D7T 



+ 



ir 2 q 2 



1 + 



<4 



6vr 2 



2 j^MS 



+ 



80vrV + W V 



-i 

T 



+O(q A ), 
48vrg 3 



(g 2 - g 2 ) (l2,u 2 + 3q 2 + q 2 ) 0(2// - q). 



(3.3.17) 
(3.3.18) 



All powers of higher order in q or qo in (|3.3.18j) are suppressed at least by a 
factor of <? 2 ff . The expansion in the first line of l|3.3.18|) in q is a priori only 
justifiable for anomalous contributions that arise from a region of small qo 
and small q. Actually, the q integral gives contributions all the way up to 
< q < 2fi and higher order contributions in q contribute to the argument 
of the leading logarithm TlnT -1 . These corrections will be suppressed by 
g 2 s at least so that we neglect them for the moment. 

Both real and imaginary part can be extended from region II between 
qo and 2/i — qo to the integration region qo to 2/i, and we should verify that 
this procedure does not introduce an error to the order of interest. In the 
following we will explicitly examine to what order the regions I and III will 
contribute. 

In region I we can argue by the maximum upper bound of the integrand: 
The arc tangent function limits the integrand to tt/2 so that the integral can 
maximally contribute with f Q go q 2 \dq = 7T(/q/6. This is a result of order q\ 

that can be neglected if we expect results of the order gol n <7o, Qb^ 3 ' or %^ ' • 

In region III the argumentation is more subtle. If we expand the exact self 
energy functions in a region 2fi — qo < q < 2fi + q and integrate them in this 
interval, we obtain as a leading order contribution — 727T(7 2 fr 7r// 2 g 2 /(4^(727r 2 + 
13g 2 ff — 12<7 2 ff In 7^)) + 0(<7q). This result is quadratic in the leading order 
of go and could therefore not be neglected as such. However, it turns out 
that this contribution exactly matches the quadratic contribution that one 
would obtain from extending the approximating formula from region II to 
region III and integrating it between 2/i — qo and 2fx. This is also suggested 
by figure l3~2*l where the triangle in region III of the exact solution is replaced 
by a rectangle of half the width of our dashed approximate expression. So by 
integrating our region II formula from qo to 2fi, we actually reduce the error 
to order 0((/q). In the following we can expect to obtain correct expansion 
terms below third order in qo by integrating out our approximating functions 
(13.3.1 7H and 13.3.1 81) from q to 2fi. 
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3.3.4 Momentum integration 

We will integrate (|3.3.1|) using our approximating functions (|3.3.17p and 
(|3.3.18|) . abbreviated conveniently as Re (containing only the explicit terms 
given in (j3,3.17|l omitting higher order contributions) and Im, using integra- 
tion by parts. 



2>M 2 Im q 3 Im 

q arctan — dq = — arctan — 
„ Re 3 Re 



V rtp. g 3 ReIm /_ ImRe ' 
qo ~J qo T Re 2 +Im 2 q - 



(3.3.19) 

The first part of this integral poses no problem, and it is straightforward to 
expand the result in terms of small qo: 

q 3 Im „2 „2 



arctan 
3 Re 



go 67r 



W ®|0(?§). (3.3.20) 



As it should be, this result is independent of whether we start the integration 
from qo or from 0. The second part of the integral is more demanding as we 
get a bulky polynomial in the denominator. (To simplify the following ex- 
pressions, we will set /2ms = 2//, thereby getting rid of logarithmic constants, 
and reintroduce them only in the final result): 

57600//vrV (Re 2 + Im 2 ) 

= q w (57600/i 4 vr 4 + U40g 2 eS fi\ 2 q 2 + 9<? e 2 ff g 4 ) 

+q 8 (75 A 2 [128^+3^-1536^] q 2 + 120g 2 eS » 2 [<4-12tt 2 ] <? 4 ) 
+q 6 (I800 ffc 2 ff [64 + g 2 cS \ fiK 2 q 2 

-20/ [480 5 2 ff 7r 2 - 2880tt 4 + 5cff("92 + 15^ 2 )] g 4 ) 
+q 4 (3600 5 4 ff /x 8 7r 2 g 2 - 600<7 c 2 ffA i 6 [l92vr 2 + 5 2 ff (-16 + 5vr 2 )] q A 

+q 2 (-7200 5e 4 ffAl 8 [-8 + vr 2 ] q* + 600 5 4 ffM 6 7r 2 g« + 100 5e 4 ff /iVc/ 8 ) 
+ (3600 5c 4 ff/ u 8 (7 6 + 600<4mV^ + 25 5c 4 fr /i 4 7r 2 g 10 ) . (3.3.21) 

Bearing in mind that this is just the denominator of the integrand, a direct 
integration seems impossible: This is a polynomial of tenth degree in q, or 
given the fact that only even powers of q appear, we have to cope with a 
quintic at least. Performing rational integrals requires knowledge of the roots 
of the denominator (see appendix IE.1| equation (jE.1.6|0 . But for a quintic, 
solutions in form of root expressions can not be given anymore in general. 
However, all we need is an expansion of the final integral in terms of qo, that 
is for small go- Can we therefore locate those terms of this polynomial that 
are important to low order terms of the qo expansion? A first naive approach 
of simply expanding the denominator in terms of small qo will certainly fail, 
because also our integration variable q can become small at the same time. 
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A possible solution might be found in the pole structure of this expression 
for small go as plotted in figure l3~3l The location of the poles reveals some 
symmetry, and presumably only a few terms from the polynomial (j3.3.21j) 
determine this structure. In the following we will see how to extract these 
terms. 

Pole structure 

We obtain (what we will call) the "low-order front" of the polynomial by first 
retaining only the leading order in go for each order of g separately, and then 
keeping only the leading order terms in q for each order of go- Applying this 
procedure to our polynomial (|3.3.21j) . we first get terms of the orders g 10 , 
Q 8 Qoi Q 6 Qoi Q 4 Qo, Q 2 Qoi an d Qo which are the relevant terms for small go- From 
the three terms proportional to gg we can omit g 8 gg and g 6 gg for small g and 
go so that we are left with only four terms that fully determine the leading 
order behavior of the integrand denominator for small q and go: 

T, 2 , T 2 _ 1 ( 10 . dts^ 4 2 S'cffM 4 ^ 2 - 8) 2 4 fi^// 6 \ 

(3.3.22) 

This "low-order front" uniquely characterize the small go and small g pole 
structure. The first three terms will not change even if we include higher- 
order terms of q or go from the beginning of the calculation in (|3.3,17p : higher 
order contributions in go are shielded by the gg term, higher order terms in 
q are shielded by the g 10 term (for small g where the small go pole structure 
appears) and all other power mixtures in between are shielded by this low- 
order front of four terms. It should be noted however that in the expansion 
of (|3.3.17|) also negative powers of q appear for higher order go corrections 
that add terms of the order (qo/q) 2n with 2n > 6. Their omission turns out 
to be negligible since we restricted our integral to the region go < g and the 
coefficients of this series quickly get smaller. 

Still, we are left with a quintic whose general solution cannot be given in 
form of root expressions. Let us see, if we can reproduce the pole structure 
suggested by figure l3~3l As a first guess we would start with 

Re 2 + Im 2 « -g (A + g 6 ) (F + g 4 ) (3.3.23) 

which would give four poles arranged in a square at a distance F 1//4 and six 
poles arranged in the shape of a honeycomb at a distance A 1 / 6 . Expanded, 
this term gives AF + ^4g 4 + Fq & + g 10 . We can read off A from the g 4 
term A = g 4 ff /i 4 gQ/(167r 2 ) which determines F from the constant term as 
F = gg- This already explains the main features of the pole structures. 
Unfortunately, our low-order front is not correctly reproduced yet. This is 
not dramatic for the additional Fg 6 -term, because it is sub-leading to the 
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Abbildung 3.3: Pole structure of the truncated absolute squared transverse 
propagator Re 2 +Im 2 at \x = //ms/2 = 1, qo = 0.05, g^ s = 1 in the complex q- 
plane. This 10-pole structure arises from the special form of truncation that 
we chose in l|3.3.17|) . For small qo the 10-pole structure naturally decomposes 
into a rectangular structure with the four poles in the middle at the order 
of \q\ (/o, while the six surrounding poles, prominently arranged in the 
shape of a honeycomb, stay at the order \q\ (g^^qo/Air) 1 / 3 . While the 
inner four poles basically contribute to the leading logarithmic contribution, 
it is these outer six poles that give rise to anomalous Fermi-liquid behavior 
of order T 5 / 3 and T 7 / 3 . Corrections to this simple picture will be discussed 
in the text. The plot is multiplied by \q\~ 3 which does not change the pole 
structure since the propagator contains an overall factor of g 6 , but merely 
makes the poles have similar heights in this plot. 
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low-order front, but the term proportional to q 2 is not reproduced at all. 
Therefore we have to extend our guess slightly: 

Re 2 + Im 2 w (A + q 6 ) (F + Gq 2 + q 4 ) = D (3.3.24) 

where we can determine G = —2(ir 2 — 8)qo/ir 2 . In this way we change the 
inner pole structure from a square to a rectangle. We can now calculate 
the integral. Since we know how to factor our denominator, we can apply 
partial fraction decomposition and reduce the integral (j3.3.19j) to doable 
simpler integrals 



where we abbreviated the numerator as N = g 3 (ReIm' — ImRe')/3, the 
original denominator as D = Re 2 + Im 2 , and its first factorization as Do 
from lj3.3.24|) . Applying this procedure, we obtain the following result: 



J qo D 127T V2 g 2 s qo) 9^/3 V 



9ls 
64 



??»t^4'] + <«^ (3-3^6) 



3V3vr 4 /3 V 72 

The result looks insofar good as we obtain the same leading order logarithm 
as from our first approach (|3.3.8j) if we take l|3.3.2f)|) into account, and we 
obtain anomalous q^ 3 and q^ 3 contributions. However, those are not com- 
plete yet, as we will see soon. This integral is only the first of a series that 
we can form: If we denote the terms we neglected in the denominator as 
5Dq = D — Do, we can write an expansion series that should contain the 
complete result: 

AT f N , fN, [N5D 0j f N SD 2 

D dq = J Do-TWo dq = jD- dq -J ~lyf dq + 1 -15T dq ~ + - 

(3.3.27) 

Again, as before, by partial fraction decomposition, all of the integrals on 
the right hand side are doable. Looking at the second term in this series, we 
again find contributions of the order q^ 3 and q^ 3 : 

NSD 0dq = 22 2 /3 (y) ^/3 , 



642V3 (geff ^) 2 /3 g V3 . ^ % 4 ff 
+ nn at aii 1 - - 777^7 + Wo m 9o)- 



D 2 H 9^ 2 /3 V 64 

64 2 1 

"27V37T 4 / 3 V" 256 16384 , 

alculate a correction contribution, 
turns out to be twice as big as the first one. Apparently the series does 



5/3 

Although we expect to calculate a correction contribution, the q term 
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not converge fast enough. The third term in the series J N 5Dq/ D^dq again 

seems promising for it starts at order 0(q^ 3 ) and does not contain a q^ 3 
term anymore, but still we would like to find a more controlled kind of 
approximation. 



Completing the hypercube 

The problem of the series lies in the correction 8Dq which is of the order 
<?o- This order is not enough to compensate negative orders of qo that are 
introduced by powers of the denominator Dq. We could ask if it is possible 
to shuffle the q$ terms from 5Dq into Dq in such a way that Dq remains a 
product of a fourth-order and a sixth-order polynomial. In this way we could 
form a new denominator D\ with a corresponding 5D\ = D — D\ = 0(g 4 ) 
and we could expect a faster converging series in terms of qQ. It turns out that 
it is indeed possible to complete the square - or actually the 10-dimensional 
hypercube - up to a given order of qo so that corrections are higher-order in 
qo by the following procedure: We start from 

Re 2 + Im 2 = 4 {a + bq 2 + cq 4 + dq 6 + eq 8 + fq 10 } 

« (A + Bq 2 + Cq A + g 6 ) (F + Gq 2 + g 4 ) = D x 

= —^M {AF + (AG + BF)q 2 + (A + BG + CF)q 4 

+ (B + CG + F)q & + (C + G)q 8 + q 10 } (3.3.29) 



where a, 6, c, d, e, and / can be read off from (|3.3.21|) . We already determined 
A, F, and G before. 5Dq contains q 2 , contributions for q 10 , q s , and q G , and 
we have three variables left: B, C, and M. This should not be too difficult, 
since the mixing terms BF, BG, CF, ... naturally turn out to be of higher 
order in qQ. We can summarize the procedure as follows: 



M 
.4 

F 
G 
C 
B 



f~l + 0(q 2 ), 



O(q A ), 



a 



AM 
b 

AM 
e-MG 



O(q 2 ), 



M 

^ r O(q 2 ) 



O(q 2 ), 



(3.3.30) 
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where the coefficients are series expanded in qo and truncated to form poly- 
nomials. In our case we get the following assignments 



.4 
B 



4 4 
fleffM „2 

16^°' 

g c 2 ff /i 2 (64 + ff e 2 ff ) 2 
32tt2 90 ' 
3 5e 4 ff + 128 5c 2 ff - 12288 



„ "J C tt 1 "Jell „z 

° 768^ *>' 



Using these in (13.3.2911 we obtain a correction term 5D\ of the order O(q' j. 
It is possible to pursue the same strategy again to get rid of the 0(g 4 ) term 
in the 5D\ correction to form a D2 result by the very same procedure as 
above, including the next order term in each variable. 5D2 would then be a 
correction of order 0(g(j). 

Let us see how this new denominator D\ changes the result. There are 
many terms involved in this calculation, so it is best to let a computer do all 
integrals each of which is doable just as before. Here is the final result: 

- r Idq = - ( i + in + 22/3 ^ )4 f g ° 5/3 fl + 

4 £>i 12vr V2 <4<W 3V3tt 2 /3 V 64 ^ 

2l/ W) 2/3 <7o 7/3 ^ 4: ^ +0 (g 3 In , ). (3.3.32) 



9>/37r 4 /3 V 32 2048 

Another calculation shows that the correction to this result f N 5D\/ D\dq 
is already of order 0(q]^ 3 ^ 3 ). Also, using the second correction f N/D2<iq 
gives the same result up to 0(§q) so that we can now trust the q^ 3 and q^ 3 
coefficients. 

Just to give an impression of how the series would continue, we give the 
(?q coefficient of the expansion (|3.3.32|) 



12tt 3 [ 3 3 320 12288 

-2(vr 2 - 16)\Ar 2 - 4 ( 7T-2arctai] , 2 ) 



3 12 30 ' 3072; ~g 2 cH q 



-8{3tt 2 - lG)!!!^ 2 ^) +0(ql 1/3 ). (3.3.33) 



TO 
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Note that this expansion results from the terms explicitly given in equations 
(|3.3.17|) and l|3.3.18j) . Including higher order contributions of qo and q in 
(|3.3.17|) will change the higher order coefficients of g e ff given here. 



3.3.5 Pressure, entropy and specific heat 

As before, we calculate pressure by integrating over qo and entropy and spe- 
cific heat from derivatives thereof. We introduce the following abbreviations 

cm = -+7fi-^C'(2), (3.3.34) 

2 2 /W) 4 / 3 r(j)C(j) g 2 

C «/3 ^Tfs (1+ 64 ). (3 - 3 - 35) 

8 2V3 (geff ^/3 r(f)C(f) g 2 s g , 
010/3 = 973V37i (1 " 32" " 2048 } - (3 - 3 - 36) 

Note that inclusion of higher order terms in qo and q in equation (|3.3.17p 
may change the higher order coefficients in g e R given here. The leading order 
in g e g for each coefficient is complete though. The results are 

P T ,n b OI^T 2 ( 327T/i \ 

~Ng~ ~ ^2^r^T +Cln J (3 - 3 - 37) 

+c 8/ 3r 8 / 3 + c 10/3 r 10 / 3 + o(r 4 inr), 



S T ,n b gj^ 2 T f 32t^ 1\ 

^ - -^K n ^ + CXn ~V (3 - 3 - 38) 

+% 3 T^ + ^c 10/3 r 7 /3 + 0(T 3 In T), 



Cy,T,n 6 _ ^ 2 r^ n 32^ + gin _3^ 



iV, 367T2 V <4 T ln 2 

+yC 8/ 3T 5 / 3 + ^c 10/3 r 7 / 3 + O(T 3 lnT). 

Figure 13.41 shows a comparison between a full numerical result and the 
series expansion for small T for the transverse ra^-part of the entropy in 
the large Nf limit. Clearly, the leading order contribution of the curves is 
accurately reproduced by the small-T expansion for different g% s of 1, 4, or 
9 in a region 

T/y. < W(2tt 2 ). This is also the region where the complete 
large- iVy result for the low-temperature entropy exceeds the ideal-gas value 
in figure 13.11 
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Abbildung 3.4: Transverse ra;,-contribution to the interaction part of the 
low-temperature entropy density in the large- Nf limit for the three values 
g 2 s = 1,4,9. The heavy dots give the exact numerical results; the full, 
dashed, and dash-dotted lines correspond to our perturbative result up to 
and including the TlnT -1 , T 5 / 3 , and T 7 / 3 contributions. 



3.3.6 Higher order corrections 

We can now apply full temperature-dependence of the self energies and keep 
the renormalization scale /2ms explicitly in our formulae as in section 13.3.21 
Our final result for the pressure reads as follows: Using the following di- 
mensionless quantities 

TT 2 T 2 

X = 1 + -^-, (3.3.40) 



3/i 2 

36/j, 2 67H hms 



Y = 1 + ^--^] I1 J^L (3.3.41) 



and the following abbreviations 



cm = ^ + ^-^C'(2) + ln|, (3.3.42) 
C8/3 " 371^ ( + MY ) {y) ' (3 - 3 - 43) 

82V3 W) 2/3 r( ^ )C( l| ) g 2 s ^ / xN l/3 



5l0 / 3 " g^SV 3 / 3 (1 "32F"2LW^ ) \Yj {3 - 3M) 
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we can write the resulting pressure as 

PT,n b gWT 2 ( 327TU ^ 

In o rr, + C ln 



N g 72.2 ^ 5 2 ffT y 

+c 8/3 T 8 / 3 + c 10/3 T 10 / 3 + <3(T 4 InT). (3.3.45) 

The temperature corrections of X and 1" only contribute to order 0(T 4 ) so 
that in a strict expansion sense, this expression reduces to (|3.3.37p up to 
ln/iMS terms. 

Let us finally note that with the tools presented here, it seems a straight- 
forward but tedious task to calculate the complete 0(T 4 ) contribution to the 
pressure (or the 0(T 3 ) contribution to entropy or specific heat). First, from 
the region II integral we expect a T 3 lnT contribution to the entropy as is 
suggested by the (still incomplete) contribution in equation Ij3.3.33j) . Then 
we have to take a careful look at regions I and III again, which could start 
contributing at this order. Finally, the temperature dependencies of the self- 
energies in equations (j3.3.12j) and (|3.3.13j) will also start contributing at this 
order via higher corrections similar to the X and Y functions above. In a 
full 0(T 3 In T) calculation of the entropy, all of these terms would have to 
be taken into account. 



3.4 Longitudinal contribution 

For the longitudinal "n^'-part from our starting equation H2.2.13J) we cannot 
neglect the real part of the self-energy Hl compared to the vacuum self 
energy as we did in equation l|3.3.2p . We have to add a term 

Ren L ( go , q) = + 0{ql) + 0(q 2 ) (3.4.1) 

which is absent in the transverse part. This contribution does not vanish 
for small go and q and actually simplifies the calculation. This time we can 
expand the arc tangent for small argument arctan x ~ x so that we can write 

lining - q + U L + II vac ) ~ g2 + (^2)^2 ■ ( 3 - 42 ) 

Performing the g-integration as before we get 

- %M( i+ is) ta(i+ S >->}■ (3 - 4 - 4) 

This integral can be performed exactly, so here and in the following step 
there is no need for any truncation so far. From this integral Q we get the 
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pressure as 

PL,n b ^ [°°dq 



N g 8vr 3 J vr e io/ T - 1 



-Q 



■^{(-^H^^-U (3.4,, 



and the entropy density as 



The logarithm can be expanded giving 

= l£? ( ln fS + + + 0(r3) - (3A7) 

3.5 Full result 

Let us take together all the results we obtained so far. We have three main 
contributions to the entropy correction 

•S — S() = S n0Q - nb + ST,n b + Sl,n b (3.5.1) 

with S Q0Q - nb from equation 1)3.2.2)1 . S T ,n b from 1)3.3.38)1 . and S Lj n b from 
1)3 .4. 7)1 . 5o is the ideal-gas value per unit volume obtained from 1)2.5.1)1 
and is given by 

NnJ^T+^tA+N^T*. 



' V 3 45 / 9 15 



Using the following abbreviations 

<*/» = - 8 tS l) (w>4/3 ' (35 ' 2) 

80 2 1 /3r(fx(J| ! ), l2/ , , , 

C " 3 = 27V^ W (3 - 5 " 3) 
we can write the final result to the entropy correction as 

+ c 5/3 T 5 / 3 + c 7/3 T 7 / 3 + 0(T 3 ln T) (3.5.4) 
and for the specific heat as 

Cy-C° v g 2 eSf j?T ( 4 5eff/ i 6 , 

+ |c 5/3 T 5 / 3 + 7 -c 7/3 T 7 / 3 + 0(T 3 lnT) (3.5.5) 
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with the ideal-gas contribution of the specific heat as calculated in appendix 



3.5.1 Extension to full QED and QCD 

Having calculated the large Nf result of entropy and specific heat, the natural 
question is how these results translate back to the full theory of QED or 
QCD. It turns out that for entropy and specific heat of QED and QCD all 
we calculated so far in the first few orders of T <C \i is actually all there is. 

In fact, our NLO expression for the pressure at large Nf (|2.2.13|1 can be 
seen as the starting point for an expansion of small g but finite, smaller Nf 
in the regime T/fj, <C g with an error of order g 4 . This is because bosonic 
loop insertions that we could omit in the large Nf limit are also negligible 
in this regime to the order of interest. The pressure would then be given by 
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(3.5.7) 



where N = 3, N g = 8 for QCD, and both equal to one for QED. The tem- 
perature T is assumed to be the smallest mass scale in the problem. Dt 
and Dl are the spatially transverse and longitudinal gauge boson propaga- 
tors at finite temperature T and (electron or quark) chemical potential n 
obtained by Dyson-resumming one-loop fermion loops, and -D V ac is the cor- 
responding quantity at zero temperature and chemical potential, just as in 
equation H2.2.13I) . This expression is sufficient to obtain results up to (but 
not including) order T 4 in the pressure or to order T 3 in the entropy for the 
regime T/fi -C g. Since the fractional powers T 5 / 3 and T 7 / 3 are included in 
these bounds they indeed give the anomalous behavior also to full QED and 
QCD. 

In order to obtain the anomalous terms through order T 7 / 3 one actually 
only needs the following terms of the real transverse part of the propagator 
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(|3 .3.1 7|) where we only keep the lowest orders in qo/q and qo/n 
ReD- 1 = q 2 (l + 0(g 2 cS )) 

+ (^-i+o( 5c W)+o(^w/^ 2 : 

+ 0(g 2 eS q$. (3.5.8) 

These are the necessary terms to give the same "low order front" as presented 
in equation ()3.3.22|) . 

For T I ii <C g the anomalous contributions we calculated for the pressure 
are actually negligible compared to zero-temperature contributions ~ g ji 
which we omitted in our starting formula ((3.5.7(1 and which can be found for 
QED and QCD to order g 4 in references [ZH] ■ However, when calculating 
the entropy density S = (dP/dT)y y/1 or the specific heat, these terms drop 
out and we obtain the correct low temperature series expansions for QED 
and QCD by the results calculated above in equations ((3.5,4(1 and ((3.5.5(1 . 

A useful way of understanding the organization of the series expansion 
is to make explicit that T/fi is the smallest scale by writing T/fj, ~ g^ s 
with 5 > 0. The terms in the expa,nsion l|3.5.4|l then correspond to the 
orders g^ S ln(c/<7 c fj), g^ ^ 5 / 3 ^ ; and gff^ 7 ^ 5 , respectively, with a truncation 
error of the order g^.g 3S ■ The expansion parameter in this low-temperature 
series is T/(g e Qfj,), which is also the scaleless parameter appearing in the 
argument of the leading logarithm of l(3.5.4p and 1)3.. 5.5(1 . but interestingly 
only after the transverse and the longitudinal contributions have been added 
together. The combination g c f[fi is the scale of the Debye mass at high 
chemical potential, whose leading-order value is mn = g e g/i/n. In fact, the 
calculation of the coefficients in l(3.5.4p required keeping the leading-order 
"hard-dense-loop" (HDL) part of the gauge boson propagator |117( 1118] . in 
particular the dynamic screening in ((3.3.18(1 . but also a HDL correction to 
the real part of the transverse self energy in 1(3.5.8(1 . The above calculation is 
therefore in a certain sense another application of HDL resummation |118j . 
which thus turns out to be necessary also for a perturbative treatment of the 
low-temperature regime T//i <C g. 

Summarizing, we have found that for the entropy and the specific heat 
(but not for the pressure) the expressions given in ((3.5.4(1 and ((3.5.5(1 already 
give the leading order contribution for full QED and QCD in a regime of 
Tjii < g with g 2 eS = g 2 N f for QED and g 2 a = g 2 N f /2 for QCD. 

3.6 Discussion 
3.6.1 Specific heat 

The specific heat Cy = Cy /V is an important quantity for potential phe- 
nomenological applications in astrophysical systems. Figure 13.51 shows the 
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Abbildung 3.5: The perturbative result for the specific heat, normalized 
to the ideal-gas value, to order T 5 / 3 and T 7 / 3 (lower and upper curves, 
respectively) for two particular values of a s in two-flavor QCD (chosen for 
comparability to reference |119| ) and g c g rs 0.303 for QED. The deviation of 
the QED result from the ideal-gas value is enlarged by a factor of 20, and 
the plot terminates where the expansion parameter {ir 2 T) / {g e sp) rs 1. 

ratio of Cy as given by l|3.5.5p to the ideal-gas value C v for QCD with two 
massless quark flavors. We compute the ratio for two values for a s which 
have been used also in reference and which correspond to one-loop 

running couplings with renormalization point 0.5 GeV (full line) and 1 GeV 
(dashed line). The upper limit of the shaded bands give the result to order 
T 5 / 3 and the lower band to order T 7 / 3 . Alternatively we can think of the 
two QCD bands as roughly corresponding to QCD with a quark chemical 
potential of 0.5 GeV and the total variation corresponding to different renor- 
malization schemes with minimal subtraction scale varied between \x and 2fj,. 
The critical temperature for the color superconducting phase transition may 
be anywhere between 6 and 60 MeV [E], so the range T/fi> 0.012 in figure 
13 .51 might correspond to normal quark matter. We will explore this region in 
more detail in the next section, looking at the anomalous correction at the 
critical temperature of color superconductivity. While it is certainly ques- 
tionable to apply perturbative results for a s > 0.65, figure 1331 suggests that 
the anomalous feature of an excess of the specific heat over its ideal-gas value 
may possibly come into play in astrophysical situations, in particular in the 
cooling of (proto-)neutron stars H20| I121j . The anomalous correction 
Cv/C v > 1 should be compared to the ordinary perturbative estimate for 
Cv/Cy based on the well-known (20] exchange term oc g 2 (which requires 
T/n > g). The latter would predict C v /C v < 0.6 for a s > 0.65. 

In figure 13.51 we also show the effect on QED, where g c g = \[^kol rs 
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i/47r/137 rj 0.303. Here the range of temperature, where the specific heat 
exceeds the ideal-gas value, and the deviations from the latter, are much 
smaller (the deviations from the ideal-gas value have been enlarged by a 
factor of 20 in figure l33l to make them more visible). As mentioned, the 
effect of the anomalous contribution to the specific heat remains small in 
QED, but it might play a noticeable role in the thermodynamics of a normal 
quark matter component of neutron or proto-neutron stars. 



3.6.2 Effect near the CSC critical temperature 

It has been argued recently that color superconductivity (CSC) will dominate 
dense quark matter long before non-Fermi liquid behavior becomes effective 
|110| II 1 1 j . The argument is based on simple dimensional analysis of the 
energy scale, which for the 2SC color superconducting gap is given by |11DI 

rr^[T2in[T2^rr^[T2inrr2Tj 

E = bb' fig~ 5 exp{-c/g), (3.6.1) 



6 = 512^(A) 5/2 , & ' = ex P (-i±4 c = ^/J^ (3.6.2) 

and for non-Fermi liquid effects is given by 

E = 6^exp(— xc/g 2 ), (3.6.3) 




24vr 2 iV c 

-3 + 7E--2('(2)), c = — -A (3.6.4) 

7T IVg 



This NFL result is obtained from equating x times the coefficient of the 
leading order T contribution to the TlnT -1 coefficient from the NFL entropy 
correction (|3.5.4|) (for the specific heat 1)3.5.5(1 . b gets another factor e). The 
parameter x gives the kind of NFL correction: for x < 1, E gives the energy 
scale for corrections of the order x; for x ~ 1 we get into non-perturbative 
NFL correction regime. Note that the pre-exponential factor b' in the color 
superconducting energy gap 1)3.6.1(1 stems from non-Fermi liquid behavior 
through the quark-selfenergy |128( 11251 H29( . 

In references |110( llll] . the energy scale for color superconductivity is 
related to the energy scale where NFL effects get nonperturbative. For a 
2SC system, the energy scale of the critical temperature of the color super- 
conductor T c is related to the gap energy 1)3.6.1)1 by T c = Ee yE /tt |127( . It 
is clear that for arbitrarily small couplings g the critical energy T c oc E will 
be larger than the NFL energy E. Thus, nonperturbative NFL effects will 
not play a role because color superconductivity will set in already at higher 
temperatures. 
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Abbildung 3.6: Non-Fermi liquid correction at the critical temperature T c 
of the color superconductor as a function of the coupling g for entropy and 
specific heat. The factor x is the rate of the next-to-leading order contri- 
bution oc TlnT -1 to the leading order (Stefan-Boltzmann) contribution 
oc T. The shaded bands give the result including T 5//3 (lower limit of band) 
and T 7 / 3 (upper limit of band) contributions in the NLO expression. For 
g < 5, non-Fermi liquid corrections are perturbative (x < 1) and the size 
of correction is given by x for the specific heat and the entropy. For g > 5 
nonperturbative NFL contributions cannot be neglected. For asymptotically 
small g — > the correction is linear x = g/(3y2). 
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On the other hand, if we are interested in when perturbative corrections 
cannot be neglected, it boils down to a choice of the size of corrections we 
wish to regard, entering via the parameter x: Indeed, if we want to study 
NFL corrections of the order of g by setting x = g in ()3, 6,3(1 . we find the 
same parametrical energy dependence oc exp(— c/g) for T c and E. Here of 
course the magnitudes of c, c, and the large exponential prefactor b start 
to play a role and one has to be more careful. Solving for the correction in 
T C = E gives 

c . 9 2 , bg 6 
x = g- H hi—; — . (3.6.5) 

From another point of view, this is just the correction x = Cy LO (T c ) / C v ° (T c ) 
evaluated at the critical temperature T c (g) = T c (E(g)) for each coupling g for 
the NLO contribution oc TlnT -1 . Figure l3~o1 shows the dependence of the 
size of the correction x on the coupling g for the specific heat and the entropy. 
According to this plot, below g < 5 we are in the perturbative NFL regime 
where our formulae 1)3.5.4(1 and 1)3.5.5(1 are well applicable, as can be seen 
from the shaded bands in the plots which show the result including T 5 / 3 and 
T 7 / 3 contributions. These are the biggest corrections that we expect to see 
from perturbative non-Fermi liquids because for higher temperatures T > T c 
the NLO corrections get smaller, and for smaller temperatures T < T c we 
are in the color superconducting phase. Since the color superconductor gap 
equation 1)3.6. ip might not be valid in this coupling range anymore, this plot 
should be treated with care. The only purpose of this plot is to show that 
from the gap energy given above, one can not readily derive that perturbative 
NFL effects could be neglected. Indeed, the energy gap formula might break 
down for g as low as g > 0.8 p.30J. 

For small g the correction decreases linearly as x = g/(3y2). This is 
larger than naively expected from the perturbative NLO NFL term oc g 2 In g 
and an effect of the exponential decrease of the critical temperature T c with 
exp(-c/#). 
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Kapitel 4 

Summary and Outlook 



In this thesis we calculated pressure, entropy, and specific heat of the high- 
temperature and high-density quark-gluon plasma in the limit of a large 
number of quark flavors. The theoretical understanding of the quark-gluon 
plasma is of current interest, as existing heavy-ion colliders start to reach en- 
ergy densities that allow for observation of this new state of matter. In the 
low-temperature but high-density region, the quark-gluon plasma features 
non-Fermi-liquid behavior with anomalous contributions to thermodynamic 
quantities. In this thesis we derived for the first time in a detailed calcula- 
tion the low-temperature series beyond leading-log accuracy and including 
anomalous fractional powers. 

We have focused on large Nf as the theoretical limit in which we let 
the number of flavors go to Nf — ► oo and the coupling g — ► such that 
the effective coupling g e g oc g 2 Nf stays of order O(l). We argued that this 
theory is of particular interest because it is exactly solvable up to next-to- 
leading order (NLO) in the 1/Nf expansion and thus provides an ideal test 
for resummation schemes that try to overcome the poor convergence prop- 
erties of strict perturbative expansions of thermodynamic quantities like the 
pressure. The NLO large Nf limit technically corresponds to resumming 
a boson loop with any number of fermion loop insertions which can be re- 
summed by the Schwinger-Dyson method. Although we cannot give a closed 
form of the NLO pressure, we can numerically evaluate the result with a 
small error introduced by a cutoff due to a Landau singularity in the cou- 
pling. Given that the renormalization scale dependence is exact to the order 
of interest, we have an exact result for all couplings and renormalization 
scales. We calculated the pressure for the whole range of temperature and 
chemical potentials where the error introduced by the cutoff is negligible 
numerically. For small values of the coupling, our numerics turned out to 
be sufficiently accurate to allow comparison to perturbation theory and to 
extract coefficients in the perturbative expansion of the pressure that had 
not been calculated analytically before. This table shows a summary of our 
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predictions: 



Coefficient 


Predicted value 


Equation 


Confirmed 


P : C 6 


+20(2) 


(|2.5.1 6| 


not yet 


x- c e 


-4.55(9) 


42.5.19| 


not yet 


d 2 X /dfx 2 : C 4 


-7.02(3) 


(|2.5.2()|) 


-6.9986997.. [M] 


Pt=0- C6,ln 


+3.18(5) 




not yet 


Pt=o- Cq 


-3.4(3) 




not yet 



Since the original publication of the predicted values in [HH] so far only one 
has been confirmed by an analytic calculation. The other four constants still 
need to be confirmed by means of an analytic calculation. 
For larger values of the coupling, we studied the convergence properties of 
perturbative series including renormalization scale dependence and the FAC 
(fastest apparent convergence) scale applied to m E (FAC-m) and to g\ (FAC- 
g) at zero chemical potential. The latter choice leads to quite accurate 
results up to g 2 s ~ 9. We further studied the scaling from pressure at small 
chemical potentials to higher chemical potentials and noticed an unexpected 
breakdown at about \x « ttT. This might indicate a generic obstruction for 
extrapolating lattice data from /i <C T to /i > T. 

In the second part of the thesis we examined the region of small T and 
large \x more closely and found non-Fermi-liquid behavior in an anomalous 
series expansion of the entropy and specific heat for small temperatures in- 
volving a leading TlnT -1 term and fractional powers J , ( 2 ™+ 3 )/ 3 . We gave a 
brief overview of the theoretical picture from the ideal gas to the non-Fermi 
liquid. While the classical ideal gas shows no temperature dependence in 
the specific heat, the ideal Fermi gas and the Landau- Fermi liquid show a 
linear dependence on temperature in the specific heat. Only the introduction 
of long-range interactions changes this behavior qualitatively: The leading 
TlnT -1 contribution in the entropy or the specific heat cannot be explained 
by classical Landau-Fermi liquid theory anymore and the system is therefore 
called a non-Fermi liquid. In our case long-range interactions are introduced 
by transverse gauge boson interactions which are only weakly screened at 
low frequencies go and not screened at all in the static limit of qo — ► 0. We 
showed that the anomalous leading logarithmic contribution leads to a range 
where the entropy exceeds the ideal gas value of the entropy for small tem- 
peratures. We calculated the leading log coefficient in a straightforward way, 
confirming a result by reference |104j . We further completed the argument 
of the leading logarithm and went beyond the leading order to find anoma- 
lous contributions to order T 5 / 3 and T 7 / 3 in the entropy and the specific 
heat. These contributions indeed come from transverse "non-n&" parts (i.e. 
contributions in the original expression of the pressure that do not contain 
a factor of the bosonic distribution function n&) and we calculated them by 
spotting a "low order front" in the denominator of the pressure integral which 
we can factor into a sixth order and a fourth order polynomial, hiding higher 
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order terms in the frequency qo in a correction term. The final result for the 
entropy and the specific heat turns out to be readily applicable to full QED 
and full QCD in a range where T/fi -C g e ft- For the entropy we found the 
result 

+JV s c 5/3 T 5 ' 3 + JV s c 7/3 r 7 / 3 + 0(T 3 InT) (4.1.1) 
from which the specific heat can be derived as 

r at at ^ T ju S^Ng^T ( 4g eSf , 6 \ 

+N g ^c 5/3 T 5 / 3 + N g 7 -c 7/3 T 7 / 3 + 0(T 3 In T) (4.1.2) 

with 

8 2 2 /3r(|K(|) 4/3 

5/3 " 9^3^/3 ( } ' ( } 

- - ^l^(-)-- 

We pointed out a possible implication in astrophysical situations, namely 
the cooling rate of proto-neutron stars, and calculated the effect of the 
non-Fermi-liquid contributions near the critical temperature of color-super- 
conductivity. We found that while non-perturbative non-Fermi-liquid be- 
havior might play a role for g > 5, there are non-negligible perturbative 
non- Fermi-liquid effects for g ~ 0(1) of the order of 10% — 20% compared 
to the leading order contribution, and that perturbative non-Fermi-liquid 
effects are of the order of g/(3y2) for small g. 

In the appendix we calculated the Feynman rules trying to match dif- 
ferent conventions in the literature, gave a derivation of the specific heat at 
next-to-leading order and presented the zero-temperature limits of the boson 
self-energy with one fermion insertion. Finally, we provided a Mathematica 
extension that is necessary in order to correctly calculate the fractional power 
series. This calculation involves a series expansion that depends crucially on 
assumptions about the variables involved that standard symbolic manipula- 
tion programs implicitly ignore. 

Possible extensions to the work presented here first of all include im- 
provement of the non-Fermi-liquid series expansion. Within the large Nf 
limit, the T 3 contribution might be of special interest. With the tools pre- 
sented here it seems straightforward to calculate it, but one has to take into 
account that disproportionately more contributions have to be considered at 
this order. Apart from higher orders in T, the series expansion could also be 
pushed towards higher orders in the coupling g e g. 
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The analytic treatment of non-Fermi-liquid behavior presented in this 
work paves the way for various possible applications: One could try and 
repeat the calculation of neutron star cooling rates with the improved series 
expansion presented here; one could study the interplay of the effects of 
the non-Fermi-liquid contributions with color superconductivity; finally, one 
could extend the analysis to the nonrelativistic case (i.e. include fermion 
masses) which may be of interest in view of the recent experimental findings 
[T7] in condensed matter physics. 



Anhang A 

Abbreviations 



2PI 


2 Particle Irreducible 


2SC 


2-Flavor Color Superconductor 


CFL 


Color-Flavor Locked 


CSC 


Color Superconductivity 


DR 


Dimensional Reduction 


FAC 


Fastest Apparent Convergence 


HDL 


Hard Dense Loop 


HTL 


iicii u xiiciiiicxi J_ j \j ' M 


ITF 


Imaginary Time Formalism 


Large Nf 


Large number of Quark Flavors 


LHC 


Large Hadron Collider 


LO 


Leading Order 


NLO 


Next-to-Leading Order 


PMS 


Principle of Minimal Sensitivity 


QCD 


Quantum Chromodynamics 


QED 


Quantum Electrodynamics 


QGP 


Quark Gluon Plasma 


RHIC 


Relativistic Heavy Ion Collider 


RTF 


Real Time Formalism 


SPS 


Super Proton Synchrotron 



85 



ANHANG A. ABBREVIATIONS 



Anhang B 

QCD Feynman rules 



B.l Motivation for a unifying approach 

The theoretical description of the quark-gluon plasma is based on quantum 
field theory (QFT). Its application to a 'Many-Body Problem' using tools 
from statistical mechanics started in the late 1950s. In the mean-time the 
theoretical tools improved dramatically, for example by reorganized pertur- 
bation theory in the 1990's, and it is no wonder that a plethora of conventions 
and notations have been invented so far. The drawback of this kind of sci- 
entific independence is that it renders comparison between different authors 
often a tedious task. In this chapter we will try to find a formal "common 
denominator" for modern literature on thermal field theory, whilst providing 
a brief introduction to Feynman rules and some basic diagrams we will be 
using later on. 

We will use the Imaginary Time Formalism (ITF) for the statistical de- 
scription of the quantum field. This formalism is based on the observation 
that the statistical density matrix for a system p{(3) = exp(—(3H) with in- 
verse temperature (3 = 1/T looks similar to a quantum mechanical time 
evolution operator exp(iHt) if one uses formally imaginary time t — > i(5. 
It is therefore tempting to simply rotate QFT formulae from zero tempera- 
ture, as taught in modern textbooks like Peskin and Schroeder |131j . in the 
complex time plane to finite temperature in ITF. Unfortunately, naive com- 
parisons quickly result in missing factors of —1 or i, so that it seems worth to 
invest some effort in locating the sources of discrepancies. In the end we will 
see that only two factors determine the sign discrepancies for Feynman rules 
between various conventions, namely the sign in the generating functional, 
which we will call 's, and a factor in the covariant derivative g. 

Let us start by looking at the gluon propagator, which basically resem- 
bles the photon propagator up to group factors. It can be found in various 
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definitions in the literature: 



G%{k) = -9„uA(k) + (A - l)k^k u A\k) (B.1.2) 
= ~3nv~j^- + (A - l)kfj,k u -^, 



^ = ^(V-(w)^)- (B.i.3) 



1 ( A h 

Equation flTTT) for T = can be found in Peskin and Schroeder |131j . 
p. 297, equation l|B.1.2|l applies to the ITF and is the definition of Blaizot 
and Iancu |132j . p. 176, while equation QB.1.3J) is the ITF definition according 
to Le Bellac |133j . p. 107. The first two equations make use of the Minkowski 
metric g^ u = diag(l, —1, —1, —1) while the third equation does not. We can 
unify the first two equations, using the following "generic" equation: 

^) = -x(<T F -(l-0— J (B.1.4) 

where we introduced a "generic constant" % whose value resembles different 
definitions. Here we can identify \ = i for the T = case in equation (|B.1.1|) 
and x = —1 for the ITF in equation (|B.1.2|1 . We cannot easily identify 
equation 1|B.1.3|1 with this method, though, because it uses Euclidean metric. 
We will see later how we can convert results obtained in this work from 
Minkowski metric to Euclidean metric. 

Another example of using different definitions is the covariant derivative. 
In [EH, p. 490, it is defined as = d^-igA ^ , while [EE], p. 171, defines 
the covariant derivative with a positive sign = + igA a ^t a . Again these 
two definitions can be easily unified, using the generic equation: 

D^d^+ggA^ (B.1.5) 

with g = — i in the first case and g = +i in the second. Itzykson and Zuber 
|134j . p. 584, use yet another convention g = — 1. 

Yet another striking difference concerns the definition of the generating 
functional. While |131j . p. 290, defines the generating functional for T = 

as 



Z[J] = J V(j)exp i J d 4 x[C + J(x)(j)(x)} 
we can find the following definition for ITF in |132j . p. 28, 

Z[j] = jV [ D(0)exp(- / dr [ d 3 x(£ E (x)-j(x)0(x)) 



(B.1.6) 



(B.1.7) 
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where the Euclidean Lagrangian Ce basically equals the negative Minkowski 
Lagrangian— C in the common Minkowski metric. Defining a "generic" position- 
integral as 

7 4~ _ / I d4x for T = 



f£ dr f d 3 x for ITF 



(B.1.8) 



we can again combine the two equations to form the generic equation 



Z[j] =Af £></>exp 



J T>(f)exp s J 



d*x (£ + jj(x)4>(x) 



(B.1.9) 



with s = i, j = 1 for T = 0, and s = 1, j = 1 for ITF. 

What we shall find in the following is that the Feynman rules can be de- 
rived basically from these three generic definitions presented here, namely the 
definition of the propagator (|B.1.4|1 , the definition of the covariant derivative 
(IB.l.o|) and the definition of the generating functional ()B.1.9|) . From these 
definitions we can derive the QCD Feynman rules. In the end the Feynman 
rules will in a simple way depend on the constants 9, and ?. 



B.2 Derivation of Feynman rules 

In the following we will present the key equations that define constants like x, 
g, or s. The aim is to unify the Feynman rules for T = and the Imaginary 
Time Formalism (ITF). We use natural units h = c = 1. 



Metric 

Equations in both formalisms can be formulated using the Minkowski metric 



9v» = 9^ = diag(l, -1, -1, -1). (B.2.1) 

In the case of the imaginary-time formalism we adopt the definition of 
Blaizot and Iancu |132j and write x^ = (rEoi x ) = (to — ^ r > x ) with < 
t < P, with the inverse temperature = T _1 , and = (ko,k) = (iu; n ,k) 
with u n = 2nirT for bosonic fields or u> n = {2n + l)irT for fermionic fields. 
This definition has the advantage that the scalar product of position and 
momentum vector still has the form of the Minkowski metric: 

k^x^ = koXQ — kx = u> n r — kx. (B.2. 2) 

We define the generic position integral to be 

,4- forr = 

I tfdTf#2 for ITF. (B ' 2 ' 3) 

Analogous "generic" definitions can be given for the position integral and 
the Dirac ^-function. 
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Lagrangian 

In all textbooks we compared, the Lagrangian takes the following form: 

1 



The covariant derivative If) = is defined as 



(B.2.4) 



(B.2.5) 



with g = i, —1, or —i. We use the commutation relation of the generators 
t a of the Lie algebra of the gauge group in the fundamental representation 



[t a ,t b ] = if abc t c 



(B.2.6) 



where f abc are the structure constants of the group. The commutator of the 
covariant derivative 



[D», D„] = gg (d^A a v - d u A° + iggf abc A b ^ t a = ggF^f (B.2.7) 

defines the field strength tensor F®. The adjoint representation is given by 
(*g)oc = if abc - In this representation the covariant derivative is given by 

(D^ ac = + iggA b ^ abc . 



The generating functional for gluons can be written as 



m = M / [VA»] exp 



s d 4 x (C + jf a (x)A a u (x) 



(B.2.8) 



(B.2.9) 



with s = i for T = and s = +1 for imaginary time formalism. We usually 
find j = 1 in both cases. For the quarks which are fermions we have to 
consider two anticommuting fields and the generating functional takes the 
form 



Z[r), rj\ = Ml VipVip exp 



d x (£ + rjrn(x)ipi(x) + rjijji{x)r}i{x)) 

(B.2.10) 

with rj = 1. For ghost fields there is a similar generating functional. 
Propagators 



Starting from the generating functional (|B.2.9p we can write the two-point 
correlation function as 



< 0|TAS(xi)^(x 2 )|0 > = 



M 
Z[0] 

1 

Z[0] {sj6j^( Xl ))(sj6j vb {x 2 )) 



[VA° e \Afc 1 )JI°{x 2 )e X p 
S 2 Z[j] 



s / d 4 xC 



(B.2.11) 



3=0 
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We can rewrite the generating functional by a shift of the gauge field 
A%(x) = A'^x) + {-~j/x) J d 4 yG^(> - y)f b {y) where the Green's function 
G a ^ v {x — y) satisfies the generic equation 

d 2 g,u - (l - |) d^dj) G vpah (x -y) = x^ ah 5^{x - y) (B.2.12) 

with x = i for T = and x = — 1 for ITF. Here we already added the 
Faddeev-Popov term — (d tJ, A^) 2 /2£ to the Lagrangian. The shift of A"(x) 
does not affect the result of the integral and we can rewrite 1 the generating 
functional (jB.2.9|) as 



Z[j] = Z[0]exp 



2x 



d 4 5? / d 4 yf a G<±(x-y)f b (y) 



(B.2.13) 



Using this result, we can calculate the two-point correlation function (|B.2.11|) 

as 



1 



<0| TA;( Xl )Al(x 2 ) |0> = ^G a ; u ( Xl - x 2 ). 



(B.2.14) 



As expected, this expression is independent of the constant j. We can calcu- 
late the Green's function for the momentum space by Fourier transforming 
equation l|B.2.12|) : 



cab 



(B.2.15) 



This is the form that we already introduced in equation (|B.1.4|1 . We can 
now write the Fourier transform of the two-point correlation function 



(0\TA-( Xl )A b u (x 2 )\0) = e -«(i-») ^G%(k) . 

momentum 

, external fields propagator 
integral 

(B.2.16) 

Written in this way, the connection to the Feynman rules is clear and we 
can readily read off the propagator that will be used in calculating Feynman 
diagrams: 



-1 

SX 



Gt(k) (B.2.17) 



•S 



1 In the literature this procedure is often referred to as "completing the square" 
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It is interesting to see that the last line is actually independent of x and 
thus independent of our choice of defining the propagat or G ab v (x). It merely 
depends on the definition of the generating functional (jB.2.9|) . 

For fermions we can derive the propagator in a similar way starting from 
the generating functional equation (|B, 2,10(1 . Again, we can shift the fermion 
field by introducing the fermionic Green's function Sp(x — y) which satisfies 

(i$> - m) S%{x -y) = a5 {4) {x - y)5 ij l AxA (B.2.18) 

with 114x4 being the unit matrix in the space of the anticommuting 7-matrices 
and a a generic constant. Similarly as above, we can write the two-point 
correlation function of two fermionic fields as 

<0| T^( Xl W{x 2 ) |0> = ^Sj(x - y) (B.2.19) 

SCT 

and we find for the propagator used in the Feynman rules 



~is£(p) (B.2.20) 



P = -} 1 5*--} ^ + m 5->. 

s ji — m s p 2 — m 2 + ie 

Finally we have similar results for the ghost field. The ghost Green's 
function satisfies 

(-dH ah ) G b g c h (x -y) = c5 ac 5^(x - y) (B.2.21) 

with the generic constant c. This Green's function gives the two-point cor- 
relation function 

<0| Tc a ( Xl )^(x 2 ) |0) = ^Gf h (x x - x 2 ) (B.2.22) 

sc a 

and the ghost propagator for the Feynman rules 

a ^ b -1 ~ -1 r) ab 

► =— G? h (k) = — - ¥ . (B.2.23) 



sc 



Vertices 

We use the field strength tensor from equation l|B.2.7|) 



FL = d^Al - d v Al + iggf abc A'Al (B.2.24) 
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to calculate the expression we need for the Lagrangian l|B,2,4|) 

C g i uon = -\F%F^ va = - l -(d^A a u d p A» a - d^A^A^)} (B.2.25) 

-id li A'%gf ab °A'*A ve 

+ ^g 2 g 2 f abc f ade A b 1 A c l ,A pd A ue . 

The first line in equation (jB.2.25|) gave the gluon propagator (IB. 2. 12ft . 
The following two lines can be treated as small perturbations with the order 
parameter g and contribute to the Feynman rules as 3-boson vertex (second 
line) and 4-boson vertex (third line). We can calculate the 3-boson vertex 
for example by calculating the 3-point function 

(n\TA a ^x)Ai(y)A c p (z)\n) = (0\AAAexp(jJ d*x(dAggfAA)\0) (B.2.26) 

= (0| AAA |0) + sggf (0| AAA J(8A)AA |0) + ... 

and applying the Wick theorem. After some careful calculation we get for 
the 3-boson vertex 

a,/i 

(Sg)gf abc {g pv {k-p) p (B.2.27) 
+g vp (p-qT 



where we use the short notation (k — p) p = k p — p p . For the 4-boson vertex 
we can start similarly from a 4-point correlation function and obtain 

&)g 2 (f ahe f cde (g pp g va -g^g" p ) (B.2.28) 

+f ace f bde {g ^ gP a_ g ^ g u P) 
+f ade f bce {g ^ gP a_ gWg ua ) \ _ 



The part of the Lagrangian (jB.2.4|l that describes the fermions is given 

by 

Cquark = A (ip - m)^ A = A (ip - m) S^j + igg^a^^fijA (B.2.29) 
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where we used the covariant derivative il/) = i^d^ + igg^^A^t 0, from equa- 
tion IB. 2. Fiji . The first part of equation ()B.2.29j) gave the fermion propagator 
(|B.2.20|) while the second part gives the quark-gluon vertex 



= (i£g)gj%- (B.2.30) 



* 3 




Similarly we can write down the Lagrangian part for the ghost fields 



Cghort = c° {-WD") c c = cf (-d^8 ac ) c c + cf (-iggf^A^j c c 

( B.2.31 ) 

where we used the adjoint representation of the covariant derivative 1|B,2.8|) , 
Note that the derivative in the second part acts on both the gauge field 
A b ^ and the ghost field c c . The first part determined the ghost propagator 
()B.2.23j) while the second part gives the following Feynman rule: 



{-89)9? f 



fi rabc 



(B.2.32) 



In this Feynman rule denotes the momentum of the outgoing ghost . 
Euclidean metric 

We can now transform our results into Euclidean metric. In deriving the 
Feynman rules we did not make use of the metric g^ v other than contracting 
momentum vectors like k^k^ or momentum and position vectors as in the 
Fourier transform k^x^. The latter is taken care of by our generic definitions 
of the integrals l|B,2.3|) and delta functions. The former just needs "undo- 
ing" the contractions. Since we started with lower indices at the beginning 
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(jB.2.26|) . resulting in propagators with lower indices and vertices with upper 
indices, we now have to "pull down" all indices. That is 

k » ^ gr'kp,, (B.2.33) 
k = k^ L k Pi > g^ k p 'k p 

such that in our formula there are only vectors k p , p p , q p , ... with lower 
indices. 

We could have done the complete calculation by using the Euclidean 
metric this way, so we can simply replace 

9»v = 9^ - -<V (B.2.34) 

and write k 2 = g^ v k p k v — > —b pv K p K v = —K 2 , where we use capital letters 
to denote vectors K 2 = K p K p = K 2 + K 2 + K 2 + K 2 in the Euclidean 
metric. 

For example the gluon propagator l|B.2.15|l becomes 
™b n„\ ~ ( „ , it 1^ kfjiv \ a ij 



-x (-V^ + « - 1)^) S" 



For vertices we have to pull down all indices, so for example the 3-boson 
vertex [IB. 2. 2711 chan ges to 

(Sg)gf abc [g^g pp '(k - P ) p , + g^g^'( P - ?V + g^g vv> {q - fc)„) 
- (£g)gf abc (+^AK-P) p + 5„ p (P-Q) p + 5 Pfl (Q-K) u ). (B.2.36) 

In the case of slashed quantities p 1 we apply the same transformation and 
write 

i> = 7% = ^'l&W - SwIlyPw = -V- (B.2.37) 
The fermion Green's function l)B.2.20p then becomes 

Sf 1 \ P ) = a—^—5 ab -5— — 5° 6 . (B.2.38) 

v ; i> - m P + m K ' 



B.3 Comparison with literature 

Values for constants 

We want to compare our "generic" results to various results of Feynman 
rules in the literature. The results are summarized in table IB. H The first 
constant s describes the sign in the Lagrangian of the generating functional 
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Reference 


M/E 


T 


s 


9 


X 


a 


c 


Peskin & Schroeder 11311 


M 


To 


i 


—i 


i 


i 




Itzykson & Zuber 11341 


M 


T 


i 


-1 


i 






Le Bellac, QSFT |135J 


M 


T 


i 


+i 


i 


i 




Huang [136J 


M 


To 


i 


+i 


i 






Le Bellac, TFT-M [133J 


M 


RTF 


i 


+i 


i 


i 




Blaizot &z Iancu j!32j 


M 


ITF 


+1 


+i 


-1 


-1 


-1 


Kapusta [T37J 


M 


ITF 


+ 1 


+i 


-1 


+ 1 


+1 


Le Bellac, TFT-E [133J 


E 


ITF 


+ 1 


+i 


-1 


-1 





Tabelle B.l: Comparison of definitions given in literature. "M/E" means 
the metric (Minkowski or Euclidean), T stands for temperature T = (To), 
Imaginary Time Formalism (ITF) or Real Time Formalism (RTF). Le Bellac 
|133| gives two sets of Feynman rules which are compared independently. 
Only s and g are relevant for the sign changes in the Feynman rules. 



Z = exp(s J C) and can be obtained by comparing equations in the literature 
to equation l|B.1.9|l . Another way is to compare the propagators of the 
Feynman rules as in equations l|B.2.17|l . I|B.2.20|I . and (|B.2.23|) . These are 
independent of the sign of the Green's function (x, 5, c) and are usually 
given in a table of Feynman rules. A first consistency check can be applied 
here to see whether the definitions of the three kinds of propagators (gluon, 
quark, and ghost) are consistent with each other. 

The next constant g appears in the definition of the covariant derivative 
Dp = dfi +99^t a l)B.2.5p . The next three constants x, 5, c for the gluon, 
quark, and ghost Green's functions respectively are not used in all books. 
They are not essential for the Feynman rules, though, as they do not enter 
the equations of the Feynman vertices and Feynman propagators. 

Vertices 

Table IB. 21 shows the prefactors of the Feynman rule vertices of the 3-gluon 
(IB.2.271) . 4-gluon l)B.2.28ll . quark-gluon 1B.2.30jl . and ghost-gluon ljB.2.3211 
interactions. Before results from literature can be compared to our generic 
results, some obvious transformations have to be made: These include rela- 
belling indices, using the antisymmetry property of the structure constants 
f abc , or flipping signs because of different orientation of momentum vectors. 

Our results agree with most of the authors. At T = in the Minkowski 
metric our results perfectly agree with Peskin and Schroeder |131| . Itzykson 
8z Zuber |134j and Le Bellac |135j . In the imaginary time formalism in 
Euclidean metric our results agree with Le Bellac |133| . We calculated the 
prefactors for Euclidean metric in table IB .21 by applying the transformation 
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Reference 


3-Gluon 


4-Gluon 


Quark 


Ghost 


Generic - Minkowski 


sg 




isg 


-s? 


Peskin & Schroeder 11311 


1 


—i 




-1(*) 


Itzykson k Zuber [134J 


—i 


i 


i 


-N 


Le Bellac, QSFT [T35J 


-1 


—i 


— i 


1 


Huang [136J 


-\£-i< 


-i^i < 






Le Bellac, TFT-M [133J 


-1 


—i 


— i 


1 


Blaizot & Iancu j!32l 


iO 


-l(t) 




-i(T) 


Kapusta j!37j 




-1 


-1^+1 < 


— i 




Generic - Euclidean 


sg 


sg 1 


-isg 




Le Bellac, TFT-E [133J 


i 


-l 


1 


-N 



Tabelle B.2: Prefactors of Feynman rule vertices. The vertex equations 
are given in (IB.2.271) . HB.2.28H . (lB.2.3f)^ . and HB.2.32II for the 3-gluon, 4- 
gluon, quark-gluon, and ghost-gluon vertices respectively. If results differ 
from literature (marked by <) then the l.h.s. of ^ denotes the generic value 
and the r.h.s. denotes the value from literature. The ghost vertex of Peskin 
and Schroeder (*) is wrong in the book (+1), but corrected in the online list 
of errata |138j . Blaizot and Iancu (') do not give the Feynman rules explicitly, 
so the values in the table are the generic ones. Peskin and Schroeder and 
Huang (*) do not give the quark indices in the graphs, resulting in a possible 
±1 factor from the antisymmetric structure constants. 
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rules from Minkowski to Euclidean metric of section IB .21 

The results of Itzykson & Zuber |134j use structure constants C a bc = if abc 
and momentum vectors pointing outwards instead of inwards. Taking care 
of these obvious changes, their results are in accordance with our generic 
Feynman rules. 

Our results do not agree entirely with Kapusta |137j in imaginary time 
formalism with Minkowski metric. There are relative sign errors —1 in the 3- 
gluon vertex and in the quark-gluon vertex. His group generators are labelled 

(~ia j.a 

In the case of Huang |136j Feynman rules disagree in a weird way for the 
3-gluon vertex, the 4-gluon vertex and the ghost vertex. In the case of the 
ghost vertex, the momentum of the boson enters the Feynman rule instead 
of the momentum of the anti-ghost as in the rest of the literature cited. 

Peskin and Schroeder |131j do not provide quark indices i, j in their Feyn- 
man graphs. A swapping of the quark indices results in a relative minus sign 
— 1 stemming from the antisymmetric nature of the representation matrices 
(t a )ij = if ta ^ = —{t a )ji. It is assumed that incoming quarks correspond to 
the second matrix index and outgoing quarks correspond to the first matrix 
index, matching the indices in the fermion interaction Lagrangian (jB.2.29|) . 

Conclusion 

Only two generic constants determine the prefactors of all Feynman rules: ? 
and g. The factor J is determined by the generating functional l|B.1.9|) . It 
is i in the case of T = or the Real Time Formalism (RTF). In the case 
of Imaginary Time Formalism (ITF) it takes the value +1, corresponding to 
the change of sign in the time-component from xq — > —it and dxo — ► — idro 
in the exponent of the generating functional. 

The second essential generic constant is g. Its choice is less motivated 
physically, but merely a matter of taste. 

Together, these two constants (s = i or +1, g = —i, —1, or +i) explain 
the prefactors of propagators and vertices used in Feynman rules. 



B.4 Summary of rules 

The constants s and g can be determined by comparison with the following 
formulae. In standard literature g assumes the values —i |131j . —1 |134j . or 
+i |139| 1133"} I135| . The constant s takes the value i for zero temperature 
field theory or real time formalism and +1 for imaginary time formalism. 
An overview of these constants can be found in the tables IB. II and IB .21 of 
section IB. 31 
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Lagrangian 

= -i (s^sm" - a^aM^) - ^(a^) 2 
-id li AZggf abc A' tb A vc + /^^MjXA^ 6 

+^ - m) 5yVj + igg^a^Affijiij 

+c a (-d^d^5 ac ) c c - iggf abc ^d ti (A b lJL c c ) (B.4.1) 

Generic integral and Generating functional 

A~-fJ d4x forT = 



d x = < „/3 , „ „ „ TrT ^ < — ► Z = D..exp s / d x£ 
\ J* dr / d 3 f for ITF J \ J 

(B.4.2) 

Covariant derivative and Field strength tensor 

L> M = a M + 55^ a F^ = a^-^ + %/ abc 4^ (B.4.3) 

Feynman propagators 

In the following, the arrow means conversion from Minkowski space result 
to Euclidean space result (Minkowski — ► Euclidean). 

a, a b,v —l-i, 

= ~~Gf y {k) (B.4.4) 



1 / 1 , //. -, \ k^k u 
s 



^/(p) (B.4.5) 
so 



p = -} 1 sa -} ^ + m sa 

s p 1 — m s p 2 — m 2 + ie 

s ^ + m s P 2 + m 2 



—Gf h {k) = - J- 2 (B.4.6) 
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Vertices 

(Minkowski — ► Euclidean) 
a,fi 



c,p 




+g vp {p-qY 

+9 m {q-kY) 



(lfg)gf abc {5^{K-P) p 
+5w(Q-K) u ] 



(B.4.7) 




(3?V (f^f^ig^g^-g^g^) (B.4. 

+f ace fb de {g ^ gP a_ g ^ fP) 

+ f ade f bce {g pv g pu — g w g va 

+f ace f bd *(S^8 P *-5^6 up ) 
+f ade f bce (S^ p(T -S w S l/a ) 




(™g)g7% 



(B.4.9) 



a c 



(-S?)5^/ o6c 



(B.4.10) 
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Green's functions 

- (l - i) 0A) - V) = xW ab 5^(x - y) (B.4.11) 

(B.4.12) 

(i# - m) - y) = a5^(x - y)S ij Ux4 (B.4.13) 

Sg(p) = a-r^—fi = a „ ^ + 9 m , ^ - -5?_i— ^" (B.4.14) 
p — m p A — m £ + ie if + m 

Gj^(x - y) = H ac 5 (4) (x - y) (B.4.15) 
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Anhang C 

Thermodynamic relations 



C.l Specific heat at NLO 

The specific heat can be calculated from the entropy by a transformation of 
thermodynamic quantities which can be found in [116J 



C v = C V /V = T{[ — \ -y^M (C.l.l) 




where S is the entropy density 

S = S/V = 

and M is the particle number density 



dP 
df 



MsN/v = {^)„ (ai - 2) 

In order to calculate the next-to-leading order (NLO) specific heat in large- 
Nf we have to expand this equation around the leading order contribution 
(|2.5.1j) . The denominator takes the form 

and the specific heat can then be written as 

L-y = t-v + Ly 
with the leading order (LO) contribution 

, (c ,,) 



v + 
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and the next-to-leading order (NLO) contribution 



/-iNLO 



T < 



* 3 



d 2 P 
dfidT 



T2 , At 
T + ^ 



+ 



d 2 P 



3 



i A 4 



+o(%: 



(C.1.5) 

with NLO-contributions of 5 and P. The ideal-gas limit of the specific heat 
contains, apart from the LO-contribution, a term proportional to N® and is 
given by 



v 



NN f 



7tt 2 T 3 n 2 T 

h - — 

15 3 



4/x 2 T 3 



T 2 



+ N a 



4ir 2 T 3 



72 



15 



(C.1.6) 



C.2 Derivative relations 



For the calculation of quark susceptibilities from the pressure in section l2.5.3| 
we evaluated the derivative with respect to a squared quantity. The deriva- 
tion with respect to a squared quantity can be obtained as follows, using 



fi 2 = t 



df^) df(Vt) 



<V 8t 
The chain rule then works as follows 



f'iVt) 



1 



2Vt 



Similarly for the second derivative we obtain 

,2 ^r/..n2 



</(/(m)) 



df(li) 
d[i 2 



2 



(dfX 



2\2 



g(M)=g"(f(fi)) 



g/M Y 



d 2 f(») 



J +9 ' (fifx)) \d^) 2 - 



(C.2.1) 



(C.2.2) 



(C.2.3) 



In the case of a quadratic lowest order term of the series expansion of /(ju) 
we can further write 



_g! 

dji- 



7(m) 



A*=0 



(C.2.4) 



This is however only valid if /(/x) = 0(fi 2 ). If /(//) contains a linear term, the 
l.h.s. vanishes, while the r.h.s. diverges. Similarly, for the fourth derivative 
we see that dV/<V = 24 but d 2 (/i 2 ) 2 / (d fi 2 ) 2 = d 2 t 2 /dt 2 = 2 so that we 
can write the following relationship if /(/i) only contains positive even orders 
in \i up to 0(fi 4 ), i.e. /(//) = c/i 2 + 0(/i 4 ): 



:/(**) 



12 



2 



(1=0 



(<V) 2 



/(/*) 



(C.2.5) 



Anhang D 

Gauge boson self-energy 



D.l Exact result at T = 



In this section we will extend the calculation of the leading order gauge boson 
self-energy from section 12. Ml to the limit of finite chemical potential \x and 
zero temperature T = 0. In this region, the fermionic distribution function 
from equation (|2,2.14|1 is given by the step-function 

n f (k, T = 0, /i) = i {B{ji -k) + 0(-n - k)) (D.l.l) 

and we can explicitly calculate the real and imaginary parts of the bosonic 
self energy for Minkowski space. 



Imaginary part 

In the equations (|2.3.25|l we can calculate the functions Fi using the T — > 
limit of the fermionic distribution function (jD.l.ljl which give 

Fi(x) = [ n f (k)dk = OZ!o(jji-x), (D.1.2) 

F 2 (x) = / kn f (k)dk = — e(fjL-x), (D.1.3) 

J X ^ 

F 3 (x) = / k 2 n f {k)dk = ^— 0(ji-x). (D.1.4) 

J X " 



The resulting functions can be further simplified. Using the following defi- 
nitions 

U(q ,q) = (2n-\q + q\)e(2fj,-\q + q\), 
V(qo,q) = (2fi-q-q )(2fi + 2q-q ), 

W(q ,q) = 6(qo + q)V(qo,q)+6(-q -q)V(-q ,-q), (D.1.5) 
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we get 

lmG(q ,q) = ^(<? 2 - 9o) MM,</) - U(-\q \,q)} , (D.1.6) 
lmH(q ,q) = {W(\q \,q)U (\q \,q)-W(-\q \, q)U(-\q \,q)} .(D. 1.7) 

From these we get the self-energies 

-9 2 eS (q 2 - ql) 

48vrg 



Imn L (g ,g) = eff ; o _„3 ° J {W(\q \,q)U(\q \,q) 



-W(-\qo\,q)U(-\q \,q)}, (D.1.8) 
lmU T (q ,q) = {[ W (\q \,q) + 6q 2 ]U(\q \,q) (D.1.9) 



- [^(-ko|,g)+6g 2 ] C/(-koU)}- 



Real part 



Also for the real part we can give the exact solution in the limit T — > 0. 
Starting from (|2,3,21|) and 1|2.H.22|) the momentum integration can be ana- 
lytically done since the distribution function reduces again to a step function 
for zero temperature. Rearranging the result, we can write 

R(qo,q) = (qo + q)ln\q + q\-(2n + q Q + q)ln\2n + q + q\,(D.1.10) 
S(qo,q) = {-12n 2 q + {2q-q )(q + q ) 2 )ln\q + q\ (D.l.ll) 
+ (2fi -2q + q ) (2^ + q + q ) 2 In \2fi + q + q \ . 

We use the following abbreviation to keep a lot of terms in dense notation 

Rt(qo, q) = R(q , q) - R(q , -q) + R(-q , q) - R(-q , -q) (D.1.12) 

and similarly for S± (just replace R by S). Then the real part of G and H 
can be written as 

ReG(q ,q) = ^ L 2 + ^^i^(« , ?)) , (D.1.13) 
ReH(q ,q) = + ^S£(«o, q)j . (D.1.14) 

Therefore the exact solution for the real part of the self-energies in Minkowski 
space is given by 

Ren L (go,g) = fa + ^±(go,g)) , (D.1.15) 

Ren T (g ,g) = {l6fi 2 q(q 2 + 2q 2 ) (D.1.16) 

+(q 2 - ql) {6q 2 R±(qo, g) - s£( qo , q)) } . 
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D.2 Approximations to second order in frequency 

We can calculate the following T — > limits in Minkowski space by using 
our exact formulae for the imaginary part of the self-energies l|D.1.8|) and 
(ID.1.9D . In the following we will use these equations in a regime where 
< qo <C q, fJ,. This is a reasonable assumption, as the bosonic distribution 
function n& will probe smaller and smaller qo if we let T — ► 0. In this limit we 
can simplify our functions l|D.l.S|l to U(qo, q) — > (2/i — qo — q)9(2/j, — q) and 
W(qo,q) — > V(q , o,g). (Actually in U(qo,q) we could have used 6(2fi — qo — q) 
which would make the following expressions slightly more complicated, but 
it naturally reduces to 9(2[i — q) as soon as we expand qo around 0.) The 
functions G and H from l|D.1.6|) and l|D.1.7|l then become 

lmG(q , q) = g ° (g ° ~ ^ 0(2/x - g), (D.2.1) 
o7rg 

Imtf(g ,g) = g ° (lV +;°" 3g2) ^- g ), (D.2.2) 

where the ^-functions come from the step-function shape of the fermionic 
distribution function at zero temperature. These are exact solutions in the 
limit of T — ► and qo — > 0. Furthermore, for the real parts of G and H we 
can provide the following expansions in small qo: 

ReG(q ,q) = ^ { 4/x 2 + q [(2fi - q) In |2/x - g| + 2glng 

-(2/i + g) In \2fi + g|]} + O(q 2 ), (D.2.3) 

Retf(g ,g) = ^2^{ 8 ^ 2 9- (2^-<?) 2 (/i + g)ln|2 / u-(/| + 2g 3 lng 

+ (fi - g)(2/i + qf In |2// + g|} + O(g 2 ), (D.2.4) 

with the chemical potential fx. 

We can combine the functions G and H to form the self-energies 11^ and 
Ht according to l)D.1.8|l and l|D.1.9|l . Again, for the imaginary part we get 
nice analytic expressions in the limit of small T and qo'- 

t tt I \ 2 Qo(-Qo + g 2 )(-lV + 3g 2 - gg) 

Imn L (g , ?) = -5efr y — - q), (D.2.5) 



) (-g 2 + g 2 )(12^ + 3g 2 + g 2 ) ; 
48vrg 3 



Imn T ( go ,?) = -ff c ff ~ % U/ A(^ " ?)■ (D.2.6) 



Since this is valid for small qo, it makes sense to expand these expressions 
for small qo. The first coefficient in the go-expansion reads 

Imn L (g ,g) = -g 2 cS M/ f + q2) q 9(2f, - q) + O(g 2 ), (D.2.7) 

o7rg 
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lmU T (q , q) = -ff c V ^ g j q Q 0{2» - q) + 0(q 2 ), (D.2.8) 
and for the real parts 

ReU L (q ,q) = { 8f i 2 q - (2fi - qf{n + g) In |2/x - q\ + 2g 3 lng (D.2.9) 

- g)(2/x + q) 2 In |2/i + g| } + O(gjj), 

Ren T (g ,g) = ^{V g + (2|i-g)(V+W + 2g 2 ) In |2 M - g| (D.2.10) 
+4g 3 lug- (2^ + g)(2/i 2 - /xg + 2g 2 ) In |2,u + g| } + O(gjj). 

For the temperature-independent vacuum contribution the imaginary part 
vanishes for small go, Imn vac (go, g) = 0{q^) while the real part takes the 
value 

ReII vac (g , q) = -<&^ ( ln ^ ~ |) + ( D - 2 -H) 
with the renormalization scale /2. 



Anhang E 

Fractional power expansion 



E.l Introduction to Puiseux series 

The first time fractional powers appear in (|3.3.8p . In this section we want to 
analyze where fractional powers come from and which terms determine their 
coefficients. 

Series involving fractional powers are also known as Puiseux series |140j . 
They appear because one expands around a singular point that does not 
have a Taylor or Laurent series with integer powers at this point. A simple 
expansion where the singular point is inherent is 

3/2 

exp(v^) « 1 + Jx~ + J + ^- + 0{x 2 ) (E.l.l) 
2 o 

which can be easily derived by variable substitution x = y 2 . Another kind 
of example is 

, T 4/3 r 7/3 

f{x) = {x + x 2 f 3 « x 1 / 3 + V - V + O(* 10/3 )- (E.1.2) 
In this case the famous Taylor series formula 

f(x) « /(0) + f(0)x + /"(0)y + ... (E.1.3) 

would give + oox — oox 2 /2 + ... and clearly fails to work. One has to pull 
out the singular part of the function and perform the Taylor series on the 
rest: 

( x + X 2)V3 = x l/3 (1 + x) l/3 = x l/3 + + 0(X 3 )) (E.1.4) 

which gives the result above. In this way, one can obtain fractional power 
series also from more complicated function compositions. 
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The simplest example of fractional powers of the kind we encountered 
appear in the following expression (for positive A, M, B, and C) 

?m " Aq 5 q d = Aq ^ Mgg^ 



M<f + Bq* + Cq*<fi 6M Bq 2 Q 

ACtT 5/3 AC 2 1T 7/3 

9^/3S 2 /3M 4 /3 9 ° 27V3S 4 /3M 5 /3 9 ° 

^(6g 3 + 9^ 2 CgLx-^ 3 gLx) 3 
36B 2 M 2 % 

5ACWq 1/3 7AC 5 irql 3/3 



3 6 V3£ 8 / 3 M 7 / 3 3 7 V3-B 10 / 3 M 8 / 3 
+ O(g 5 ) (E.1.5) 

Without the q 2 q^ term in the denominator (that is C = 0) we would just get 
the leading logarithmic contribution and integer powers of qo (as one can also 
see from the series by setting C = 0). It is this additional term Cq 2 q 2 that 
determines additional fractional power terms. The integral can be calculated 
using the following formula |141j for square-free g(x) and deg(/) < deg(g): 

l^lda = Y^r\n(x-a) (E.1.6) 

where a is the set of all roots of g(a) = 0. In our case we get a sum over 
roots of a cubic equation which first have to be expanded for small qo. We 
can motivate the exponents by 

(\ X / x*^ ^ x*^ x^ 
x + x 1 ' 3 ) = -lnx + x 5/3 1 hO(i 4 ) (E.1.7) 
/ 3 2 3 4 

but we cannot explain factors of ir/y/3 appearing in the coefficients of l|E.1.5|) : 
They only appear together with pure fractional coefficients like T 5//3 and 
T 7 / 3 , but not in the T 3 coefficient. We can give an idea of their appearance 
by this simple sum over roots 

(a 2 +a)\na = -= (E.l. 8) 

a£{x\x s +l=0) * 

but to fully calculate the complete coefficients of the series in l|E.1.5p it seems 
unavoidable to plug in the roots of the denominator into (|E.1.6|) . 

In the course of the qo expansion, terms from the roots of the denomi- 
nator of HE. 1.511 contain expressions like -9BVM + ^3 V / 27B 2 M + 4C 3 ^ 



which upon symbolic series expansion reduce to -9B^M + 9\/Wm. This 
expression is zero for positive B and non-zero for negative B which would 
result in two completely different series. It turns out that from the beginning 
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of the expansion, B and M have to be assumed to have a certain sign, in 
our case to be positive quantities, in order to get the correct series. Unfortu- 
nately, symbolic manipulation programs do not readily allow for this choice 
(in the course of an expansion, B is not considered equal to which 
in our case means the implicit assumption that B is negative.) A possible 
way to accomplish the series expansion also for positive B in Mathematica 
is presented in the next section. 

For the sake of completeness we also give the other integrals with the 
same denominator as (|E.1.5|) but different odd powers of q in the numerator: 

rgmax Aqq 3 ^ An 5/3 ACtt 7 /3 

J qo Mq^ + Bq 2 + Cq 2 q 2 9 ~ 3V3B 2 / 3 MV 3% 9^B^M^ % 

A(3B s + 6B 2 Mqi ax -C 2 qt iax ) 3 

i2Bm q ^ % 

bACW 1/3 7AC^qf /z 



3 5 V3B 8 / 3 M 4 /3 3 6 y35 10 / 3 M 5 /3 
+0{ql) (E.1.9) 



dq 



7/3 



30 



Mq G + Bq 2 + Cq 2 q 2 H Sy/^B^M 2 /^ 

A(3B + C 2 q 2 max ) 3 
— -% 



QBMq 2 

AC 2 7rql 1/3 4ACW 3/3 



3 3 V3J3 5 / 3 M 4 /3 3 B V3J3V3M 5 / 3 
+O(q 5 ) (E.1.10) 

Note that equation l|E.1.9jl contains the same fractional power coefficients as 
equation 1|E.1.5|I up to a factor — C/(3M). This is because the combination 

is a series in integer powers of go- 



E.2 Workaround for Mathematica 

Assume positive 

The calculation of fractional powers crucially depends on a series expansion 
which in the current version of Mathematica |142|I143] is not straightforward 
to calculate. As a simple example, let us calculate the series of 



-A + VA 2 + x 

(E.2.1) 



-A + VA 2 + 2x 
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Expressions of the form of the denominator or numerator naturally appear 
when solving quadratic, cubic, or quartic equations. The series of s delicately 
depends on whether A is positive or negative. We get the following series 
expansion for \x\ < A 2 /2: 



1 x x 2 19x 3 25x 4 t. „ 

2 + ^-^ + 12^-lW + ° (s) forA> °< 

l-^ + ^-^ + ^ + ^x 5 ) forA<0. 
4A 2 16A A 16A 6 256A S v ; 



(E.2.2) 



Trying to calculate the series symbolically results in subexpressions of the 
form —A + VA 2 which is not automatically reduced to 0. Therefore Math- 
ematica does not recognize that for positive A l'Hospital's rule has to be 
applied to calculate the series elements as x — > 0. Unfortunately, even 
in the most recent version of Mathematica (which by the time of writing 
is Mathematica 5.0 |143j ) commands like Assuming[{A>0}, ...] or options 
like Limit[s ,x— >0, Assumptions— >{A>0}] are applied only before or after 
a series calculation and not in the course of the calculation, and there- 
fore only provides a solution which is valid for —A + yA 2 / which is 
satisfied by negative A. Trying to plug in positive A in the resulting se- 
ries unfortunately gives the wrong result, as can be checked numerically or 
by plotting the functions. We get a complete nonsense result by the sim- 
ple Mathematica command s+0[x]~3//PowerExpand — ► 1+ComplexInfin- 
ity x+Indeterminate x"2+0[x]"3 since series expansion (implicitly) assumes 
A < while Power Expand assumes A > 0. 

In order to get the correct result, we have to modify the standard be- 
havior of the internal Power[..] function: Whenever we find an expression 
(x a ) b with positive x, we immediately expand it to x ab . Our rules also in- 
clude (x a y) b — > x ab y b for x > and its standard-case for a = 1. We use 
Length[Cases[ListOfPositives,x]] > to check whether x is part of varlist of 
positive variables. The attribute HoldRest is set first, so that expressions are 
only evaluated after the power rules have been modified. After calculation, 
the power rules are reset to their original state, and the new result is re- 
turned. The following routine also expands the Log[..] function accordingly. 

SetAttributes [AssumePositive , HoldRest] ; 
AssumePositive [varlist_ , expression_] := 

Block [{ListOf Positives = varlist}, 

(*Andreas Ipp, Sept 12, 2003*) 

(*First we change Power Rules*) 

Unprotect [Power , Log]; 

Power [Power [x_ /; (Length [Cases [ListOf Positives , x]] 

> 0), n_.]*y_, m_] := Power [x, n*m] *Power [y, m] ; 
Power [Power [x_ /; (Length [Cases [ListOf Positives , x]] 
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> 0), n_] , m_] := Power[x, n*m] ; 
Log [Power [x_ /; (Length [Cases [ListOf Positives , x]] > 0) , 

n_.]*y_] := n Log[x] + Log[y] ; 
Log [Power [x_ /; (Length [Cases [ListOf Positives , x]] > 0) , 

n_]] := n*Log[x] ; Protect [Power , Log]; 
(*There is a bug if we try to use FullSimplif y [ . .] , 

so we turn the Error Message off *) 
Off [Pattern: :"nodef"] ; 
(*Then we evaluate expression*) 
{expression, 
On [Pattern: :"nodef"] ; 

(*and finally we change rules back to original state*) 
Unprotect [Power , Log]; 

Power [Power [x_ /; (Length [Cases [ListOf Positives , x]] 

> 0) , n_.]*y_, m_] =. ; 

Power [ Power [x_ /; (Length [Cases [ListOf Positives , x]] 

> 0) , nj , mj =. ; 

Log [Power [x_ /; (Length [Cases [ListOf Positives , x]] 

> 0), n_.]*y_] =.; 

Log [Power [x_ /; (Length [Cases [ListOf Positives , x]] 

> 0), nj] =.; 
Protect [Power , Log]; 

}[[l]](*but we take the expression from before - 
in this way we don't have to introduce a new variable*) 

] 

Now it is straightforward to calculate the series above: Use AssumePos- 
itive[{A}, s + 0[x]"5] for the upper series and AssumePositive[{B}, s + 
0[x]"5 /. A — ► -B] /. B — > -A for the lower series in (|E,2.2|) . AssumePos- 
itive can take any list of variables, e.g. AssumePositive[{a, 6}, V ab 2 cd 2 ] — > 
y/a bV cd 2 . Possible extensions to this routine might include Abs[x] — ► x or 
UnitStep[x] — > 1 for x > 0. 
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